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PREFACE 


This volume contains papers which, have been accepted for 
publication by the Program Committee of the Conference on the 
Origins of Planetary Magnetism. Papers were solicited which 
address one of the following major topics: 

I. Geomagnetism in the context of planetary magnetism 

II. Lunar magnetism 

III. Dynamo theory and non-dynamo processes 

IV. Comparative planetary magnetism (terrestrial and 
outer planets) 

V. Meteoritic magnetism and the early solar magnetic 
field. 

The Program Committee consists of L. J. Srnka, Convenor (Lunar 
and Planetary Institute'), P. R. Brett (National Soienoe 
Foundation) , F. H. Busse (University of California) , 

R. B. Merrill (Lunar and Planetary Institute), N. F. Ness 
(UASA Goddard Space Flight Center), S. K. Runcorn (The 
University, Newcastle upon Tyne), and C. T. Russell 
(University of California ) . 

Logistic and administrative support for this conference has 
been provided by P. P. Jones (Administrative Assistant, Lunar 
and Planetary Institute) . This abstract volume has been 
prepared under the supervision of P. C. Robertson (Technical 
Editor, Lunar and Planetary Institute) . 

Papers are arranged alphabetically by the name of the first 
author. A subject and author index are included. The subject 
index, based on key words returned by the authors, does not 
represent the total number of papers in this volume. 

The Lunar and Planetary Institute is operated by the Universities 
Space Research Association under contract No. NSR 09-0S1-001 with 
the National Aeronautics and Space Administration. 
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MAGNETIC PROBING OF PLANETARY INTERIORS; E. R. Benton, 
Univ. of Colo., Boulder, colo. 80309 

This paper addresses the following general question: what 

can be learned about a planetary interior from measurements at 
the surface of the global planetary magnetic field? The specific 
issues considered in detail are those prompted by pressing prob- 
lems in the dynamo theory of geom agneti sm , but most are clearly 
of interest for other planets as well . Attention is restricted 
to the idealized but representative situation of a spherically 
symmetric planet containing a highly conducting liquid core 
surrounded by a nearly insulating solid mantle, the rest of 
space being vacuum. Core conductivity is presumed constant, 
while mantle conductivity, a function of radius only, is every- 
where small in comparison, falling monotonically from its largest 
value at the core-mantle boundary to zero at the planetary 
surface . 

The planetary magnetism in question then arises entirely 
from dynamo action in the liquid core, being devoid of external 
sources. This is probably too primitive a model to permit direct 
application to such objects as Jupiter or the sun, where external 
currents are also important, but the model should be relevant for 
planets such as Mercury, Venus, and perhaps Saturn, as well as 
Earth. Alternatively, if through spherical harmonic analysis, 
and/or time averaging, the effects of external sources can be 
removed from the observations, then the model is relevant to 
planets which are both internally and externally excited. 

The global surface magnetic observations required are 
spherical harmonic expansions, truncated beyond degree N {that 
is N^+2N Gauss coefficients) for the scalar magnetic potential 
as observed at two epochs separated by a time short compared to 
the ohmic free-decay time of the planetary core, yet long com- 
pared to the time required for secular variations of the magnetic 
field to penetrate through the mantle. 

For such an idealized planet, one can construct a systematic 
perturbation procedure for analyzing the magnetohydrodynamics of 
the planetary interior. The procedure, though intricate, is 
straightforward; it will be reported in detail elsewhere. In 
essence, three small parameters, proportional to core resistivity, 
mean mantle conductivity, and variations away from mean mantle 
conductivity, are introduced to describe perturbations about a 
ground state in which the core moves like a perfect conductor and 
the mantle is a uniform insulator. Information about each of the 
following properties of the planetary interior is then accessible 
and will be discussed: 

(1) Radius of the liquid core (a result of Hide) . 

(2) Depth of substantial vertical motion and intense 
electric currents in the core. 

(3) Presence or absence of unsteadiness and turbulence 
in the upper reaches of the core. 

(4) Horizontal fluid motion adjacent to the core-mantle 
boundary. 

(5) Azimuthal field strength in the upper reaches of 
the core. 
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(6) Current system in the mantle and strength of 
electromagnetic core-mantle coupling. 

This work is supported by the Earth Sciences Section of the 
National Science Foundation under NSP Grant DES75-03960. 



ON THE GROWTH RATE WITH DEPTH OF HORIZONTAL MAGNETIC 
FIELDS AT THE SURFACE OF EARTH'S LIQUID CORE; E. R. Benton and 
L. A. Muth, Univ. of Colo., Boulder, Colo. 80309 

A key parameter arising in self-consistent hydromagnetic 
(as opposed to kinematic) models of the geodynamo is the strength 
of the toroidal magnetic field developed in the liquid core 
compared to the poloidal field there. Because it is easy to un- 
derstand how differential rotation (either with depth or in the 
horizontal) would shear an initially weak poloidal field into a 
strong toroidal field, it has been conventional to assume that 
relatively strong azimuthal fields exist in the core; but their 
associated electric currents, being precluded from flowing in 
the nearly insulating mantle, implies that such fields are con- 
fined to the core, being undectable (directly) at earth's 
surface . 

On the other hand, Busse (1975) has made considerable 
progress with the convectively driven hydromagnetic dynamo 
(in cylindrical geometry) on the assumption that toroidal fields 
are no more than comparable with the poloidal fields . In this 
situation, the Lorentz force acts as a perturbation in the momen- 
tum balance, thereby permitting analytic progress to be made. 

In this paper, an effort is undertaken to determine how 
rapidly B 0 and B(|, are apt to increase in magnitude with depth 
at the outer surface of earth's liquid core. The model is based 
on a nearly insulting mantle, a nearly perfectly conducting 
inviscid core of uniform density, spherical symmetry, and gravi- 
tational driving of unaccelerated flow in the momentum balance. 
Under these conditions the radial torque balance is magneto- 
strophic (Coriolis torques counterbalance Lorentz torques) . At 
particular locations on the surface of the core it then becomes 
possible to evaluate components of the horizontal electric cur- 
rent and subsequently, to infer values for 3 B 0 / 3 r or 3B^/3r. 
Results obtained from geomagnetic field models at epochs 1965 
and 1975 will be presented and discussed. This work is supported 
by the Earth Sciences Section of the National Science Foundation 
under NSF Grant DES75-03960. 
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ANCIENT MAGNETIC FIELDS ESTIMATES FROM METEORITES : A CAVEAT AND 
HOW-TO GUIDE; Aviva Brecher, Dept, of Physics, Wellesley College, Wellesley, 
MA 02181 and Dept, of Earth and Planet. Sci. , M.I.T., Cambridge, MA 02139 

In the two decades elapsed since the exciting discovery of fossil extra- 
terrestrial magnetism in meteorites, a large amount of data has been col- 
lected and variously interpreted [e.g,, reviews 1,2,3 j. The promise and po- 
tential which meteorites hold for recovering information on ancient' magnetic 
fields [2-7] are becoming more fully realized, having received renewed im- 
petus from the re'l’ated riddle of lunar magnetism [7-9]. Since meteorites 
are extremely well-characterized with respect to equilibration, metamorphic 
or shock-reheating temperature, their radiometric ages can be determined, 
and their likely place of origin reasonably well-established, the temporal 
evolution and spatial configuration of ancient magnetic fields could also 
be conceivably reconstructed from paleointensity determinations on various 
classes of meteorites. Their origin in asteroidal, mostly undifferentiated 
objects is widely accepted, so that any magnetic fields they recorded are 
probably external, extended solar fields. These objectives can be met pro- 
vided that stringent sample-selection and data-reliability criteria are 
imposed and a suitable methodology is followed as outlined below. 

The present paper will provide an updated review of and subjective 
practical guide to paleointensity determinations on meteorites. 

(1‘) Select for magnetic study meteorites for which detai led i nformation 
on chemical-petrological classification, shock exposure and various radio- 
metric ages already exist. Ideally, data on metal phase structure, abun- 
dance, morphology and thermal history should also exist. These will aid in 
deciding on the appropriate experimental methods and in interpreting magnet- 
ic results. 

(2) Obtain large samples, typical of the whole (often inhomogeneous) 
meteorxte. It is desirable to examine mutually oriented samples from the 
same hand specimens, and several from the strewn field, in order to ap- 
praise the uniformity of magnetization. 

(3) Determine the bulk magnetic properties (e.g. , natural remanence NRM, 
magnetic susceptibility) and compare with other meteorites of the same 
class, to find out if one deals with a unique or representative meteorite. 

For example, the ureilite meteorites are extremely coherent magnetically as 
a group, whereas the achondritic breccias are very dissimilar [2]. 

(4) Analyze the NRM by standard (cold) methods (e.g., AF cleaning) to 
isolate a stable and possibly primary component. General criteria for the 
identification are: hard coercive spectrum of NRM; clustering of NRM dir- 
ections under AF cleaning, on stereonet projection; and vector component 
(Van Zijl) analysis, requiring smooth and linear convergence to the origin 
with progressive AF cleaning [10]. 

The method of NRM analysis should be appropriate for the magnetic car- 
riers present; thus one could exploit the magnetic memory effect of magne- 
tite in carbonaceous chondrites [2, 10]. Thermal cleaning, as for lunar 
samples, is counterindicated since heating incurs changes in magnetic tex- 
ture, phase-structure, and raineralogical-chemical alterations [3, 5-7]. 



5 


ANCIENT MAGNETIC FIELDS ESTIMATES FROM METEORITES 
Brecher, A. 


(5) It is desirable to saturate the sample and use the AF coercivity 
spectrum of IRM as an index of the intrinsic capability for stable rema- 
nence and as a comparison standard for any alterations due to the process 
of modelling the NRM in the laboratory. 

(6) While realizing that the physical process by which the initial NRM 
may have been imprinted cannot be duplicated adequately in the laboratory 
(e.g., sedimentation in low-g, for depositional remanence (DRM);or cooling 
rates of 0.1-100°C/m.y. with simultaneous phase change, for thermochemical 
remanence (TCRM)), select the experimental method best suited to the physi- 
cal history of the meteorite and most likely to reproduce its NRM (e.g., 
rapid post-shock cooling). This is quite subjective! In any case, one must 
ascertain that the selected method for NRM modelling does not produce sub- 
stantial irreversible magnetic alteration, or else the paleointensity esti- 
mate is unreliable. For these reasons we have adopted the Van Zijl method 
for paleointensity determination on meteorites [2, 4, 11], which requires a 
single heating to a temperature preselected to match information available 
under (1) above. Other workers prefer the double- Thellier method L5-7]; or 
other variants involving repeated heating cycles Llj3]. 

(7) The caveats of modelling experiments are: (a) preconceived ideas re- 
garding the nature of NRM (e.g., thermal) and biasing an open-minded ap- 
praisal of results. Thus, every assumption (e.g., NRM ^ a TRM; a TEM — 
and by implication, the NRM — is proportional to the ambient field at 
cooling) must be tested experimentally. Such precautions lead to surprising 
and unexpected results; for example, many lunar rocks and meteorites are 
not even capable of acquiring a TRM when cooled on laboratory time-scales 
and fields. Similarly, iron meteorites and many chondrites possess sub- 
stantial magnetic moments (TRM ) even when cooled from above the Curie point 
in field-free space [12]. Hence, even though iron meteorites can well ac- 
quire a TRM, the similarity of NRM and TRM to spontaneous magnetization 
(TRM ) precludes meaningful paleointensity estimates, (b) preconceptions in 
the interpretation of magnetic paleofields; is "primary" synonymous with 
"primordial" magnetization; is a "secondary NRM" always soft, or vice versa? 
For example, the striking fact is that soft magnetization in carbonaceous 
chondrites [2, 10] and in some chondrites [2] is often in the same direction 
as the high-coercivity NRM, so that only a single magnetizing event need 

be invoked. 

(8) Adhere to more stringent experimental tests of the reliability of 
a paleofield estimate than the often-pleaded "internal consistency"; the 
case for strong primordial fields is often excellent, until closely ex- 
amined. Thus, iron meteorites [12] can, and probably do, carry a paleo- 
remanence, and TRM vs. NRM comparisons can be obtained, with impressive 
"internal consistency," from comparing the NRM with zero-field moments. If 
similar field strengths are obtained for a group of meteorites by diverse 
experimental methods (e.g., "^l oe for carbonaceous chondrites [4-7]), these 
can be considered trustworthy, divergent interpretations notwithstanding. 
Zero-field cooling offers a promising new standard for judging the suitabi- 
lity for samples for paleointensity estimates, since it permits one to sep- 
arate internal vs. external-field contributions to the overall magnetiza- 
tion. 
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From meteorite data at hand, a plausible — but not airtight — case 
can be made that strong, extended nebular fields were present during the 
cold condensation- aggregation stages of meteorites (carbonaceous chondrites) 
and during the colllsional— shock modification of ureilites. In contrast, 
fields lower by one to two orders of magnitude are obtained from (shock)- 
metamorphosed ordinary chondrites. This may pose a temporal constraint on 
field decay over the metamorphism interval (''^10^ years) . 


References : [1] Stacey (1976) Ann. Rev, Earth and Planet. Sci, p. 

147; [2] Brecher (1977) in Comets , Asteroids , Meteorites ; Interrelations , 
Evolution and Origins (ed. A.H. Delsemme) , p. 256; [3] Guskova (1972) The 
Magnetic Properties of Meteorites (NASA-TT-F-792) ; [4] Brecher (1972) in 
On the Origin of the Solar System (Nice, CNRS), p. 260; [5] Banerjee and 
Hargraves (1972) Earth and Planet. Sci. Lett. 17 , p. 110; [6] Butler (1972) 
Earth and Planet. Sci. Lett. 17 , p. 120; [7] Lanoix et al. (1977) Proc. 
Lunar Sci. Conf . 8th, p. 689; [8] Banerjee and Mellema (1976) Nature 260 , 
p. 230; [9] Brecher (1977) Nature 266 , p. 381; [10] Brecher and Arrhenius 
(1974) J. Geophys. Res. 79 , p. 2081; [11] Brecher and Ranganayaki (1975) 
Earth and Planet. Sci. Lett. 25 , p. 57; [12] Brecher and Albright (1977) 

J. Geomag. Geoelectr. 29 , p. 379. 



ON THE OBSERVED RELATION BETWEEN MAGNETIC FIELDS AND 
ROTATION IN ASTRONOMICAL OBJECTS; A. Brecher, Department of 
Physics, Wellesley College, Wellesley, Mass., 021B1 and K. 
Brecher, Department of Physics, M.I.T., Cambridge, Mass., 02139. 

It is believed that dynamo action resulting from fluid 
motions in the cores of rotating planets and stars can give rise 
to magnetic field generation. The quantitative dynamo theories 
presented to date are mainly kinematical, aiming primarily to 
show that magnetic fields can exist in rotating fluids, though 
some dynamical arguments do lead to scaling "laws" involving 
various physical parameters. Without delving into the dynamo 
problems directly, it is our purpose here only to re-examine a 
phenomenological relation which connects the angular momenta and 
magnetic dipole moments of planets (as well as for stars and 
galaxies) which we believe should be recoverable from any 
theory of the origin of magnetism in astronomical bodies which 
claims to be correct. 

That the magnetic field of the Earth is caused by its ro- 
tation was suggested at least as early as 1891 by Sir Arthur 
Schuster and, independently, by Lord Kelvin. Schuster (1) and 
later P.M.S. Blackett (2) pointed out that a simple relation- 
ship exists for astronomical bodies between their magnetic 
moment y = aBR^ and angular momentum J = gMR^cj of the form 

y = [1] 

where M, R, B and u are the mass, radius, average surface mag- 
netic field strength and angular velocity of the object, and 
a, $ and y are constants of order unity. The presence of 
Newton's constant G in "Schuster's hypothesis" suggested to 
Blackett a possible deep new law of physics. The universality 
of such a new fundamental law was not, however, borne out by 
subsequent laboratory studies of rigidly rotating conductors or 
by observations of the radial variation of the Earth ' s magnetic 
field (3) . 

What do measurements of planetary magnetic fields reveal 
about Schuster's relation? In figure 1 we have plotted the best 
available data on magnetic moments versus angular momentum for 
several planets, the Moon and Sun. Many of the values of B are 
quite uncertain, but are probably correct to within better than 
an order of magnitude; no error bars have been assigned in the 
figure since each would be extremely subjective. The solid line 
has slope dJ/dy =1. As it can be seen, (Fig. 1) there is a 
general trend extending over ten orders of magnitude suggesting 
a linear increase of y with J. Is this result meaningful? 

Most recently, in a review of the same data, Russell (4) finds 
a similar relation, but concludes that it is not significant. 
First, he suggests, because there is no physical justification 
for the result. However, a similar argument could have been 
made about Kepler's laws of planetary motion before Newton 
showed how to derive them from an inverse square law force, or 
about Kirchoff's observation that eoy^c^ = 1 before Maxwell's 
equations were formulated. Since Schuster's relation does not 
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FIGURE 1 

obviously violate any known law of physics, we do not take the 
absence of a theoretical basis for the phenomenological relation 
to be an argument against its validity. Second, Russell ar- 
gues that plotting any two planetary extensive variables against 
one another would produce an equally appealing result. However, 
plotting, for example, the magnetic energy Ej^ versus the ro- 
tational energy E^ for each object yields no such simple scaling 
which is independent of planetary properties. (In fact, for the 
many astronomical objects which rotate near their breakup speed 
so that Gp, one finds using relation (1) that Em/Ej- varies 

roughly as (uR/c) ^ . ) Finally, it is worth noting that if one ex- 
tends the present result to stars and the Galaxy (as has also 
been done by Warwick (5)), one finds a similar result (figures 
2 and 3) extending over some 40 orders of magnitude in J and y! 



FIGURE 2 


FIGURE 3 
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Brecher, A. and Brecher, K. 

What is the significance of Schuster's relation? The 
Moon, Mercury, Neutron Stars and White Dwarfs do not fit very 
well; all these are objects in which current dynamo action is 
probably non-existent. It seems reasonable that the plot is tel- 
ling us something about objects in which dynamo action is oc- 
curing over a significant fraction of their volume, so that the 
gross macroscopic parameters considered here reflect well some 
average over their internal "microscopic" fluid motions, heat 
flows, etc. We think that the Schuster relation bears the same 
connection to dynamo theory that thermodynamic relations bear to 
statistical mechanics. We may never know anything certain about 
the detailed internal features of stars and planets and, there- 
fore, may never fully determine what must go into a detailed 
dynamo model. However, any dynamo model giving rise to 
Schuster's relation will probably be a good bet for the "correct" 
solution to the origin of cosmic magnetic fields. 

References : (1) Schuster, A. (1912), Proc. Phys. Soc. Lond., 

24, 121; (2) Blackett, P.M.S. (1947), Nature, 159, 658; (3) 
Blackett, P.M.S, (1952), Phil. Trans. Roy. Soc. London, A245, 

309; (4) Russell, C.T. (1978) , Nature, 272, 147; (5) Warwick, 

J.W. (1971), Phys. Earth Planetary Interiors, 4, 229. 



CONSTRAINTS ON THE COMPOSITIONS OF THE CORES OF THE TERRESTRIAL PLANETS 
Robin Brett, National Science Foundation, Washington, DC 20550 

Knowledge of the composition of planetary cores is important: l) to 

establish bulk compositions of planets and thus test condensation and accre- 
tionary models, 2) to understand thermal regime and history, 3) test model for 
core formation, 4} understand planetary dynamos. The existence and size of a 
core can be established by data including bulk density of the planet, coeffi- 
cient of moment of inertia, existence of a dipole magnetic field, and seismic 
observations. The composition of a core may be constrained by density consider- 
ations from seismic studies, geochemical and cosmochemical abundances, and in 
the future, by electrical conductivity measurements. 

There appears to be no relation between the size of the core and distance from 
the sun or size of the planet. The earth is the only planet for which we have 
seismic data that clearly indicate a metallic core. Density considerations 
indicate that the bulk of the core is metallic iron, and the molten outer core 
has a density 8-20 percent less than that of Fe. Cosmochemical considerations 
suggest some Ni in the inner core. Several suggestions have been made for the 
light element in the outer core; S and 0 are presently the most popular candi- 
dates. No quantitative geochemical arguments exist proving the presence of 
either element in the outer core, however if the earth accreted with appreci- 
able sulfur, it would be unlikely if the bulk of this did not enter the core. 
Potassium has also been postulated to enter the core in substantial quantities, 
thus providing a source of heat deep within the earth, but the hypothesis must 
remain speculative until the correct critical experiments have been done. 

Estimates of the compositions of cores of the terrestrial planets are closely 
related to our knowledge of condensation, accretion and coreforming processes. 

If the accretion was a heterogeneous process, v/ith the core forming first, then 
it is difficult to balance the thermal budget to get a molten core early in 
planetary history, and it is difficult to understand how appreciable quantities 
of light element entered the earth's core. Homogeneous accretion with subse- 
quent differentiation to form a core provides a ready heat source since Flasar 
and Birch have shown the temperature of the earth is raised by 2300 K by such 
a process. Such a process of core formation readily provides a light element, 
in fact the presence of a light element alloying with iron considerably lowers 
the temperature of core formation. 

Models exist (e.g., those of J. L. Lewis) to explain the compositions of 
planets including their cores in terms of condensation theory. For example. 
Mercury is proposed to have accreted from high temperature material and thus is 
low in volatiles and has a Fe-Ni core. No meteorites (apai^T from some rare 
Allende inclusions), including the refractory enstatite chondrites are without 
appreciable S and other volatiles, nor is the earth or refractory moon. This 
therefore suggests that the accretion of planetary and asteroidal bodies 
involved both high and low temperature condensates in varying proportions. The 
cosmochemical models therefore represent a very necessary beginning but may be 
somewhat unsophisticated. Geophysicists concerned with thermal history calcu- 
lations for planets should therefore not be constrained at present by the 
cosmochemists, and all core compositions and sizes consistent with the other 
constraints should be considered. 



GENERATION OP PLANETARY MAGNETISM BY CONVECTION, by 
P. H. Basse, Institute of Geophysics and Planetary Physics, UCLA, 
Los Angeles, California, 90024. 

Thermal convection is the predominant source of motions in 
the metallic hydrogen cores of Jupiter and Saturn. Convection 
driven by thermal and/or chemical buoyancy is likely to occur in 
the Earth’s core and possibly in the core of Mercury. Thus a 
dynamo driven by convection in rotating spherical shells repre- 
sents a suitable model for a general theory of planetary magnet- 
ism. While a numerical analysis of this problem is in progress, 
the experience in the non-magnetic case has shown that the 
physical properties of the problem can be best understood by 
considering the simpler geometry of a rotating annulus. The 
annulus model of planetary dynamos (Basse, 1975 j 1976) has 
recently been extended to include in higher degrees of approxi- 
mation the effects of the Lorentz force. In particular the 
reinforcement of dynamo action, the nonlinear oscillation of the 
magnetic energy about its equilibrium value, and the effect on 
the convection of the Taylor-shear generated by the deviations 
from the Taylor condition (Taylor, 1963) will be discussed. 

Because the toroidal field of the annulus model is roughly 
of the same order of magnitude as the poloidal field the model 
can be used for a comparison of the theoretically inferred field 
strength with the observed values of the planetary magnetic 
moments. This comparison is based on the hypothetical upper 
bound on the magnetic field strength for which convection in the 
form of thermal Rossby waves can exist. In addition an assump- 
tion about the maximum scale of motion derived from numerical 
solutions of kinematic dynamos in a sphere (Gubbins, 1973) 
enters into the comparison formula. This and other proposed 
scaling laws for planetary moments will be briefly discussed in 
relation to the observations. 


Basse, P. H., Geophys. J. 
Busse, P, H., Phys. Earth 
Gubbins, D. , Phil. Trans. 
Taylor, J. B., Proc. Roy. 


Roy. Astron. Soc. 437-459, 1975. 
Plan. Int. 1^, 350-358, 1976. 

Roy. Soc. London A274, 493-521, 1973. 
Soc. London A274 ,~^-283, 1963. 
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GLOBAL GEOMAGNETIC ANOMALY STRUCTURE. Joseph C. Cain, U.S. 
Geological Survey, Denver, Colorado 80225. 

The results obtained by Ptusslan and U.S. satellites have been reviewed 
for geomagnetic field components of spherical harmonic degree and order above 
el^t. Resolution is limited to scale sizes of 500 km. The spatial spectra 
of the average power and energy for each degree harmonic implies that a 
transition from core to crustal sources occurs between n = 8 and n = 12, 
Anonalies that have wavelengths in this range are not uniformly distributed 
and do not appear to be related to known tectonic elements. Their presence 
could arise either from small-scale eddies on the core surface or, in some 
instances, from a spatial modulation of secular change by lateral conductivity 
contrasts in the mantle. Anoiralies corresponding to n > 12 (wavelengths 
below about 3000 km) appear to be more related to crustal structure, although 
few have yet been investigated adequately to make detailed Interpretations, 
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ON THE POSSIBILITY OF USING LUNAR FINES TO EETERMINE THE INTENSITY 
OP THE ANCIENT MAGNETIC FIELD OF THE MOON; D. W. Collins on. Institute of 
Lunar & Planetary Sciences, University of Newcastle upon Tyne, England. 

1. Introduction 

Reliable determinations of the intensity of the ancient lunar magnetic 
field are difficult experimentally partly because of restrictions imposed by 
the lunar sample research programme. Among these are the small sample mass 
available for investigation, which sometimes leads to low or very low total 
moment for raeasin’ement, lack of samples on which to conduct concurrent tests 
for mineral alteration diiring heating, and the availability of samples whose 
magnetic properties suggest successful palaeointensity investigations. 

Bs' imparting a depositional remanence (DRM) to lunar fines it has been 
demonstrated that they contain particles which possess an NRM of hi^ stab- 
ility (Runcorn et al 1970), and it is of some interest whether this NRM can be 
utilised for palaeointensity determinations and whether any field value so 
obtained can be interpreted in a useful manner. Lunar fines have rather 
strong magnetic properties and afford the possibility of comparatively strong 
intensities of magnetization with which to work and the availability of ample 
material for idle testing of thermally-induced changes. 

2. Method 

The magnetized particles in the fines, normally randomly directed, are 
aligned by allowing them to settle throu^ a liquid into a shallow quartz dish 
in an ambient magnetic field. After drying out the fines possess a strong DRM, 
the magnetic properties of vdiich can be studied in the usual way. A Thellier- 
Thelller palaeointensity determination can then be carried out by comparing the 
loss of DRM(NRM) with the gain of PTF5M in successively hi^er temperature 
intervals . 


3. Results 

Two independent Tnellier determinations were carried out on < 1 mm fines, 
sample 7508I, from Camelot crater at the Apollo I7 site. Mineralogical 
changes were monitored by measuring low field IRM on a separate sample after 
each heating. There was evidence of changes taking place above 500°C, and 
measurements were then discontinued. Pig. 1 shows the results obtained from 
the two samples, giving apparent ancient field intensities of 0.9 Oe and 
1.2 Oe. 

^1-, Discussion 

The interpretation of the above result hinges on two main points, 

(a) the origin of the iron particles in the soil and the nature of their NRM, 
and (b) the validity of the DRM method. 

(a) It is well-known that limar soils are not simply an erosion product 
of igneous rocks and breccias, and among the Important compositional differ- 
ences is the excess iron metal (Fe ) in soils relative to that in igneous 
rocks, by a^factor of 5-10 (Housley et al 197O) . In addition, therefore, to a 
basaltic Fe contribution in soils there is other iron metal present the two 
main sources of which are meteoritic and Pe° contained in agglutinates. 
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Pig. 1 

Plots of NEM lost against 
PTRM gained for two samples 
of 75081 fines. The 
applied fields were 0 .^^ Oe 
and 0.67 Oe for runs 1 and 
2 respectively. The points 
correspond to approximately 
100 C temperature Intervals. 
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Agglutinates are the glass-welded aggregates of rock, mineral and glassy 
fragments which occur in lunar soils. They are believed to be formed by micro- 
meteorite or solar wind particle impacts into the regolith (Housley et al 1973, 
Gardiner et al 1977) and they contain fine-grained Fe° in the form of droplets 
on and inclusions in the glass. The metal is believed to be precipitated out 
of the silica-bes, either by the action of reducing gases trapped in vesicles 
at the time of impact (Housley et al 1973), by subsolldus reduction at elevated 
temperattires 700° - 900°C) in, for instance, an ejecta blanket (Pearce et al 
1972), or by a preferential sputtering process (Gardiner et al 1977)* 


Thus, the iron metal in lunar soils can occur in particles of basaltic 
rock, in glass-welded aggregates in which grainsizes spanning the superpara- 
niagnetic to multidomain range occior, and as raeteoritic iron particles. The 
likely origin of any MRM possessed by the different types of iron is difficult 
to assess, and also the contribution of each type of iron to the DRM of the 
fines. The simplest hypothesis is that the rock fragments and meteoritle iron 
possess their original NRM, and that the agglutinate iron has also acquired a 
remanence, either a C!RM or TFEM, at the time of its formation. However, in view 
of the meteoritle bombardment history of the regolith, the effect of shock on 
the NRM of all iron particles is likely to be an Important factor, potentially 
more important than may be the case with rock samples (Cisowskl et al 1973)* 

(b) The assumption that the DRM acquired by a fines sample can be used as 
the starting point for the Thellier technique appears to be a valid one 
providing certain conditions are met. The DRM should be at or near saturation, 
i.e. the particle magnetizations should be approximately parallel, and it is 
necessary that there should not be a substantial contribution to the DRM from 
isolated single-domain grains, since they cannot individually record the 
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Collinson, D.¥. 

ancient field intensity when they become -magnetized. ®ie FIRM should be 
imparted parallel to the DEM direction in the sample. 

5 . The result from 75Q81,15 

The origin of the iron metal in 75081 fines is not clear. Goldstein et al 
( 197 ^) report that the composition of the most magnetic particles from this 
soil is consistent with a minaly basaltic origin for the Pe , with ««' 20^ of the 
Pe^ content of meteoritic origin. However, an approximately equal mixture of 
meteoritic and solar wind reduced Pe° is suggested by Chou and Pearce (I 976 ). 

—■5 

The saturated DRM intensity of ihe two samples was <\/2.0 x 10 emu/g. 

This strong magnetization suggests that much of the iron metal in the fines is 
magnetized and contributing to the DRM: with 0.5^ by wei^t of Pe in the 

fines, the average intensity of magnetization of the iron is hi^, about 
0.4 emu/g, which, because of the presence of SPM iron, is a minimum value. It 
also seems likely that some of the agglutinate iron may not be contributing to 
the IBM, because the combined moment of iron and glass is too weak for good 
alignment to occur during DRM acquisition. 

The results obtained with the 7508l fines suggest that a palaeointensity 
of /-^l.O Oe has been determined, but it must be admitted that without further 
investigation the result is difficult to interpret and it is not possible to 
assign the value to any age on the Moon. The thermal and a.f . demagnetization 
curves of the DRM do not suggest more than one NRM component, and the HRM lost/ 
PTRM gained plot mi^t be expected to show evidence of two magnetic components 
of different blocking temperature spectra^ if they were present. Tiere are 
reasons for believing that 1.0 Oe may be an underestimate of the field value, 
for instance because of the possible presence of unmagnetized Pe , 

If the question of the nature and origin of the ancient lunar magnetic 
field is to be clarified, it will be necessary to make maximum use of all 
available lunar material. These preliminary results suggest that the fines may 
be of -use in elucidating the problems of the lunar magnetic field, and they 
offer the further advantage of economy of sample required (-^100 mg per deter- 
mination) and strong, easily measurable magnetizations. At the present stage 
of lunar palaeointensity research the question of age is not of overriding 
importance, since a well-defined, reliable palaeointensity at any stage in the 
Moon's history is of considerable interest. 
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ORIGIN OF LUNAR MAGNETISM BY SURFACE PROCESSES 
William D. Daily, Eyring Research Institute, Provo, UT 84601 
Palmer Dyal, NASA-Ames Research Center, Moffett Field, CA 94035 

The first lunar samples returned by Apollo astronauts ei^ht year^ ago 
were found to have a thermo^ emanent type magnetization of 10 to 10 emu/g 
acquired in fields from 10 to 1 Oe. Subsequent in situ magnetometer 
measurements reveal that much of the lunar crust is uniformly magnetized over 
scales of 10 to 100 km. Many hypotheses have been proposed to explain the 
magnetic remanence but no single theory has been entirely satisfactory. We 
review here, and schematically depict in the figure following this review, 
the principal hypotheses which do not Involve global effects but which require 
localized or near surface processes. 

(a) Shallow Fe-FeS Dynamo . Pearce et al. (1972) suggested magnetohydrody- 
namic dynamo action in pockets of molten metal just beneath the lunar surface 
to produce the fields necessary to magnetize the cooling lunar crust. Murthey 
and Banerjee^(1973) have show that if these dynamos are Fe-FeS eutectic 
liquid 10 lunar mass they may be capable of a regenerative dynamo action. 
Such dynamos may have a symmetry with the lunar spin axis capable of explain- 
ing the observed preferential magnetization alignment perpendicular to the 
spin axis. Their variety of lifetimes and sizes could also explain the four 
decade range in pallointensity values measured for samples which span nearly 
10 years in age. However, it is doubtful that these dynamos would produce 
the fields up to 1 Oe. which have magnetized some lunar samples. In addition, 
these Fe-FeS concentrations (fescons) are expected to be preferentially 
located beneath maria. Therefore, the maria should be more highly magnetized 
than highlands, and this is contrary to results of surface and orbital 
measurements . 

(b) Local Unipolar Dynamo . Enhancement of the solar wind field is possible 
by induction currents in highly conducting lava pools on the lunar surface. 

As the solar plasma and magnetic field _B are convected past the moon with 
velocity j/, a motional electric field ^ = ~ Y. ^ — drives currents through the 
lava which are closed jn the solar wind. Nagata et al. (1971) estimate 
enhanced fields of 10 Oe. at the surface for pristine solar wind conditions. 
This field may account for remanence in some samples, but it should be noted 
that fields generated by this process vary idth the solar wind field and are 
even absent at night and in the magnetotail when the plasma return path is 
not present. 

(c) Comet Impact Magnetization . Collisions of comets mth the moon may mag- 
netize the crust in a tt?o stage process proposed by Gold and Soter (1976). As 
the comet nucleus approaches, its coma engulfs the moon and the ambient solar 
wind field is compressed to about 1 Oe. between the partially ionized coma 
streaming toward the surface and the coma ions reflected from the surface. 
Coincident with the field enhancement, the comet nucleus impacts, and the 
second stage is shock remanent magnetization in the high field region. While 
this mechanism explains many properties of lunar remanence, the observed lack 
of correlation between craters and magnetic anomalies is difficult to explain. 

(d) Thermoelectric Dynamo . Generation of electrical currents and associated 
magnetic fields by the Seebeck effect in the early cooling crust has been 
suggested by Dyal et al. (1977) as a mechanism for crustal magnetism. In this 
model lava basins are electrically connected below the solid crust by magma 
and above by the solar plasma. Differential cooling results in large tempera- 
ture differences between the basins which act as junctions in a thermoelectric 
circuit and for probable electrical parameters of magma and plasma, calculated 
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fields are as large as 10 Oe. in the cooling crust between the basins. In 
this mechanism the thermoremanent magnetization and field generation processes 
are both naturally connected to cooling of large regions during crustal forma- 
tion. However, many of the model parameters such as the relative Seebeck coef- 
ficient of plasma and magma are unknown and require empirical determination. 

(e) Impact Induced Piezomagnetism . The most studied surface process for 
explaining crustal remanence is subsurface shock loading during crater forma- 
tion. Cisowski et al. (1973) have measured shock induced remanence of 10 
emu/g in lunar soil samples shocked to 50 k bar in the earth’s field. Asguming 
linear dependence of induced remanence with field, a magnetization of 10 
emu/g (comparable to that of lunar samples) may be expected for possible lunar 
fields of 10 Oe. Hypervelocity ejecta from explosive charges have also been 
used by Martelli and Newton (1977) to show that shock produced cratering in a 
basalt target results in a transient magnetic field during crater formation 
and an induced remanence in the target material. These results, as well as 
those of other experiments, are encouraging since this mechanism explains the 
crustal remanence as a natural consequence of an acknowledged lunar process. 

The mechanism also explains the observed larger remanent fields in heavily 
cratered highlands even though there is no observed correlation between indi- 
vidual craters and magnetic anomalies. 

(f) Impact Ejecta Field . Magnetic fields generated in the ionized ejecta of 
hypervelocity lunar impacts were first postulated by Hide (1972) . Srnka (1977) 
calculated ^odel dependent estimates of magnetic fields in crater ejecta of 

10 to 10 Oe. lasting for tens of seconds. Laser- target plasma experiments 
have produced measured fields up to 10 Oe. Piezo- or thermoremanence may be 
acquired by ejecta and nearby rock from these large fields. This mechanism 
explains the small present global magnetic moment as well as the larger rema- 
nent fields measured over heavily cratered highlands and smaller fields over 
sparcely cratered maria. However, there is little evidence that magnetic 
sources are correlated with individual craters. 

(g) Volcanic Ash Flow Magnetic Fields . Several schemes for generation of mag- 
netic fields by surface processes wer^ investigated by Cap (1972) . In the 
most plausible model fields up to 10 Oe. were calculated from convection 
currents in ionized volcanic ash flows. The model generates highly transient 
fields, however thermoremanent magnetization requires a relatively constant 
field as the crustal material cools from the iron Curie temperature. 
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HYDROMAGNETIC INTERACTIONS IN A FLUID 

DIELECTRIC 

Bibhas R. De 

Lunar and Planetary Institute 
Houston. Texas 77058 


The subject of magnetohydrodynamics deals with the dynamic behavior of 
a conducting fluid permeated by a magnetic field ~ and forms the basis of the 
planetary dynamo theory [ 1 ]. Here we discuss a parallel development dealing 
with the dynamic behavior of a dielectric fluid having a finite molecular pol- 
arizability, when permeated by a magnetic field. Whether or not this deve- 
lopment has implications for the possibly fluid dielectric interiors and exte- 
riors of giant planets is not clear at this time. 

The basic equations of magnetohydroelec tries ( so termed because of the 
imoortancp of the eleotrostatic energy which is negligible in magnetohydro- 
dynamics ) are first the four Maxwell’s equations ( in S.I, units ): 


V.B = 0 


0 ) 

^•X«o^ ~ ~ Pp 


(2) 

VxB = + P^oE 


(3) 

VxE = - B 


(4) 


where E is the electric field in the laboratory frame, B is the net magnetic 
field, J is the polarization current, Pp is the polarization charge density, 
and X, £3 , and p. are respectively the molecular polarizability* the permi- 
ttivity of free space, and the magnetic permeability of the fluid. Note Onat Eq. 
(2) - a modified form of Maxwell's equation - is written in terms of One pol- 
arization charge since this is tbe only charge that can arise in tJne body of the 
moving fluid, any free charges arising in a bounding conductor. The dielect- 
ric permittivity of the fluid is 

e = (1 +X)€„ (5) 

and is assumed here to be a scalar quantity. To these equations we add the 
Newton’s law 

P y == d X g (6) 

(_ P= fluid mass density, u = fluid bulk velocity ), and the equation relating 
the polarization current to the electric field [ 2 ] 

J = x«o[| + It P) 


*Referred to the bulk fluid 
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HYDROMAGNETICS OF FLUID DIELECTRICS 
B. R. DE 

The set of equations (1) - (7) completely describe the ideal magnetohydro- 
electric behavior. Note that there is no energy dissipation involved in these 
equations. 

Upon eliminating and E from Eq.(7) by means of Eqs.(3) and (4), and 
using Eq.(1), we have 

§ " Cux B) + ~ B (8) 

In this equation the two terms on the right hand side have magnitudes; 


vx ( ux B ) 


Xgp u B 
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H-e 


B 
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^ e L 


where L and t are the characteristic length and time scales of the problem. 
The ratio of the first of these terms to the second is 

X = uL^Xe^/T (9) 

If this dimensionless parameter X is much greater than unity, then Eq. (8) 
reduces to 


i = 7^ fr 00) 

Integrating once with respect to time, we have 

I = 7^ vx(uxB) + C (11) 

Where the integration constant C must be set equal to zero, since otherwise 
it indicates a monotonic change in g with time even when u = 0 everywhere. 
Under the condition X^^1> E'q.(ll) reduces to 

B = V X ( u X B ) (12) 

which IS the familiar condition for a state of "frozen flow" to exist in a mag- 
netized fluid : the magnetic flux through any contour attached to the fluid 
stays constant with time . 

If on the other hand X«1 , then Eq.(8) reduces to 

B = o/ V^B (13) 

—1 /2 

where = (H-e) is the velocity of electromagnetic waves in the dielec- 
tric. This IS a three dimensional wave equation - indicating that the distur- 
bance IS now purely electromagnetic in nature . 

We thus find that the magnetohydroelectric interactions are in principle 
very similar to the magnetohydrodynamic interactions, except that the former 
interactions are energy~conservative so that a situation with a diffusion of 
the magnetic field lines does not arise ( as does in the latter case ). 
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It has been shown elsewhere[2] that analogous to the Alfven waves of mag- 
netohydrodynamics, magnetohydroelec trie waves - consisting of oscillations 
of magnetic, electrostatic and kinetic energies — are theoretically possible. 
Such waves propagate with a velocity 


V 


c 


-1 2 2 
[ X ' + V Vc*^ ] 

a 

_-i " 6 

[ 1 + X ' + V /c 
a 


1/2 


1/2 


(14) 


1 /2 

where c is the velocity of light, and v = B/(M-p) . Note that v would 

be the Alfven velocity in the medium had the medium been conducting. The 
effect of introducing a finite conductivity in the fluid has also been discussed 
in the above reference. 

It can be shown that the induced magnetic energy, the electrostatic energy, 
and the kinetic energy can reach approximate equipartition when c^ 

«c. However, even when this condition does not prevail, magnetohydroelec- 
tric interaction may still be of interest. In the possibly superdense core 
fluid of giant planets, X may be quite large ( since X a number density of 
the molecules ), On the other hand, the velocity u of the fluid may be quite 
small and t large, making X smaller than unify. Thus the importance of 
magnetohydroelec trie effects in the planetary dynamo theory is questionable, 
but may possibly be worth bearing in mind. 
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MAGHEPIC FIELDS lU THE VICIHITY OP YENUS 
ACCORDIEG TO "VEHERA" AHD "MARINER" DATA 
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The nightside Venus cavity is the magnetosphere in which 
the magnetic energy exceeds that of the solar wind, and the 
field regularity distinguishes the magnetosphere field from the 
magnetic field within the magnetosheats. 

At distances of 2,5+6 Rv from the planet center the ni^t- 
side magnetosphere topology is characterized by stretching out 
of lines of force in radial direction; it has two lobes of mag- 
netic lines of force: "from the planet" one (northward of equa- 
tor), and "toward the planet" one (soudhward of equator), which 
are divided by a layer with the magnetic energy of deep mini- 
mum, which is quite analogous to the earth magnetosphere neut- 
ral layer 

During four of measurement seances from "Venera- 10" there 
were complete magnetosphere crossings at distances of 'viooOO* 
■*■25000 km. During other seances a boundary region was revealed 
only with During 10 seances of measurements from "Vene- 

ra-9,10" a neutral layer was observed when the ecliptic plane 
was being crossed at distances of Zs'e = 8000+23000 km. 

At closer distances from the planet the magnetotail is be- 
coming narrow. At distances of '^bRvthe magnetosphere dia- 
meter exceeds that of the geometrical umbra, and close to the 
planet the magnetotail is inside the region of the planet * s 
geometrical umbra. The magnetosphere sizes depend on conditions 
both in the magnetosheath zone and solar wind. 

Assuming a circular shape of the magnetosphere cross-sec- 
tron it has been estimated that the magnetosphere radius at a 
distance of = -34000 km is 7000 km, and at a distance of 
Xyg = -10000 km it is 'v 4500 km (Pig.1). 

The intensity of field Bx radial component in the tail is 
of about 17+15 (f’ without occurence of its diminution moving 
away from the planet which should have been expected due to a 
growth of magnetosphere sizes with the growth of distance from 
the planet. This is probably connected with a non-conservation 
of flux constancy in the venerian magnetosphere due to a perma- 
nent course of "reconnection" processes of the magnetosheath 
field with the fields of either north or south lobe of the mag- 
netosphere field. 

This "reconnection" is manifested by a diminution of the 
field Bjc radial component intensity in the tail over the eclip- 
tic with a solar-pointed field Bx in the transition zone, while 
observed under the ecliptic - with a antisolar field Bx. In 
these cases, the field Bx is decreasing to zero rather quickly, 
and then a field Bx occure directed either from the planet 
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(over the ecliptic) or ^toward the planet” (under the ecliptic). 

The best way to explain the totality of data is to assume 
that this field is the planet intrinsic field blown off to the 
night side and amplified by th^fields of surface currents gene- 
rated by solar wind £2], as it takes place in the earth magneto- 
sphere. 

The above-mentioned magnetosphere boundaries at close dis- 
tances have been determined by the field B abrupt diminution 
commencements. ^ 

The magnetosphere field characteristic feature is a small 
value of "southward” field Bg, although sometimes a weak field 
component of inverse sign has been revealed which is apparently 
due to the external field activity. 

A complicated topology of the field at distances of 
^1500*2500 km can be most clearly explained by the model of 
the planet intrinsic field interaction with the external field 
of solar origin. The intrinsic magnetosphere field is seen when 
the orbit er penetrates inside the geometric umbra up to 'Zs'a ^ 
3000 km with a = ( y/z < 4500 km. 

With greater values of Zse. some solar-origin fields are 
observed, B:*: component of which coincides by sign with the field 
Brc in the transition zone. Sometimes in such cases there is 
also observed a considerable component Bz contrail to the mag- 
netosphere field in which B 2 is small. 

At greater distances (3+5 Ru) the tail topolo©r also de- 
pends on the field intensity in the magneto sheath zone and 
retroactive. The existance of an inner magnet ot ail-like wake 
structure was supposedly suggested on the base of the spatial 
variations observed by "Mariner-10" £6 J in the distant magneto- 
sheats up 750 Rr behind Venus. According to f13J these varia- 
tions were probably due to temporal changes. 

Uearly in all measurements performed at distances of 1500^- 
2500 km the component By is directed toward (-Y^f). There has 
been discussed a probable cause of this by a deviation error 
occuring during the satellite work on the "pericentre" programme 
as well as a probable effect of currents flowing along the lines 
of force in the tail. 

Contrary to £3} it has been concluded that a probable 
effect in the day-time ionosphere of a unipolar induction mecha- 
nism from Bi component is not clearly manifested in the night- 
side tail topology. At the same time, during some seances meri- 
dian crossings small effects have been observed. 

r2 

Equation Hj = 4H0 — ^ Sin^ ^ 

Rt 

which application is limited in case of venusian magnetosphere, 
permits to infer that the tail might be formed with Ho = 10/^ 
and ^ = 31® (that is, with Gv * Gm of Mercury) as well as 
with Ho = 5^ and G = 46° which corresponds to values 
M = (1, 2+2,5) lO^^Gs-cm^ 

^ In any case, this field is not an obstacle to solar wind 
in according with the conclusions on the obstacle being of iono- 
spheric origin made on previonaly performed experiments £^4,5,6j 
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and conclusions [ 7 , 8 ], 

Initial assumptions and model, from which it had been con- 
cluded [ 9 ] that Venus had a field with the magnetic moment 
M = 6,5'10^^ GS'cm^, contradict the ”Venera-9> 10" experimental 
data, although there is some agreement between the suggested 
and observed field signs, that is, the venusian field is oppo- 
site to the earth field by sign. 

The model of Venus magnetosphere shown on Pig. 1 explains 
the observed £l0,11j peculiarities of plasma distribution. 

The observed field being very small, its origin can be 
most naturally explained by the operation of a dynamo-mechanism 
in the venusian liquid core. The Venus field appears to' be simi- 
lar to the earth field if to take into consideration the pecu- 
liarities of its rotation and precessions C12J. 
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MAGNETIC EIELDS OP PLANETS AND PECULIARITIES 
OP THEIR ROTATION 

Dolginov Sh. SL. IZMIRAN, Akademgorodok, Moscow, USSR 


It has not justified the supposition of the planetary mag- 
netic fields proportionality to the velocities of their angular 
rotation which appeared to be natural £1]. At the same time it 
has been revealed the close correlation between the changes of 
the Earth’s rotation velocity (the accelerations) and the chan- 
ges of the geomagnetic field parameters £2,3J. 

The permanent acceleration of the Earth that is associated 
to the planetary precession 

W = s: 2,25 10“^^rad/sec^ (the Poincare acceleration [A]) 

is 10^-frlO times as much as the others known accelerations of 
the ^Earth £5j6]e Naturally, that the possible role of the Poin- 
care acceleration in the planetary magnetic fields generation 
4 should be studied regardless of the problem of the preces- 
sion effiency as the source of the dynamo-process energy £7J. 
This work as well as previously published ones r8,9J is related 
to -that problemo 

The following experimental facts and physical principles 
were initial to the establishment of the analytical dependence 
H = E([cOi<SlJ) 

1 . The scaling of some parameters of the planetary magnetic 
fields f8,9J. 

2. The principle of scaling law that is read: if two phy- 
sical phenomena are similar the characteristics one of them may 
be obtained according to the others by means of the simple re- 
calculation that is similar of the measurement units transforma- 
tion [10J. 

3 . The induction equation 

^=nv^ti+'^d[V>‘ll] ( 1 ) 

dt ^ 2 

where the first term is proportional to / = R^ - (Rc and^ 
are size and magnetic viscosity in the fluid core; and the se-_ 
cond term is supposed to be proportional to k-Pp (Pp -J^CoPxSiJxr 
is the Poincare force £4J), "k" is the factor that depends on 
the relation of the Poincare force and the other ones acting in 
Gore,jO - the core density, r - radius vector). 

According to the principles mentioned above the planetary 


field is 


Ui~T a Fp_~ KJ’jC dSl 6'<n cL/ll 


( 2 ) 

where <L is angle between and<^. In comparison to the geomag- 
netic field the planetary field will beY '' 
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(3) 


S? ^ 3 £ 

Q\ flcE <^£5?e i^indE fie: 

where Rn/Rc allows to compare the dipole planetary fields 

on the surface at equator. Under the perfect similarity and 
accuracy of the initial parameters k = 1. 

There are the values of factor ”k" in the table 1 that are 
calculated according formula (3) for the investigated planets 
and the Moon, the initial data accepted for the calculation and 
references. 




Table 1 


p. 

Ho 


/- % 


Wb 

/np 


e-itf 

ym 

•‘k» 

j 

1 


170 


[ii] 

s 

(1,2+0, 9) 10“^^ 

[15] 

i 

1+2 

0,74 

8,0 

9 

1,2 

o 

5-i-IO 

N 

R 

(6, 0+2,0) 10“^ 

[16] 

1+2 

0,53 

10,0 

7 

1,3 

0 

30950 



S 

3,27 10“^ 

[17] 

23,5 

0,53 

11, 0| 

5 

1,0 

d 

60 

i 

[13] 

R 

' 5,24 10“^ 

[17] 

25,2 

0,37 

7,5 1 

4 i 

1,2 


419200 


[14] 

R 

6,66 10"^ 

[17] 

3,0 

0,70 

1,9 

7 

0,4 

3 

1 



- 

3,07 10“^ 

[18] 

1,5 

0,20 

1 

7,3 

5 

0,8 


The dispersion of the factor "k" is unsignificant ^ spite 
of the immense difference of parameters and is in th^ limits of 
the accuracy of modern knowledge of Ho» t K ando . It in- 
dicates that at least the precession effects on the mechanisoi 
of the planetary magnetic field generation and in last one the 
external gravitation forces the important part play. 

Rote that the values of "f" for Mercury and Venus determi- 
ned in £153 and fl6j were predicted in [21,93 on base of formu- 

1 3 ^ ( ^ ^ 

The fact of a small dispersion of the factor "k'» leads with 
the necessity to the equality 


where "k” - 
We are 


where \d\ - 
units, ^ - 


actor is the same for all the planets, 
apable to determine the dimension of "k" 


IKhC 


% 


\3\ 



the dimension of the electric induction in CGSM 
the permittivity. 


^ = °^^ho = "^ho 


where C - the velocity of light, \o - the velocity of the hyd- 
roaaagnetic oscillations. 


V 
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Then (4) can be expressed as 

[heJ = o' .rg/om^aeo ( 5 ) 

where 1? - a dimensionless coefficient. 

Thus the physical sense of the formula (3) is that the pla- 
netary magnetic fields are proportional to the density of energy 
flux (a flux of the Poynting's vector) carrieng by the hydromag- 
net ic waves across the fluid core boundary. 

It is necessary the new experimental data to precise the 
data of table 1. Prom formula (3) and modern models of the in- 
terntal structure of Jupiter and Saturn 19 the intensity of Sa- 
turn's magnetic field should be Hq 14000/* . 

If it occured that the polarity of Saturn's dipolar field 
is opposite of Jupiter's one then it should not consider as a 
chance phenomenon the observed alternation of the dipolar field 
signs of planets (see table 1) that is disturbed where the 
Bode's Law of the planetary distancies is disturbed. 
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DOES MARS HATE Alf ACTIVE DYNAMO? 

Sh, Sli. Dolginov, A. N* Pushkov, IZMIRAN, Moscow Region, 

Akademgorodok, USSR. 

A negative answer to the question, mentioned in the title, 
is given by Russell f1] who argue it in the following way: 

1 , Mars spins at a rate very similar to the Earth and much 
faster than Merctiry spins. It's size is intermediate between 
Mercury and the Earth. Thus, one might expect Mars to have a 
magnetic moment intermediate between that of the Earth and Mer- 
cury. However, the reported moment is much less than that of 
the Barth and that of Mercury. 

2. The magnetic field measured near Mars f2] is induced in 
its ionosphere by the perpendicular to the solar wind flow B 
component of the interplanetary magnetic field. The version 

of Mars magnetic field origin suggested by Russell ,_allows him, 
in his opinion, to eliminate the established fact that observa- 
tion magnetotail field sign is independent on the interplane- 
tary magnetic field sign, which contradicts another version of 
Mars magnetic field origin, suggested earlier by Russell E3»4]. 

Reference the reader to [" 3-6 J in which the reality of the 
martian intrinsic magnetic field is discussed, let us consider 
the validity of the new Russell's arguments, in favour of the 
induced origin of the martian magnetic field, 

1 . In accord with the unipolar induction theory T 7 ] the 
electrical field of the solar wind E = Consequ- 

ently if this mechanism is effective, the component must be 
also effective. Further, it is known that the 3-hour average 
By and components of the IMP are in 90-95% coincide in their 
sign change. The daily averaged B^-component of IMP if smell [SJ , 

It was shown in f 11j that a probedile effect in the day- 
-time Venusian ionosphere of a unipolar induction mechanism 
from Bg-component is not clearly manifested in the night side 
tail topology at the distances 1500+2500 km above the surface. 



29 


Does Mars have an active dynamo? 

Sh, Sh, Dolginov 

The same may he observed even to a greater extent in the mar- 
tian night side field. 

2. The magnetic fields of planets are not proportional to 
their const emt spin rate but they seem to depend on aiyr accele- 
rations in the spin rate. As the Poincare acceleration 
connected with the precession is the largest one of the accele- 
rations positively known on the Earth, it defines (besides the 
other parameters ) the proportionality between the pla- 
netary magnetic field intensity and the value W but 

not the value ct) f9j. 

The discovery of the magnetic field of the planet Mars 
which has the values and similar to that of the Earth, but 
the values ^ different, allowed to reveal the mentioned depen- 
dence ilOJ. 

The magnetic field of Mars -5, being completely similar to 
the terrestrial magnetic field, can not be of large intensity. 

, This probably indicates the magnetic field of 
Mars to be the field of an active <^namo. 
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SURFACE MEASUREMENTS OF LUNAR MAGNETIC FIELDS 


Palmer Dyal, NASA- Ames Research Center, Moffett Field, CA 94035 
William D. Daily, Eyring Research Institute, Provo, UT 84601 

Magnetic fields originating from two sources have been measured on the 
lunar surface; (1) remanent fields associated with magnetized crustal 
material, and (2) Induced lunar fields associated with the solar wind and 
terrestrial magnetic fields. 

Six magnetic field experiments have been deployed on the lunar surface 
by the Apollo and Lunokhod programs; one each on the Apollo 12, 14, 15, and 
Lunokhod 2 landings and two on the Apollo 16 lunar landing. The magnitudes 
of the remanent fields vary from 3 gamma at the Apollo 15 site to *v 300 
gamma at the Apollo 16 site and the time dependent induction fields vary from 
0 to 100 gamma. Measurements of these fields have been used to determine 
the vector properties and scale sizes of the remanent fields, the electrical 
conductivity and temperature profile of the lunar interior, and the global 
permeability of the moon. 

The three stationary lunar surface magnetometers were deployed by astro- 
nauts on the Apollo 12, 15, and 16 missions and measured the vector field 3.3 
times per second continuously over a period of years. A lunar portable mag- 
netometer was developed for the Apollo 14 and 16 astronauts to obtain measure- 
ments at preselected locations along their lunar traverses. Dolginov et al, 
(1977) developed a magnetometer for the automated Lunokhod 2 lunar vehicle 
which obtained vector field measurements every 5 seconds during several 
traverses in the Le Monnier Bay region. 

A local remanent magnetic field of 38 ± 3 gamma was measured by the first 
magnetometer (Apollo 12) deployed on the eastern edge of Oceanus Procellarum, 
We attributed this field to local subsurface magnetized material since analy- 
ses of Apollo 11 lunar samples first demons^rated^that crustal material had 
a natural remanent magnetism of lO” to 10“ G cm'^ g"" (Fuller, 1974). Simul- 
taneous Apollo 12 magnetometer and solar wind spectrometer measurements have 
been analyzed and show that the lunar crust at this site is uniformly mag- 
netized on a scale of 18 km (Dyal et^ al. j 1978) . Astronauts Shepard and 
Mitchell used the Apollo 14 portable magnetometer in the Fra Mauro region to 
measure fields of 103 i 5 and 43 ± 6 gammas at two sites located 1.1 km apart. 
Later in the Apollo program a steady field of 3.4 ± 2.9 gammas was measured 
near Hadley Rille by the Apollo 15 stationary magnetometer. At the Apollo 16 
landing site in the Descartes highland region the stationary magnetometer 
measured 235 ± 4 gammas and astronauts Young and Duke obtained four additional 
measurements of 112 ± 5, 113 ± 4, 189 ± 5, and 327 ± 7 gammas at locations 
along a 7 km traverse. These traverse measurements and subsequent solar wind 
plasma interaction measurements both show the local magnetic field scale size 
to be 5.7 km. The Lunokhod 2 magnetometer measured fields from 2 to 30 gammas 
in the Le Monnier Bay region (Dolginov ^ al, , 1977) . This automated lunar 
rover measured fields associated with craters and plains which varied in scale 
size from 0.5 to 7 km. 



31 


SURFACE MEASUREMENTS OF LUNAR MAGNETIC FIELDS 
Dyal, P. and Daily, W. D. 

The Apollo 12, 15, and 16 stationary and Lunokhod 2 rover magnetometers 
measured the induction of global lunar fields by time varying extralunar 
(solar or terrestrial) magnetic fields. These data were used to investigate 
the electrical conductivity of the lunar interior. Toroidal (V x B) induction 
measurements by simultaneous Apollo 12 solar wind spectrometer andgtnag- 
netometer experiments indicate an electrical conductivity of 10~° mhos/m for 
an average global crust thickness of 'v 80 km. Poloidal (eddy current) induc- 
tion measurements from different combinations of simultaneous orbital 
Explorer 35, Apollo subsatellite, and Apollo surface experiments have been 
used to determine the interior electrical conductivity profile. Both the 
time-dependent transient response technique (Dyal £t , 1976) and the 
frequency-dependent Fourier harmonic technique, (Wiskerchen and Sonett, 1977) 
yield^a conductivity profile rising from 'v 10 mhos/m at 200 km depth to 
10 mhos/m at 1000 km depth. This conductivity profile corresponds to a 
temperature profile which rises rapidly with depth from the surface to 
1100 °K at 200 km and 1800 °K at 1000 km depth for a moon composed principally 
of olivine (Duba et al. , 1974) . 

Simultaneous Apollo surface and Explorer 35 orbital field measurements 
were also used to construct a whole moon hysteresis curve, from which the 
global lunar relative magnetic permeability was determined to be 1.012 ± 

0.006 (Parkin et al., 1974). Assuming the lunar composition to be olivine 
with a ferromagnetic component of free metalic iron the total iron abundance 
was calculated to be 5.5 ± 1.2 wt %. 

In summary, the origin of the magnetized crustal regions on the moon 
remains unknown and a central problem in lunar magnetism. A high resolution 
( s 3 km)mapping of lunar remanent fields to correlate field topology with 
surface features, crustal thickness, and mineralogy would resolve the problem. 
The possible existance of a chemically or physically differentiated lunar 
core is one of the most important unknowns in studies of the thermal history 
and evolution of the Moon. Induction experiments involving two simultaneous 
lunar orbiting magnetometers and three widely dispersed surface instruments 
with 0.1 gamma resolution would be able to resolve this very important 
question. 

References 


Dolginov Sh.Sh., Eroshenko E.G., Zhuzgov L.N., Sharova V.A., Vnuchkov G.A., 
Okulessky B.A. , and Ballevsky A.T. (1977) Magnetic Field in Le Monnier Bay 
According to Data of Lunokhod 2. NASA SP-370, 433-441. 

Duba A., Heard H.C., and Schock R.N, (1974) Electrical Conductivity of 

Olivine at High Pressure and Under Controlled Oxygen Fugacity. J. Geophys . 
Res. 79, 1667-1673. 

Dyal P. , Daily W.D., and Vanyan L.L. (1978) Crustal Evolution Inferred from 
Apollo Magnetic Measurements. Proc. Lunar Sci. Conf. 9th , in press. 

Parkin C.W., Daily W.D., and Dyal P. (1974) Iron Abundance and Magnetic 
Permeability of the Moon. Proc. Lunar Sci. Conf. 5th . p. 2761-2778. 

Wiskerchen M.J. and Sonett C.P. (1977) A Lunar Metal Core? Proc, Lunar Sci . 
Conf. 8th , p. 515-534. 



LUNAR PALEOINTENSITy DETERMINATIONS - A REVIEM, M. Filler, Dept, of 
Geological Sciences, University of California, Santa Barbara, California 93106 

The determination of the Intensity of the ancient lunar fields, in which 
the natural remanent magietizatlon (NFM) of the Apollo sanples was acquired, 
has proved to be an extranely difficult task. Moreover, this task is not yet 
satisfactorily completed. In addition to such obvious difficulties as pau- 
city of sanple, lack of sample history, no record of Initial orientation of 
the sanple, inability to heat sanples without degradation, there now appear 
to be more subtle difficulties which further restrict the use of classical 
intensity determination techniques. Intensity determinations have by now 
been attenpted on both breccias and crystalline rocks. The NRM of the brec- 
cias is, however, so poorly understood, and their history potentially so 
coiiiplicated, that it is probably wise to restrict our discussion initially to 
crystalline rocks. Even the Intensity determinations for mare basalts and 
liipact melt rocks, which might be expected to have the slnplest histories, 
give inconsistent results for sub-sairples from the same rock and for rocks of 
the same age. These results suggest that the NRM is not simply a primary NRM 
of thermal origin acquired homogeneously as the rocl<s initially cooled on the 
lunar surface. Nevertheless, most workers would probably agree that fields 
of as much as one oersted are required on the lunar surface at some time to 
account for the observed NRM. Whether the NRM records the ambient field at 
the time of the formation of the rock, or that associated with some later 
event in its history is not clear. In this paper the methods of intensity 
determination are reviewed and the results discussed, with special enphasis 
on the mare basalts. 

The principal methods of intensity determination were developed for ter- 
restrial sanples, whose NRM was known to be a thermoremanent magnetization 
(TRT4). The intensity estimate is obtained from the ratio of NRM to TRM ac- 
quired in a known field, assuming the process of TRM acquisition to be linear 
in field. The method was refined by the Thelliers (1) to take advantage of 
the additivity property of pTRM, so that many Independent determinations 
could be obtained for each sample by utilizing ccmparlsons between NRM and 
TRM blocked in different tenperature Increments. This pennits one to com- 
pare the blocking tenperature distribution of NRM and TRM and can reveal 
changes in the sanple brought about by the heating. Because lunar sanples 
are severely degraded by the necessary heating to generate 'TRM, a microwave 
equivalent of the Thelliers method was developed. In this technique expo- 
sure time plays an analogous role to tenperature in the original experiment 
( 2 ). 

Single heating methods have been developed with independent checks for 
sanple degradation by Shaw (3) and by Rigotti (4). In the former method the 
demagnetization of ahhysteretic remanent magnetization is used to monitor 
the effects of heating, and In the latter ARM acquisition. Both of these 
methods can in principle be modified to correct for sanple degraxiation. like 
the Thelliers technique they use a comparison between NRM and TRM as the 
basis of the intensity estimate. 

Other methods utilize nonnallzatlon of the NRM with a remanent magneti- 
zation induced at room taiperature, so that difficulties associated with 
heating the samples can be avoided. In the simplest technique saturation 
remanent magnetization (IRM ) is used (5)* This method has the advantage 
that it is applicable to many sanples for which other data are not available. 
However, it is a rather crude method, relying only on the alternating field 
(AP) demagnetization characteristics of the NRM for interpretation of its 
nature. A more sophisticated technique is to normalize with ARM, using com- 



LUNAR PALEOINTENSITY DETERMINATIONS 


Puller, M. 

parisons in discrete increments throughout the microscopic coercivity range 
of the sanple. However, to obtain an absolute intensity estimate, from 
either this technique, or the IRMg normalization method requires a calibra- 
tion of ARM or IRMg in terms of TFM. Ohis is of course precisely the exper- 
iment, which we have all been unable to carry out, without irreversibly 
changing the sample, so it is not surprising that the various calibration 
constants are in dispute. Probably it is best to regard both of these meth- 
ods as giving rough indications of the ancient field intensity, good to per- 
haps a factor of 5, with the ARM method the more reliable of the two. 

A successful Thellier-Thellier determination was completed using 62235, 
a highland inpact rock and yielded a field value of 1.2 oe (6). No success- 
ful Thellier-Thellier intensity determinations have been conpleted on mare 
basalts. It was indeed attempts to carry out such determinations which led 
to the observation of the aberrant partial TRM acquisition (7,8). Similar 
behaviour has also been seen now in 12022 (9). A microwave Thellier-Thellier 
experiment was also attempted with this rock but a reliable Intensity deter- 
mination was not achieved. 

A convincing Shaw (3) analysis has been reported by Hoffman et al. (10) 
for the fine grain Apollo 11 basalt 10017 which gave a field value of 0.71 
oe. Shaw (3) and Rigotti (4) analyses with correction estimates have been 
carried out for 10049, a fine grain Apollo 11 basalt and gave mutually con- 
sistent results for a single sub-sanple of 0.09 and 0.1 oe (11). 

IRMg normalization has been carried out for several mare basalts with 
calibration by TRM experiments with 12053 and 10048 (12) . Estimates of 
approximately one oersted were obtained from mare basalts 10017, 10020, 12022 
and the highland crystalling rock 62235, but 10049 gave an estimate of an 
order smaller. 

ARM noimallzatlon has been used extensively by Stephenson et al. (6) and 
led them to postulate a major decrease in the field during the time span of 
the formation of the various rocks studied. It remains to be seen whether 
this analysis will stand in the light of new intensity determinations. 

Despite the apparent similarity of the Intensity estimates quoted above, 
taken as a whole, the mare basalt intensity estimates present a picture of 
remarkable variability. For example, for almost all of the results quoted 
above, other determinations from different sub-samples of the same sample 
have given different results. In fact, individual sub-samples generally 
give widely divergent results, as do individual rocks of conparable age, f2?om 
similar localities on the moon. It therefore seems unwise to assume that the 
NRM is simply a primary TRM acquired at the time of formation of the rock. 
Indeed this concern is reinforced by the variability of the AF demagnetiza- 
tion characteristics of sub-sanples from a single rock. The mare basalts, 
which give the hipest intensity determinations tend to have a fine grain 
size. However, some vitrophyres give low estimates, so that the rapid ini- 
tial cooling inplied by fine grain size does not always correlate with high 
intensity estimates. Anomalous pTRM acquisition behaviour is exhibited by 
several of the mare basalts, which give the high intensity estimates. These 
samples exhibit a remarkable increase in their capability to carry remanence 
and in their magnetic hardness upon heating. This makes one wonder if some 
of the variability in the NRM could be due to inhomogeneous heating of the 
samples on the moon. One event in the sample’s history which might give such 
inhomogeneous heating is Impact related shock. Until the variability 
of the NRM from individual sanples is better understood, one can only con- 
clude that the fields required to account for NRM appear to be variable and 
to include fields of the order of magnitude of an oersted. When the magneti- 
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zation was acquired is not yet clear. Neither is it clear what the nature 
of the fields inay have been. 
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VEtTUS: a) MAGNETIC FIELDS, and b) lONOPAUSE ALTITUDE 

B. E. Goldstein and R. S. Wolff, Jet Propulsion Laboratory, 
Pasadena, California 91103. 

Does Venus have a permanent magnetic field? Russell (1976) 
interpreted the increase in magnetic field measured by Venera-4 
with decreasing altitude as due to a planetary field; Russell and 
Goldstein (1976) interpreted the reversal in the component 
of the Mariner V field upon entering and leaving the magnetotail 
as also indicating a planetary field. However, the reversals can 
also be explained by assuming a prior change in the solar wind 
transverse magnetic field and that ionospheric field had not 
adjusted to the more recent solar wind draping direction. More 
recently, Dolginov et al. (1977) have reported on Venera-9 and 
-10 data in the wake of Venus. The magnetic field observations 
as a whole indicate a B^ component directed away from the sun in 
the Northern hemisphere and a solar directed component in the 
Southern hemisphere, in agreement with the previous estimate of 
direction. However, Eroshenko (in Dolginov et al.) shows that in 
the near wake of Venus the B component is directed as would be 
expected from draping of the transverse component of the solar 
wind field. 

These observations pose several interesting problems. On 
theoretical grounds, the induced magnetic fields are due to 
currents that close through the dayside ionosphere (Daniell and 
Cloutier, 1977). Field lines passing into the nightside tail 
must pass above or through the ionosphere at the terminators and 
close through the dayside ionosphere. However, the observed flux 
in the tail seems to be considerably more than can pass through 
the dayside ionosphere; which has a much smaller cross-sectional 
area than the tail. If this is in fact a real problem, we must 
conclude that currents are also flowing through the nightside 
^ ionosphere and allowing flux to pass through the equator beneath 
the ionosphere. The conductivity of the nightside ionosphere is 
at least a factor of ten less than the dayside due to the lower 
ion densities. Additionally, in the model of Daniell and Cloutier 
the solar wind makes effective electrical contact with the iono- 
sphere only at low altitudes near the subsolar point; electrical 
contact near the terminators is inhibited by poor conductivity 
across field lines at high inopause altitudes. Alternatively, 
we suggest that field lines crossing the equatorial plane beneath 
the ionopause may be connected electrically to the solar wind at 
higher latitudes (an open ionospheric model rather than a closed 
bag-like field model) . Such a model would allow a significant 
portion of the current to close through the nightside ionosphere. 
However, even in this case it is clear that most of the flux will 
pass through or above the terminator ionosphere. 

In the opposite extreme, we note that thermospheric winds 
will induce large currents if a permanent magnetic field is 
present; a dipole moment aligned parallel to the solar wind 
direction would be completely confined beneath the ionosphere by 
a toroidal current system. A perpendicular dipole can not be 
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completely confined; instead, currents will flow over the poles 
due to VxB drift and return at lower latitudes due to the 
consequent electrostatic field. The magnetic flux will be 
pulled to the nightside and emerge from the ionosphere in a 
large "spot” where convection balances diffusion. 

From these considerations, we believe that directly behind 
Venus and near the planet induced magnetic fields will dominate 
near the terminator and a permanent magnetic field should be 
better observable closer to the anti-solar point. However, 
Eroshenko notes that the best evidence for dominance of induced 
magnetic fields is in the region where we believe that changes 
for detecting a permanent field should be favorable. The 
observational and theoretical uncertainties in interpreting 
the Venus magnetic field data make it difficult to reach a 
conclusion as to the presence or absence of a permanent magnetic 
field at the present time. 

Additionally, we have analyzed Radio Frequency Occulation 
(RFO) estimates of ionopause altitudes (Fjelbo et al., 1970; 
Ivanov- Cholodney et al., 1978) by taking into account upstream 
solar wind dynamic pressure, the dependence of solar wind pressure 
upon the zenith angle, and the variation of ionospheric pressure 
with the solar ultraviolet flux. The atmosphere pressure as a 
function of altitude obtained in this fashion is somewhat uncer- 
tain. However, a reasonable estimate of the pressure scale height 
between altitudes of 260 km to 320 km is about 23+5 km; this 
figure is in good agreement with density scale heists obtained 
from direct RFO measurements of density. In addition, the scale 
height between 300 km to 450 km appears to be substantially 
larger, with a value of 100 + 40 km. The increased value of the 
scale height depends upon two ionopause observations at zenith 
angles of 74®. The observation suggests a higher temperature in 
the upper portions of the Venus ionosphere; regions at high 
altitudes in some cases may be on magnetic field lines that 
connect to the solar v?ind so that cooling by electron interchange 
with the lower ionosphere would be inhibited. Alternatively, 
other interpretations such as a decreasing molecular weight or an 
effect of a solar wind boundary layer flow are possible. 
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ORIGIN OF MAGNETIZATION IN LUNAR BRECCIAS, Wulf A. Gose, Geophysics 
Lab., Univ. of Texas Marine Sci. Inst., Galveston, TX 77550 

The origin of magnetization in lunar breccias has been the subject of 
controversy. Part of the difficulty in interpreting magnetic data on 
breccias stems from the fact that there are no suitable terrestrial analogs 
with which to compare them. Secondly, none of the returned samples is an 
oriented bedrock sample and thus no information can be obtained from the di- 
rection of magnetization. Thirdly, the breccias contain at least two, mostly 
more generations of metallic iron, the dominant magnetic mineral. 

The most suitable samples for magnetic studies are the samples from the 
boulders at the Apollo 17 landing site. These samples are unique in that a 
series of handspecimens were collected from one (or a few related) boulder 
and for some of these handspecimens it was possible to obtain relative orien- 
tations. We have reported the results on 26 samples from seven hand speci- 
mens from the station 6 boulder (Gose aL . , 1978). Based on these data we 
propose that the natural remanent magnetization in these breccias is the vec- 
tor sura of two magnetizations, a pre-impact magnetization and partial thermo- 
remanence acquired during breccia formation. The relative contribution of 
the two components is controlled by the thermal history of the ejecta, which 
in turn is determined by its clast population (abundance and size). Depend- 
ing on the clast population, the NElf can be a total thermoremanence , a par- 
tial thermoremanence plus a pre-impact magnetization, or a pre— impact magne- 
tization. 

This model of thermal overprinting seems to be applicable to all lunar 
breccias of medium and higher metamorphic grade. In the low metamorphic 
grade breccias the magnetization might well be controlled by shock effects as 
argued by Cisowski ^ (1974). A proper understanding of the origin of 

magnetization in lunar breccias is important in that the breccias provide the 
possibility of obtaining information outside the time range of mare basalt 
crystallization ages. In addition, they are the prime candidates as a source 
rock for the magnetic anomalies measured with subsatellites. If the above 
model is correct then it should be possible to map ejecta blankets or rays 
from orbit thereby significantly contributing to our understanding of lunar 
geology 
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:EiIET?'6-Y S’OURCiS FOR -THE EARTH'S MAGNETIC AN.D. 

consequehi&IsfF^^'The ojher, planets 

David (JubMns, Dept. Geodesy & Geophysics, Cambridge, 
England 


We know more about the Earth's magnetic field than we will 
ever know about that of the other planets, and geomagnetism is 
a showpiece for studies of magnetism of the terrestrial planets. 
The poor electrical conductivity of liquid iron means that a 
constant energy supply is needed to maintain the field, the 
rate of supply being quite close to the observed heat fltix from 
the interior. There are several different possible sources of 
energy: radioactive heating, tidal dissipation and th^.^nergy 

release of earthquakes are all just about adequate. Numerical 
estimates of the energy required, that allowed, and that 
available are all close to one another. If we lived on another 
planet, it would be quite easy to deduce that Earth had no 
dynamo driven magnetic field at all. The terrestrial planets, 
plus some sattelites like the Moon and Gannymede, are not that 
dissimilar and the calculations for the Earth serve as a warning 
against firm conclusions being drawn for other planets. Jupiter 
and the outer planets are quite different: they have no energy 

shortage and in many respects resemble the sun more than the 
Earth. The central problem here is one of explaining the 
temporal and spatial variations of the magnetic field rather 
than its very existence. This discussion will be restricted 
to the Earth, and comparisons with the terrestrial planets. 

Dynamo theory assumes magnetic fields to be generated by 
fluid motions. In steady state the work done by the fluid 
against magnetic forces is equal to the Ohmic heating, which 
is one estimate of the energy requirement. It can be computed 
from the electrical conductivity (5X10^ mho.m”^) and the 
magnetic field. The latter is only known at the surface of 
the core, and the values inside must be guessed. Comparing 
numerical values based on lower boimds and idealised dynamo 
models, the Ohmic heating for the Earth's core must be between 
Q XP 

10 and 10 watt, depending on the spatial length scale of 
the total field and on how large the toroidal field is. 

A power source is required to drive the fluid motions, not 
all of which goes to Ohmic heating. Viscosity and diffusion 
of different components of the liquid mixture also dissipate 
heat, but these effects are probably small. A lot of power 
is lost without ever being converted to kinetic energy of 
the fluid motions. Two adequate energy sources are the slow 
loss of the Earth's ro •Rational energy, and seismic energy. The 
mechanisms available can only drive oscillatory motions directly 
which are ineffective for magnetic field generation, except for 
turbulent effects which would probably involve much larger 
Ohmic heating anyway. Thermal convection can be driven by 
■radioactive heat sources or by gradual cooling of the core. 

Heat is' converted to electrical energy with an ideal thermo- 
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-dynamic efficiency of only 209^* More heat is conducted down 
the adiabat leaving only a 59^ conversion. Heat mi^t generate 
a modest magnetic field without exceeding the observed heat 
flux, but the numerical values are uncomfortably close. 

Gravitational power has received a lot of attention recently. 
Cooling will shrink the core which involves a loss of gravit- 
-ational energy. This is particularly large where core liquid 
free 2 ses because of the density contrast of about D5Mg.m-3 
between solid and liquid iron at core pressures. Some work 
is done in deformation and some appears as heat, helping to 
fuel a thermal dynamo. The large density jump at the inner 
core boundary indicates that there is a compositional difference 
between the inner and outer cores. In this case freezing of 
the liquid will lead to freezing out of a heavy, iron-rich 
component, leaving a light fraction in the outer core. The 
light material rises up, driving convection directly. Some 
light material may leak away by molecular diffusion, but most 
of the available energy is dissipated in magnetic fields. 
Hiindreds of Gauss can be sustained while keeping within the 
Earth’s thermal budget, and this makes gravitational povi^er 
an attractive contender for the Earth's dynamo. 

A similar idea is that the core originally formed near the 
melting point and solid particles have been settling out ever 
since, stirring the liquid. This problem is bound up with that 
of core formation. Further complications like that of phase 
changes have been studied lead to interesting dynamical 
consequences. But the controlling factors are the heat loss 
through the mantle and the initial thermal state of the core, 
which v/ill determine the cooling and differentiation rates. 
Comparisons with other planets must incorporate these tw'o 
features accurately. The other terrestrial planets all have 
small magnetic fields that could be due to magnetised rocks. 

Only the demonstration of a liquid core or the discovery of 
secular variations would prove a dynamo driven field. The 
cooling of the planetary core is determined by the complete 
energy budget discussed above, and the uncertainties are too 
great to make sound deductions about whether particular planet- 
-ary cores are molten or not. 



A mw liBTHClD FOR DEO!Ea?iaWl«& THE RADITJS OF THE ELSCTRICAIiY- 
CONDUCTIWCJ- FLUID CORE OP A PIANET PROM EXTERHAL MAODETIC 0BSERVATI0N8: 

RECENT DEVEIOPMBNTS. 

Raymond Hide, G-eophysical Fluid Dynamics Laboratory, Meteorological Qfrice, 
BracknelJ., Hericshire, U.K. 

The new method Mature 64.O-6AI ; exploits the 

applicability of Alfven’ 3 "frozen field'* concept when dealing with 
changes in the planetary magnetic field over intervals that are much 
less than the time-scale of ohmic decay. The outer radius of the 
electrically-conducting fluid core is determined from magnetic 
observations near the surface of the planet r=rg by finding from those 
Observations that sxuTface r=ro within the planet at which the time 
variation of the total number of lines of magnetic force K^rJ ^,equal 
to times the average magnitude of the vertical component of the 

magnetic field) vanishes. '\Yhen applied to the Eeirth (Hide, R., loc. cit.; 
Hide, R. , Lowes, P.J. and Malin, S.R. C,, in preparation), for which 
rg = 6371 km, the method gives for a value which diners insignificantly, 
by less than 2fo, from the seismological value of 3486 ? 5 km. This 
determination is based on one of the best available models of the 
geomagnetic secular variation^ where determinations of the main geomagnetic 
field (namely that field obtained when contributions due to ionospheric 
and magneto spheric currents have been eliminated by taking annual mean 
values of the observations) for the period 1965 to 1975 were fitted by 
a spherical harmonic series up to and including terms of degree 8. The 
apparent success of the method ^rfien thus tested on terrestrial data 
indicates that it might prove valuable in the investigation of the 
internal structure of magnetic planets with unknown values of r^., such 
as Jupiter, Saturn and Mercury, provided that those responsible for 
planning .future missions with space probes ensure that sufficiently 
detailed measurements are made of the spatial configurations of the 
magnetic fields of these planets and of secular changes. 
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HOW STRONG IS THE MAGNETIC FIELD IN THE EARTH’S LIQUID CORE? 

R. Hide, Geophysical Fluid Dynamics Laboratory, Meteorological Office, Brack- 
nell, Berkshire, England and P. H. Roberts, School of Mathematics, University 
of Newcastle upon Tyne, Newcastle upon Tyne, England. 

The geomagnetic field in its region of origin, the Earth's liquid metal- 
lic core, consists of a poloidal part B„, with lines of force that emerge 
from the core and penetrate to the Earth's surface and beyond, as well as a 
toroidal part B , of unknown strength, with lines of force that have no 
radial component and are thus confined to the core . Toroidal fields are 
readily produced by inductive interactions between poloidal fields and hori- 
zontal fluid motions and it is likely but by no means certain that B greatly 
e^ceed^ B in magnitude: since the volume integral over the whole core of 
(V X (B^ ^ Sp))2/y2a is a lower limit to the power requirement for the geo- 
magnetic dynamo (where a is the electrical conductivity and p the magnetic 
permeability) any discussion of the energetics of the dynamo process must in- 
volve inter alia the estimation of <B^>, the average magnitude of B^. 

We review here various methods that have been used to estimate > and 
we propose a new method, which exploits certain dynamical considerations of 
hydromagnetic flow in a rotating fluid. The determination of the electric 
field at the surface of the Earth that would be associated with B in the core 
is in principle the only direct method, but the electric field is^so weak as 
to be undetectable. The earliest "eddy" models of the geomagnetic secular 
variation require ~ 3001 , whilst the interpretation of the geomagnetic 

secular variation, including the westward drift, in terms of hydromagnetic 
oscillations (strongly modified by rotation) of the core requires > - lOOF. 
The various kinematic dynamo models that have been applied to the Earth give 
a range of values of <B^>, from the 5F associated with "a^-dynamos" to 1750r 
associated with "aw-dynamos" . The new method explored here, exploits deter- 
minations of features of the flow field in the upper reaches of the core on 
the basis of the "frozen field" hypothesis . When combined with dynamical 
considerations of the radial component of vorticity an expression can be ob- 
tained for the vertical gradient of B in terms of the flow field and B . 
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BEHAVIOR OF THE GEODYNAMO DURING THE MATUYAMA-BRUNHES REVERSAL, 

K.A. Hoffman, Physics Department, California Polytechnic State University, 

San Luis Obispo, CA 93407. 

Paleomagnetic records of polarity transitions are seen to display a 
dependence on the site of observation when paths of the virtual geomagnetic 
pole (VGP) are compared (1,2). For example, nearly all available reverse- 
to-normal (R-^) transition paths from sites in the northern hemisphere are 
found to reside on the hemisphere centered about the site meridian, the so- 
called near side. In particular, data corresponding to the Matuyama-Brunhes 
reversal, for which there presently exist records from five site localities 
scattered about the northern hemisphere, well illustrate this claim when all 
paths are plotted with respect to a common site longitude (Fig. 1). The 
absolute path locations along with the associated sites (SI through S5) are 
shown in figure 2. 

The behavior depicted in figure 1 indicates that the field configu- 
ration during this R^N transition was controlled by low-order zonal com- 
ponent(s). Moreover, such a situation may be roughly modelled by considering 
the reversal process to start either at low latitudes or in the southern 
hemisphere of the core and subsequently extend axi symmetrically until the 
entire core has been affected (1,2). 

These zonal "flooding" models predict for the Matuyama-Brunhes (or, for 
that matter, any R-^ transition) VGP paths which run along the site meridian 
for the case when the site is located in the northern hemisphere. They 
cannot explain, however, why most paths (although on the near side) are not 
found in close proximity to the site meridian (see Fig. 2). 

Site 



Figure 1: VGP paths corresponding to seven Matuyama-Brunhes 
reversal records from five northern hemisphere 
site localities (see Fig. 2), and plotted with 
respect to a common site longitude. The "near 
side" lies between -90° and +90° meridians. 
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In order to account for this non -axi symmetric behavior, the flooding 
approach has been extended so as to incorporate both azimuthal as well as 
latitudinal proliferation of magnetic field of opposite sense throughout the 
core from a point of origin. That polarity transitions may initiate locally 
in the core has been previously suggested (3). 

For the present model cylindrical geometry of the magnetic source region 
has been chosen for computational simplicity. Starting from a point on the 
periphery of the cylinder, flooding is assumed to occur symmetrically in the 
azimuthal +4> and -({) as well as along the +z and -z directions. The five 
variables associated with the model include the dimensions of the source 
region, the coordinates of the point at which the reversal originates, and 
the ratio of the flooding velocities. By constraining the transition to 
start at the equator of the core and applying two other reasonable assump- 
tions reduces the number of fitting parameters to two; namely r, representing 
the apparent equatorial depth of the magnetic sources, and (}»p, the longitude 
at which the reversal originates. 

The predicted paths associated with the site localities SI through S5 
are shown in figure 3 for the case when r=0.13 Tq (rg=radius of the earth) 
and i|)n=338°E. Comparison of figure 2 with figure 3 reveals a surprisingly 
good fit of the observed Matuyama-Brunhes reversal data. Taking as a quan- 
titative measure of the goodness of fit for each case the absolute value of 
the angular distance between the observed and predicted equatorial crossings, 
one finds a range from 0*^ to 39° with an average of only 23°. 

Validity of the modelled solution also depends on its ability to accu- 
rately reflect the observed variation in field intensity. In this regard, 
only one acceptable intensity record (from site S3 (6)) is presently avail- 
able and the predicted intensity variation for this site is found to be in 
close agreement with that observed. 

Hence, these results show that the non-axisynmetric flooding approach 
to the geomagnetic reversal process provides a close simulation of observed 
world-wide transitional field behavior, at least for the presently available 
records of the Matuyama-Brunhes. 
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COMBINING PALEOINTENSITY METHODS: A DUAL-VALUED DETERMINATION ON 
SAMPLE 10017,135, K.A. Hoffman* J.R. Baker and S.K. Banerjee, Dept, of Geol. 
and Geophys., Univ. of Minnesota, Minneapolis, MN 55455. 

The Shaw method (1) for determining paleointensity involves the 
comparison of demagnetization characteristics of a sample's natural remanent 
magnetization (NRM) with that of a laboratory induced thermoremanent 
magnetization (TRM). In contrast, the method developed by Stephenson and 
Collinson (2) relies on a similar comparison of NRM with an anhysteretic 
remanence (ARM). A single heating of the sample to above its Curie point 
is required, however, by both methods since for the latter case the ratio 
of TRM to ARM efficiency f must be known. Accounting for the alteration 
of the magnetic carriers during heating is accomplished for the Shaw method 
through the comparison of ARM's induced prior to and after heating, however, 
no such correction procedure has been proposed for the "ARM-method". 

We suggest that the NRM versus TRM plot necessary for the Shaw method 
determination may be used to make the needed correction to f‘, provided that 
an unaltered coercivity range is recognizable. Hence, if this condition is 
satisfied, two primarily independent paleointensity values may be determined 
for a particular sample, one serving as a check of the other. 

Data corresponding to sample 10017,135 are shown here for illustration. 
The NRM for this fine-grained vesicular basalt sample is found to be stable 
and well-behaved in direction during AF-demagnetization similar to the 
behavior previously reported for other chips (3,4). The plot of NRM demag- 
netization versus that of a 0.5 oe TRM is shown in figure 1. The straight 
line drawn through the data corresponding to AF steps between 0 and 500 oe 
strongly suggests at most a minimal amount of alteration to this coercivity 
range occurred during heating. 
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Experiments conducted on synthetic multi -domain Fe indicate that ARM is 
invariably more sensitive to alteration than is TRM. The plot of an 0.5 oe 
ARM induced in 10017,135 before and after heating (Fig. 1), however, indicates 
behavior which closely resembles the NRM-TRM plot thus confirming the lack of 
perceptible alteration to the lower coercivity range upon heating. The Shaw- 
type determination renders a paleointensity of 0.71 oe, similar to a previously 
reported value (4). 

The ARM-method determination may be made by noting that an extrapolation 
of the linear segment to the TRM intercept (Fig. 1) makes possible the 
recovery of the magnitude of TRM which would have been acquired by the sample 
if nc alteration took place. (Indeed, the NRM-TRM plot indicates that the 
principal effect of heating this sample was to produce highly stable carriers. 
Hence, the total TRM value, if uncorrected, would result in an inflated value 
of f.) With the value of ARM acquired in an infinite alternating field, 
determined from a plot of 0.5 oe ARM acquisition versus 1/H, f is found to be 
equal to 0.94. Theoretically, this value should be no less than unity (5); 
however, a study of magnetite-bearing samples (6) casts doubt on this predic- 
tion for the case of multi-domain grains. 

Figure 2 shows ARM acquisition versus NRM lost for the indicated AF’s 
for sample 10017,135. The plot is remarkably linear and its slope coupled 
with the determined value of f‘ render a paleointensity of 0.93 oe. Thus, the 
average paleointensity corresponding to this 3. 6-3. 8 AE lunar basalt from 
these two primarily independent determinations is 0.82 oe, a value which lends 
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support to the earlier suggestions of relatively strong paleofields experienced 
on the lunar surface. 
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APOLLO SUBSATELLITE MAGNETOMETER DATA: CONSTRAINTS ON THE 
INTERPRETATION OF LUNAR MAGNETISl^ L. L. Hood, Institute of Geophysics 
and Planetary Physics, University of California, Los Angeles, California 
90024. 

The major result pertaining to lunar fossil magnetism of Apollo orbital 
magnetometer experiments is the detection of widespread magnetic anomalies 
(Coleman et al . , 1972) whose invariance on successive, closely spaced orbits 
and with externally applied field strength and orientation suggests crustal 
sources of coherent magnetization. Additional results established to date 
include: ^ 2 

(a) An upper limit for the global permanent dipole moment of '^0.2 y— & (lY = 
10“5 Gauss; 1 Rj^ = 1 lunar radius = 1738 km) (Russell et al ., 1975) . Thus 
although coherent remanent magnetization must exist on a local scale to 
explain measured surface anomalies ranging from 3 to 300y in magnitude, 

a uni-directional magnetization of anomalous source regions (or any 
significant portion of the crust) is ruled out. 

(b) No general correlation of orbital anomaly maxima with prominent surface 
morphology such as the larger craters (Coleman et al ♦ , 1977b; Hood et 

al. , 1978b). By default, compositional contrasts rather than topographic 
variations are the indicated causes of most observed anomalies. Several 
lines of evidence suggest that concentrated deposits of crater and basin 
impact-generated ejecta materials may be expected to possess a stronger 
than usual remanent magnetization and could explain the majority of 
anomalies detected from orbit (Strangway et al . , 1973) . The linear mag- 
netic anomaly apparently associated with the rille Rima Sirsalis 
( Anderson et al . , 1977) is difficult to explain on this basis, however. 

(c) A depletion in number and amplitude of anomalies detected across the 
front side maria relative to those detected across the far side highlands 
(e.g. Coleman et al . , 1977a). This is consistent with returned sample 
studies which show that mare basalts typically possess a lower stable 
remanence (''^ 10~6 e.m.u./g) than some classes of breccias ('v> 10 “^ e.m.u./g). 

(d) A correlation of remanent magnetism and crustal geology on the central 
far side (Russell et al . , 1977) . In particular, a strong dominance of 
Imbrian aged crater melt and ejecta debris is present in the Van de 
Graaff-Aitken region roughly coincident with the occurrence of large 
magnetic anomalies. One interpretation is that the ancient magnetizing 
field was significantly stronger during the Imbrian period than during 
earlier epochs. An alternative interpretation is that highly susceptible 
ejecta was preferentially deposited at that location (nearly antipodal 

to the present Imbrium basin) at the time of the Imbrium Impact event 
(Hood et al . , 1978b) . The second possibility is rendered more plausible 
by the existence of a second region of strong magnetic anomalies visible 
on the high resolution electron reflectance maps (frontispiece, Proc . 
Eighth Lunar Sci. Conf . , Vol. 1) and located nearly antipodal to the 
Orientals basin. Further, some lunar geologists (e.g. Stuart-Alexander , 
1976) have independently noted the antipodal deposition theory in dis- 
cussing the origin of peculiar grooved terrain found in both areas as 
well as an area antipodal to the Caloris basin on Mercury. "Gardening" 
of older magnetized deposits by subsequent bombardments could also have 
contributed to the relative magnetic strength of younger Imbrian and 
Orientale ejecta materials. 
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(e) Inferred directions of magnetization for major source regions across the 
central far side are generally highly inclined with respect to one 
another but seldom possess a strong north-south component ( Hood et al , , 
1978a, b). The latter characteristic is difficult to explain in terms of 
entirely local, non-global causes of the magnetizing field and implies 
that any assumed internal magnetizing dipole or external linear magne- 
tizing field must have been oriented mainly perpendicular to the present 
spin axis. The lack of alignment for nearby Inferred magnetization vec- 
tors suggests that reorientation with respect to the Moon of the magne- 
tizing field may have tjrpically occurred on time scales less than the 
mean time betxjeen the formation of nearby source regions. This charac- 
teristic may also represent an explanation for the negligible global 
dipole moment (Hood et al . , 1978a) . 

Only a fraction of usable magnetometer data (that collected in the geo- 
magnetic tail lobes) has been mapped. A large amount of excellent unmapped 
data exists for times when the Moon was in the solar wind and the subsatellite 
was at a relatively low altitude in the lunar wake. We are in the process of 
mapping all usable data using methods recently developed for mapping the 
geomagnetic anomaly field. An exact accounting for variations of the sub- 
satellite altitude is one feature of the new mapping procedure. Figure 1 
shows a prototype map constructed from v30 Apollo 15 lunar wake orbit seg- 
ments. The large anomaly is evidently due to a relatively isolated source 
region whose bulk magnetization vector is nearly radially inward. Its ampli- 
tude {'u2.7y at 88 km) is comparable in magnitude to the Van de Graaff-Aitken 
anomalies and its location (235° E. Long., -23° N. Lat.) is approximately 
antipodal to that of the Crisium basin. 
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A Tldal-Precessional Dynamo for the Terrestrial Planets H.HOTJBEN , 
Laboratory for Planetary Studies , Cornell U .- Recent developments in the ex- 
ploration of the terrestrial planets including the discovery of Mercury’s 
intrinsic magnetic field (Ness et al., 1975) and the possibility of an ancient 
lunar magnetic dynamo (Runcorn, 1975) require the reexamination of theories of 
planetary magnetism, including the geodynamo. In this light it is interesting 
to note that a d 3 niamo driven by the combined action of diurnal tides and pre- 
cess ion- similar to that proposed by Boudi and Lyttleton (1953) and Mallaus 
(1968) -can be shown to be viable for the terrestrial planets and to predict 
correctly the magnetic moments of Mercury, the ancient Moon and (with certain 
provisos) the Earth. 

The magnetic dynamo equation is converted into two equations with very 
different timescales: l)an equation for the diurnally varying magnetic field 
(for which the frozen field approximation is valid) and 2)an equation for the 
slowly varying average field (with time variations on the order of the free 
decay time of the field) in which the VXH term is the time-average product of 
the diurnally varying fields. The problem is solved for a spheroid of small 
elllpticity in conjunction with the fluid dynamic momentum equation for an in- 
compressible fluid. A large toroidal field-satisfying the boundary conditions 
for a conducting sphere surrounded by an insulating mantle-is assumed to exist. 
It is shown that such a field can be sustained by dynamo action in the sphere. 

A steady solution is found and the magnitudes of resultant toroidal and polold- 
al field components are calculated. 

The following are important intermediate results. 

1) If the toroidal fields have the correct radial dependence to satisfy 
the boundary conditlons-sphericalBessel functions (Backus, 1958) -the Lorentz 
force cannot be balanced by pressure gradients and a differential rotation is 
driven (cf. Busse, 1975). When a pololdal field is present, this differential 
rotation (acts to regenerate the toroidal fields (Gubbins, 1974). The energy 
for this process comes from the planetary rotation, as the spinup time to uni- 
form rotation is short in conducting coves. 

2) The diurnal tides and precession are nearly resonant with the spher- 
oid's toroidal Spinover mode (Greenspan, 1968). The strange motion of this 
mode does not lead to dynamo action by itself, because it has no radial com- 
ponent (Rochester et al. , 1975). However, the motion does lead to a diurnally 
varying magnetic field. As a result, there are diurnally varying Lorentz 
forces in the conducting core. These Lorentz forces, in turn, can drive 
motions which have radial components and which are candidates for dynamo action. 

3) When the frozen field approximation is valid, the interaction of 
oscillatory velocities and an axisymmetric magnetic field leads to an oscilla- 
tory magnetic field which is parallel to the velocity field (Mallaus, 1967). 
Thus, there can be no dynamo action from the VXH terra. However, in the pres- 
sure of a non-axisymraetric toroidal field, there can he a net VXH term. 

4) A steady solution to this problem can be found for the eight magnetic 
field components given by combinations of poloidal-toroidal, dipole (il=l) - 
quadrupole (2.=2) , and axisymmetric (m=o) - antisymmetric (m=l) . In fact, these 
are two such solutions given by the relative poloidal field strengths at the 
core surface listed in Table 1, 
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Table 1. 

s'] 


S2 


Solution 1 

1.0 


0.3 


0.4 


Solution 2 
0.23 


0.8 


0.21 



The first solution corresponds to the normal arrangement of planetary magnetic 
fields-an axial dipole tilted by 16 with a fairly large relative quadrupole 
moment. The second solution is a nearly equatorial dlpole-whlch may be appro- 
priate for the ancient moon (S.F. Dermott, personal communication). 

The results of the calculations for the terrestrial planets are given in 
Table 2. 


Table 2. 


R core (cm) 

Mercury 

1.5X10® 

Earth I 
3.5X10® 

Earth II 
3.5X10.® 

Moon 

3X1o7 

Vr (cm/sec) 

lO""^ 

2X10“® 

1 . 3X10“^ 

2X102 

o 

surface g| (G) 

0.015 

0.005 

0.23 

100. 

Maximum toroidal 
field (G) 1.6X10^ 

180. 

9X10® 

2X10^ 


The results for Mercury show good agreement with those of Ness et al. 
(1975). The results for the Moon (at an assumed distance of 10 Earth radii) 
are also quite good (Runcorn, 1978). In the case of the Earth, good agree- 
ment is only obtained if it is assumed that the Earth's core response to the 
precesslonal force is enhanced by a factor of about 700 (Earth II). This is 
possible because the Earth's inner core, being rigid (Jacobs, 1975), should 
be quite spherical since it is loaded by a fluid of almost equal detSity 
(Houben, 1978). The near resonance of the spinover mode is then enhanced by 
several orders of magnitude over that for a spheroid of ellipticity 1/400 
(Earth I) . 
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THE EARTH'S CORE 

J. A, Jacobs, 'Dept. Geodesy and Geophysics, Univ. Cambridge, 

England. 

The major source of our information about the Earth's interior comes 
from seismology. Travel-time curves can be used to construct velocity- 
depth curves and from them the variation of density and other physical 
properties within the Earth. More recently, travel-time curves have 
been supplemented by free oscillation data. Considerable advances have 
been made in dealing with this inverse problem. Many Earth models 
have now been constructed and bounds placed on the variables. However 
the broad division of the Earth into crust, mantle and core has been known 
for some time - Oldham deduced the existence of a core in 1906. It 
must not be forgotten that seismology only gives us a "snapshot" of the 
interior of the Earth as it is today and gives no information about its 
structure in the past nor of its evolution. 

In 1936 , Miss Lehmann suggested that the liquid core contained an 
inner core which is solid. More recent evidence that the inner core is 
solid IS reviewed. The problem of the origin of the core is then 
discussed - which question cannot be divorced from the much broader and 
more difficiilt problem of the origin of the solar system and the 
constitution of all the planets. This leads naturally to the question of 
the thermal evolution of the Earth and a discussion of possible sources of 
heating. Apart from heating due to accretion and the impacts of large 
falling bodies, possible heat sources include the decay of long and short 
lived radio nuclides, adiabatic compression, tidal dissipation, induced 
electric currents and the release of gravitational energy due to the 
formation of the core. In addition to initial conditions and heat sources, 
the thermal regime of the Earth depends on the adiabatic and melting 
(or liquidus) temperatures. Our knowledge of both these temperatures is 
very rudimentary. 

Additional information about the Earth's core comes from its magnetic 
field. It is generally believed that the geomagnetic field is generated by 
some form of dynamo action in the core. Since rocks as old as 
-'^3, 000 m. yr have been found possessing remanent magnetism, a fluid 
(outer) core must have been in existence at least that long ago. A 
fundamental question in geomagnetism is what is the driving mechanism, 
of the geodynamo. If some form of convective motion in the outer core 
is demanded, this places further constraints on the thermal regime of 
the core and its chemical composition. Since dynamo theory is 
considered in other sessions, it is not discussed in any detail here. 
However the stability and possible stratification of the core is important 
in this connection and this question is examined more fully. Finally, 
the chemical composition of the core is briefly reviewed. 



THE ELECTRICAL AND THERMAL CONDUCTIVITY OF VARIOUS TERRESTRIAL CORE 
MATERIALS TO 1.8 MB’i'by R. Norris Keeler, George M. Matassov*, A. C. Mitchell 
and William J. Nellis, University of California, Lawrence Livermore Labora- 
tory, Livermore, CA 94550 

Electrical conductivity measurements have been extended to 1.8 Mb for 
the core materials Fe, Fe-Si and Fe-Ni. New measurements have been made for 
Fe-0, Fe-S and Fe-C; All conductivities fell within the currently accepted 
Magnetic Reynol ds - Number criterion for dynamo action. The wide limits and 
uncertainty of this criterion points to the need for better quantitative 
theories of the earth's dynamo to make high pressure conductivity measurements 
more meaningful in selecting the most probable core materials. 

*George M. Matassov is currently at the Air Command and Staff College, Mont- 
gomery, Alabama 

^Work performed under the auspices of the U.S. Department of Energy by 
the Lawrence Livermore Laboratory under contract number W-7405-Eng-48. 



FUTURE PLANETARY MISSIONS; Margaret Galland Kivelson, Department of 
Earth and Space Sciences and Institute of Geophysics and Planetary Physics, 
University of California, Los Angeles, CA 90024 

For almost four centuries after Gilbert recognized that the magnetic 
field of the earth approximated that of a magnetized sphere, the study of 
planetary magnetism was limited to the study of the earth. Terrestrial 
measurements demonstrated the complex and unsteady nature of magnetic proc- 
esses and defined an observational framework for theoretical analysis. The 
formulation of a complete and self-consistent theory of planetary magnetism 
in the future will undoubtedly be stimulated by comparative studies based on 
knowledge of the magnetic properties of other objects in the solar system 
including such diverse bodies as major planets, their satellites, comets and 
asteroids. Simultaneously, problems connected with remanent magnetism should 
be elucidated. From closeup spacecraft observations, we have already begun 
to establish the properties of Mercury, Jupiter, the Moon, Venus and Mars, 
although the interpretation of the magnetic properties of the latter two 
planets has been somewhat controversial. Indeed, there are many questions 
whose resolution must rely on future extraterrestrial missions. 

The table below lists missions which are expected to provide data on the 
magnetic properties of extraterrestrial solar system bodies in the near future. 
The list is chronologically ordered by launch date, with missions near the 
end of the table not yet approved but well into the planning phase. 


LAUNCH 

SPACECRAFT 

TARGET (S) 

COMMENTS 

1973 

Pioneer 11 

Jupiter (12/74), Saturn (9/79) 

Flybys 

1977 

Voyager 1 

Jupiter (3/79), Saturn (11/79) 

Flybys 

1977 

Voyager 2 

Jupiter (7/79), Saturn (8/80) 

Flybys 

1978 

Pioneer Venus 

Venus (12/78) 

Orbiter 

1982 

Galileo 

Mars (82) , Jupiter orbit (85) 

Satellite 

flybys 

1983 

VOIR 

Venus (83) 

Magne- 

tometer? 

1983 

Solar Polar (2) 

Jupiter (84) , Solar pole (86) 

Flybys 

1985 

Comet Encke 

Encke (87) 

Flyby 

1986 

Mars Geochemi- 
cal Orbit er 

Mars (87) 

Orbiter 

1987 

Saturn Orbiter 
Dual Probe 

Saturn (92) 

Titan, 

rings 

1988 

Mars Sample 
Return 

Mars (89) 

Magne- 

tometer? 


All of these missions will contribute to our understanding of planetary 
magnetism. Flights by Saturn will test the models we use to infer the sur- 
face field strengths of remote objects. Similarities and differences between 
the fields of Saturn and Jupiter will provide valuable Insight into the prop- 
erties of planetary dynamos. Orbits of Mars and Venus should help settle the 
question of whether either or both of these planets have intrinsic magnetic 
fields and clarify the nature of their interactions with the solar wind. Only 
peripherally noted in the table are the Galilean satellites which are objects 
of intense interest to anyone studying the formation of the solar system. 

These near Moon— sized bodies form a mini-solar system circling Jupiter. Close 
flybys planned for Voyager and later for Galileo will provide data on their 
interactions with the Jovian plasma and their intrinsic magnetic properties. 
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Ultimately, similar studies of the moons and rings of Saturn will supplement 
our knowledge of such bodies. However, it should be recognized that detailed 
mapping of surface anomalies familiar from Apollo investigations require 
multiple passes and will not result from flyby measurements. 

Not entered in the table are several missions on the "wish list" which 
have already been studied carefully. Examples includes a multi-asteroid 
reconnaisance mission which would look for remanent magnetic fields, a 
Galilean satellite lander, a flyby of Uranus whose magnetic dipole is expected 
to be almost in its orbital plane resulting in unique interactions with the 
solar wind, and a Saturn orbiter with probes. Finally, support for a lunar 
polar orbiter which could provide magnetic maps of the entire surface of the 
Moon is still strong in the scientific community and the European Space Agency 
is now considering the support of such a mission. 
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DY¥AMO MODELS FOR THE EARTH 

BASED ON MBAN-EIBLD MAGNET OHIDEODXNAMICS 

E. Krause and K.-H. Radler 

Zentralinstitut ftir AstropRysik, DDR - 15 Potsdam 

Some new results on dynamo models are represented whicii 
can be applied to the Earth, The models consist of a sphere of 
electrically conducting matter surrounded by vacuum. As usual 
in mean-field magnetohydrodynsmics both the magnetic field and 
the motion are understood as a sum of a mean and a fluctuating 
part. According to the assumptions on the mean motion two dif- 
ferent types of models are considered: one with a mainly rigid 
body rotation and one with a non-uniform rotation. As far as 
the mean electromotive force caused by the fluctuating motions 
is concerned in addition to the oc -effect also other effects 
are included. 

Models with rigid body rotation (which operate due to an 
oc^-process) make allowance for mean magnetic fields with dif- 
ferent symmetry properties with respect to the rotational axis 
and different time behaviour. In addition to axisymmetric sta- 
tionary fields there are uon-axisymmetric fields in the form of 
waves travelling around the equator. Also superpositions of 
different fields may occur. This leads to a possibility of not 
only explaining the existence of the dipole-like magnetic field 
of the Earth but also the deviation of this field from symmetry 
with respect to the rotational axis and its westward drift. 

In models with non-unifom rotation (which operate due to 
Oco or cjto -processes) axisymmetric magnetic fields are fa- 
voured, which may be either stationary or oscillatory. This 
provides for an interpretation only of that part of the magnet- 
ic field of the Barth which is symmetric with respect to the 
rotational axis. The behaviour of the remaining part can, e.g., 
be understood in terms of MAC-waves. 
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PLANETARY DYNAMOS— THE KINEMATICAL THEORY, E. H. Levy, Department 
of Planetary Sciences and Lunar and Planetary Laboratory, University of 
Arizona, Tucson, Arizona 85721. 

The primary magnetic fields of most large astrophysi cal objects, inclu- 
ding planets, are thought to be generated by a hydromagnetic dynamo process. 
In this process motions of an electrically conducting fluid interact with a 
magnetic field to generate electrical currents which produce the field. The 
dynamo phenomenon is described physically and mathematically by Maxwell's 
equations and Newton's laws, along with appropriate constitutive relations to 
describe the material properties of the fluid. Complete solution of the pro- 
blem, including simultaneously both the electromagnetic and mechanical be- 
havior, is a formidable task which has so far resisted definitive attack. 
However, tentative advances are being made against this problem based on the 
known solutions of simpler fluid dynamical problems and the results of kine- 
matical dynamo analysis of magnetic field generation. 

Essentially all of our understanding of the kinds of fluid motions which 
are capable of generating magnetic fields comes from the kinematical theory. 
In this approach, the problem is truncated by assuming that the fluid velo- 
city is given and fixed. Thus Newton's laws are dropped from consideration. 
Of course, most interest centers on fluid motions similar to those occurring 
in natural bodies. 

We review kinematical dynamo theory and its impact on our understanding 
of the magnetic-field generation process. Included is a discussion of recent 
and ongoing work, in this context we briefly review our ideas about the re- 
versal phenomena displayed by magnetic astrophysi cal objects. 



MEASUREMENTS OP LUNAR MAGNETIC FIELDS BY THE ELECTRON 
REFLECTION METHOD by R. P. Lin and K. A. Anderson, Space 
Sciences Laboratory, University of California, Berkeley, 
California 9^720. 

The electron reflection method for remote sensing of 
lunar surface magnetic fields depends on the fact that charged 
particles tend to be reflected from regions of increased magne- 
tic field strength. Charged particles in the cis-lunar environ- 
ment are guided to the lunar surface by whatever ambient 
magnetic field is present. If the field strength does not 
increase toward the surface, the particles will be absorbed by 
the surface material. If fields of near-surface origin are 
present, the total field strength will increase as the particles 
approach the surface, causing a fraction of the particles to be 
reflected back with an intensity that increases with the total 
field strength. Apollo 15 and l6 subsatellite measurements of 
the angular distribution of 0.5 and l4 keV electrons incident 
on and reflected from the lunar surface have been used to obtain 
high spatial resolution ('\/ 10 - 50 km) and sensitivity 0.1 y> 
1 Y = 10 ^ Oersted) maps of lunar surface magnetic fields. It 
is found that many patches of increased field strength ranging 
in size from < 10 km to > 200 km are distributed over the lunar 
surface. We have searched for correlations of these magnetic 
features to surface geology. In general there is no obvious 
correlation. The exceptions are as follows: 

1) Surface fields of > 100 y are present in a region of 
order 10 km wide and > 300 km long centered on and paralleling 
the long linear rllle Rima Sirsalis. The magnetic feature 
extends 60 km beyond the place where the rille disappears 
beneath the lava flows of Oceanus Procellarum. These magnetic 
results imply that the rill© has a strong magnetization 

(> 5 X 10~^ gausscm^ gm ) associated with it, either in the 
form of intrusive, magnetized rock, or in the form of a gap in 
a uniformly magnetic layer of rock. 

2) In lunar mare regions a statistical correlation is 
found between the surface field strength and geologic age of the 
surface as determined from crater erosion studies. 

3) However there is a lack of correlation of surface field 
with impact craters in the mare. This suggests that the mare 
do not have a strong uniform magnetization and that these 
impacts do not in themselves generate strong magnetization 
coherent over v lo km scale size. 

4) There appears to be a correlation of strong surface 
field regions on the lunar backside with ejecta from large cra- 
ters and basins. 

The implications of these findings on the origin and his- 
tory of lunar magnetism, and the application of this measurment 
technique to other planetary bodies will be discussed. 
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ARE PLANETARY DYNAMOS DRIVEN BY GRAVITATIONAL SETTLING? 
D. E. Loper, Geophysical Fluid Dynamics Institute, Florida State 
University, Tallahassee, FL 32306 and P. H. Roberts, University 
of Newcastle upon Tyne, Newcastle upon Tyne NEl 7RU, England. 

Several of the planets are known to possess magnetic fields. 
It is generally believed that these fields are sustained by 
dynamo action within the planetary interiors . One of the most 
challenging problems of planetary physics concerns the energy 
source for these dynamos. The two energy sources mentioned most 
often in the literature are thermal convection and precessionally 
induced motions, but each of these mechanisms encounters diffi- 
culties in efficiently converting the available energy into 
mechanical motion: in a metallic fluid such as the Earth's 

core, thermal convection may be short-circuited by thermal con- 
duction while the energy of precession appears to be dissipated 
in boundary layers. In this paper we shall elaborate upon a 
third possible energy source which is not subject to these pro- 
blems of efficiency: the energy released by gravitational set- 

tling. This idea was first proposed by Braginsky (1963) and has 
recently been discussed by Gubbins (1977) and Loper (1978 a, b) . 
A complete theory of the operation of dynamo driven by gravita- 
tional settling is currently being developed (Loper and Roberts, 
1978) . 

The basic features of graviationally driven convection may 
be summarized as follows. Consider an initially homogeneous 
self-gravitating sphere of molten binary alloy representing a 
planetary core, the alloy being composed of a heavy metal 
(principally iron) and a light nonmetal in the case of a ter- 
restrial planet or of helium and metallic hydrogen in the case of 
a Jovian planet. As the planet cools over geological time, a 
phase change is likely to occur, either to a new liquid phase 
immiscible with the first or to a solid phase. Assuming the 
adiabatic gradient to be smaller than the liquidus gradient, the 
phase change first occurs at the center of the sphere. It is 
known from metallurgy that the new phase which forms from a 
binary melt does not in general have the same composition as 
the liquid. To simplify the discussion we shall assume the melt 
to contain a larger fraction of the heavy component than does 
the eutectic. In this case the new phase has a larger fraction 
of the heavy component and hence is more dense than the initial 
liquid even in the absence of a density change associated with 
the phase change. The denser phase will accumulate at the 
center, forming an inner core. Associated with the growth of 
this inner core, there is a net motion of heavy material down- 
ward and light material upward, releasing gravitational poten- 
tial energy. More specifically as the inner core grows by 
accretion of the denser phase, a residue of light material is 
left in the liquid near the inner core. Since diffusion of 
matter is ineffective over planetary length scales, this excess 
of light material must distribute itself throughout the outer 
core by means of convective motions driven by compositional 
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bouyancy. The fact that the material dxf fusion coefficient is 
typically much smaller than the thermal diffusion coefficient 
accounts for the inherent efficiency of the gravitationally 
driven dynamo first pointed out by Braginsky (1963) and quali- 
fied by Gubbins (1977) . 

The power available from gravitational settling depends 
upon the density contrast Ap between the two phases and may 
by crudely estimated by the formula 

P = ApVgL/t 

where V is the inner core volume, L is the outer core radius 
g is a typical value of the local acceleration of gravity and 
t is the time which has elapsed since the inner core began to 
form. Typical estimates for the Earth’s core give P = 4 x lO^W, 
which is sufficient to sustain a large field (see Loper, 1978a) . 
In light of the uncertainty in composition and structure of the 
planetary cores, it is difficult at present to estimate the 
viability of gravitational settling as a power source for the 
other planetary dynamos. 

If gravitationally driven convection occurs, the outer 
core is in effect mechanically mixed and Loper (1978b) has 
shown that several interesting thermal regimes are possible. 

For example, if the liquidus gradient is less than the conduction 
gradient (i.e., that gradient which would occur in the absence 
of motion) , a slurry must occur in the fluid at the bottom of 
the outer core. A general thermodynamic theory for the motion 
of a binary-alloy fluid containing a slurry has been developed 
by Loper and Roberts (1978) . Also, compositionally driven 
convection can occur even if the actual temperature gradient is 
less than the adiabatic gradient, making the fluid thermally 
stably stratified. In this case the fluid motions convect heat 
radially inward — up the temperature gradient. Hence the heat 
conducted down the adiabat cannot be taken as a lower bound on 
the heat transferred out of the core. 

Braginsky, S.I., Structure of the F layer and reasons for con- 
vection in the earth's core, Doklady Akad. Nauk. SSSR , 

149 , 8-10, 1963. 

Gubbins, D. , Energetics of the earth's core, J. .Geophys. , 43 , 
453-464, 1977. 

Loper, D.E., The gravitationally powered dynamo, Geophys. J. Roy, 
astr . Soc . , to appear 1978a. 

Loper, D.E., Some thermal consequences of a gravitationally 
powered dynamo. J . Geophys . Res . , to appear 1978b. 

Loper, D.E. and P.H. Roberts, On the motion of an iron-alloy 
core containing a slurry. I. General theory, Geophys . 
Astrophys. Fl. Dyn . , 289-321, 1978. 



"MACRODYNAMICS OF DYNAMOS IN ROTATING STRATIFIED FLUIDS" 

Willem V.R. Malkus, Dept, of Mathematics, M.I.T., Cambridge, Mass. 02139 

Convective motions in the sun are presumed responsible for its magnetic 
fields, yet models for the solar stability require that these motions are 
confined to the outer sixth of its radius. In contrast, the magnetic fields 
and electric currents prevade the entire sun, and may be largest in the 
stable layer just beyond the lower boundary of the convection zone. Recent 
models for the stability distribution in the earth's core also suggest that 
the dynamo producing motions are confined to a small fraction of the 
electrically conducting volume. Theories for the geodynamo, e.g. Busse 1978 
Moffatt 1978, deal with magnetic fields generated in and, in turn, limiting 
the fluid motions responsible for them. No current theory deals with the 
role of an adjacent stably stratified electrically conducting fluid as the 
significant amplitude determining factor in a dynamo. Three earlier papers 
Greenspan 1974, Malkus and Proctor 1975, Proctor 1977) do discuss the 
equilibration of global magnetic fields by the large scale fluid motions 
resulting from the Lorentz forces. When ohmic loss is the principal 
dissipative mechanism and inertial forces play a secondary role, the cylin- 
drical integrals of the geostrophic components of the Lorentz force must 
vanish in a homogeneous fluid (Taylor 1963) . For magnetic fields resulting 
from averaged small scale motions (the -effect) , the Taylor constraint 
leads to an unusual eigenvalue problem in which the determination of a finite 
amplitude velocity field precedes the determination of the eigensolution for 
the magnetic field S . However the Taylor constraint is both inappropriate 
in a stratified flow and suspect in any intense regions of generation such 
as solar convection. Here the following idealization is proposed: consider 
a thin spheroidal shell with intense fluid motions characterized by a local 
; determine the appropriate reformulation of the Taylor-like constraint 
for a rotating stable stratified fluid; formulate the eigenvalue problem for 
0 ( o( ) and for the two-dimensional velocity field (eigenflow) required by 

the Taylor-like constraint; formulate the finite amplitude problem, including 
the reduction in the static stability by the finite 0 fields and also by 
large thermal conduction. First conclusions are presented for the finite- 
amplitude eigenflow and the structure of axisymmetric initial S fields for 
particularly simple . Eigenflows scale with the magnetic diffusion velo- 
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city as anticipated from previous studies. The amplitude dependence of 6 

on stability and conductivity is not yet determined. 
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"TURBULENT CORE FLOW DUE TO PRECESSIONAL AND TIDAL FORCES" Willem 

V.R. Malkus, Dept. o£ Mathematics, M.I.T., Camhridge, Mass. 02139 

The secular tidal deceleration of the earth is observed to he approxi- 

1 2 

mately three milliseconds per day per century ’ . In consequence, the 

secular loss of rotational energy by the liquid core of the earth is greater 

than 5.10^® ergs/sec.. Can the magneto-fluid dynamical tidal s^in-down, or 

3 4 

the interrelated precessional spin-over ’ , of the core fluid transfer this 
rotational energy to the fields and currents of the geodynamo? Laboratory 
studies to be reported here suggest that either core spin-dovm or core spin- 
over separately could power the geodynamo if the flows are fully- turbulent. 
How turbulent are the magneto-fluid dynamical processes in the earth's core? 
We may never establish the amplitude or time dependence of the small scale 
flows, however, the large scales exhibit considerable disorder. An example 
is the observed ten to thirty year random transfer of momentum between the 
core fluid and mantle. This process involves an exchange of energy^ with 
an r.m.s. value ten to one hundred times larger than the secular loss- of 
rotational energy. Hence this turbulent exchange would have to be Tuore than 
ninety percent (or ninety-nine percent) efficient to dissipate less than 
10^^ ergs/sec.. An electromagnetic core mantle coupling sufficient to 
explain these observed non-secular changes in the length of the day- remains 
a matter of debate^, to be briefly assessed here. The central purpose of 
this study is the determination of turbulent onset and transport "efficiency" 
in flows induced by precession and tides in laboratory rotating spheroids. 

The external and fluid parameters determining the transition of thes-e motions 
from laminar to turbulent are found for boundaries of varying roughness. 

The turbulent momentum transport is observed to relax the constraints on the 
flows imposed by rotation and the ellipsoidal cavity. In homogeneous or 
convectively unstable fluids the induced turbulent regimes occupy the entire 
spheroid. Initial experimental results on cylindrical ly stratified fluids 
suggest a host of new instabilities hut reduction in the r.m.s,. turbulent 
amplitudes. Over much of the parameter range estimated for the earth's 
core fluid there appears to be little doubt that both tidal and precessional 
forces could sustain turbulent flows akin to those observed in the 
laboratory. More complete dynamo theory, or an operating laboratory- dynamo, 
may be needed to establish that these flows actually act as energy^ sources 

for the geodynamo. 
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AN OUTLINE OF SOME OF THE PROPERTIES OF THE GEOMAGNETIC FIELD 

R.T. Merrill and M.W. McElhinny 
Research School of Earth Sciences 
Australian National University 
Canberra A.C.T. 2600 

The magnetic field measured at the earth’s surface originates 
from three sources: external sources that arise mostly from ionospheric 

currents, crustal remanent magnetization sources, and hydromagnetic sources 
in the earth's core. This talk deals with the latter, which can be well 
described outside the earth's core in terms of a magnetic scalar potential. 

Ip given by: 

ip = aZ E (cosS) (g? cos m(p + h? sin m<p) 

i=l £m=0 ^ ^ ^ 

x<;here x^e use spherical coordinates (r, 6, (p) , a is the earth’s mean radius 
(6371 km) , P^ (cos6) are the Schmidt polynominals and g^ and h^ are referred 
to as the Gauss coefficients. The magnetic field can be obtained, by 
taking the negative gradient of ip. 

Primarily through the use of measurements from magnetic 
recording stations, one finds that approximately 90% of the present 
magnetic field at the earth's surface can be described by a geocentric 
dipole inclined 11 - 5 ° to the axis of rotation, while 10% of the field is 
described by higher order (f->l) non-dipole field components. 

The field constantly changes x<fith time, x^ith the dominant 
feature being a x<restward drift of the non-dipole field at approximately 
0.2° per year (Bullard et dl. 1950). It is clear that this westward drift 
does not account for all the secular change in the field, since zonal 
harmonics (m = 0) also change significantly with time. A considerable 
amount of controversy exists concerning hox^ to best describe the changes 
in the field that are not drifting x^estward. Part of the problem is 
certainly the fact that the historical record represents only a very small 
time xflindow through x^hich the magnetic field can be observed. 

Palaeomagnetic studies allows one greatly to extend this 
window in time, but at the expense of considerable loss in resolution. 

Even xi7hen one confines xffork to rocks of the last 5 million years, for 
which tectonic problems are not too serious, one must deal xjith the fact 
that the palaeomagnetic record is discontinuous in both time and space. 

In spite of these problems, a crude first order, approximation 
for the "mean palaeomagnetic field" can be obtained. To do this we exclude 
times when the field is undergoing transitions between normal and reversed 
polarity states or times when the field is undergoing large changes in 
direction that are not associated with polarity transitions (excursions). 

The primary reason xAy one can obtain a mean palaeomagnetic field is that , 
to a first order approximation, the time averaged field turns out to be 
axisymmetric, probably reflecting the averaging done by the westx^ard drift 
of both the non-dipole and dipole fields. Because of this, one can average 
measurements from sites of a common latitude to surmount the problems that 
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result from having incomplete time averages of the field at any one given 
site 

The time averaged palaeomagnetic field obtained in this way in 
described by the zonal harmonics given in Table 1, as taken from Merrill 
and McElhinny (1977). This field exhibits north-south asymmetry. In 
addition, the time averaged normal polarity field is significantly different 
(at the 95% confidence level) from that of the reversed polarity field. 

Phillips (1977) has used an entirely different approach that 
also leads to the conclusion that there are significant differences between 
the normal and reversed polarity states. Phillips has shown that both the 
lengths of normal and reversed polarity intervals can be modelled by a 
gamma distribution, but that the shape of this distribution is significantly 
different (at the 95% confidence level) for the two polarity states. 

One can model the difference between the normal and reversed 
polarity states non-uniquely by a standing zonal magnetic field. One can 
then test the adequacy of such modelling by seeing if this standing field 
produces systematic virtual geomagnetic pole (VGP) paths during transitions 
between polarity states. Hoffman (1977) pioneered work on the systematics 
of VGP paths during transitions and showed that the observed systematics 
are consistent with kinematic dynamo models for polarity transitions 
developed by Parker (1969) and Levy (1972). Hoffman and ‘Fuller '(1978) 
have subsequently examined 8 quadrupole and octapole models that in a 
generalized sense are consistent with the Levy-Parker reversal mechanism. 

Of these, 2' models appear to be most consistent with the VGP transition 
data. The first is a Quadrupole model in which reversal occurs fxrst in 
the southern hemisphere and the second is an Octapole model in which 
reversal occurs first at low latitudes. These are precisely the -two models 
that are predicted from the use of the standing zonal field discussed 
earlier; 




POLARITY 


TERM 

0 

§2 

NORMAL 

REVERSED 

PRESENT FIELD 

“o 

0 

§3 

0 

§1 

+0.050 

+0.083 - 

+0..063 

-0.017 

-0.034 

-0.043 


Table 1 . Gauss coefficients for the time averaged palaeomagnetic field. 
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METEORITE MAGNETISM AND THE EARLY SOLAR SYSTEM 
MAGNETIC FIELDS ( A REVIEW ) : Takes! NAGATA, National 
Inst. Polar Res., Kaga 1-9-10, Itabashi, Tokyo 173, Japan. 

Since the first systematic studies on meteorite magnetism by Stacey et al. 

{ 1961 ), a number of research works on the magnetic properties of meteorites 
and their implication for the early solar system magnetic fields have been review- 
ed, particularly by Brecher ( 1971 , 1977 ) and Gus'kova ( 1972 ). The main 
problems in meteorite magnetism may be classified into (a) basic magnetic pro- 
perties in relation to their chemical and petrographical compositions and (b) 
characteristics of natural remanent magnetization ( NRM ) of meteorites in asso- 
ciation with the early solar system magnetic field. 

(a) Basic magnetic properties of meteorites 

Intrinsic magnetic properties of meteorites such as saturation magnetization 
( Is ) and magnetic transition temperature ( @ ) represent contents of ferromag- 
netic constituents and their compositions. Corresponding to the Urey-Craig 
diagram, Is values of ordinary chondrites and E-chondrites decrease in the 
order of sequence expressed by E->H-^L-»LL, whereas Is values of achondrites 
are extremely small ( Is ^ 1 emu/ gm ) because of a very small content of metallic 
phase. Ig values of C-chondrites are mostly due to spontaneous magnetization of 
magnetite or magnetite plus a small amount of native iron ( e.g. Larson et al. , 19 
74; Watson et al. , 1975; Herndon et al. , 1976 ). 

The composition of ferromagnetic phase in meteorites is reasonably well re- 
presented by thermo magnetic characteristics which separately give magnetizations 
of thermally irreverssible Oi-phase, a thermally reversible ^-phase and a 
thermally unstable (CX+ /)-phase of Fe-Ni-Co-P alloy grains and thermally re- 
versible magnetites. Corresponding to the Prior rule, ratio of saturation mag- 
netization of Ot-phase component, Is( Ol), to the total Ig value decreases in the 
order of sequence given by E-5>H-»L-?»LL-^4!J . 

In Fig. 1 , observed 
values of Ig and Is( ^ )As 
are plotted on an Ig vs 
Is( <5^ )As diagram, which 
gives rise to a reasonably 
clear magnetic classification 
of stony meteorites corres- 
onding to their chemical 
classification . 


The structure-sensitive 
magnetic properties of 
meteorites such as magnetic 
susceptibility ( yt^o), satu- 
ration remanence ( Ir ) , 
coercive force ( He ) and 
remanence coercive force 
( Hrc ) also have been exam- 
ined. A significant charac- 
teristic of meteorite suscep- 
tibility is its considerably 
large anisotropy { -iZ^max/ 



Fig. 1 
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^mm - 1.5 ''^2.0 ), which suggests a possibility of their aggregation in the early 
solar system ( e.g. Stacey et al. , 1961 ; Brecher and Arrhenius, 1974 ). 

( b ) Natural remanent magnetization 

Among various kinds of meteorites, C-chondrites generally maintain stable 
NRM (Banerjee and Hargraves, 1972; Butler, 1972; Larson etal,, 1973; Brecher 
and Arrhenius, 1974 ). Since it is believed that C-chondrites were formed in the 
early stage of solar system formation, NRM of C-chondrites may represent the 
interplanetary solar wind magnetic field at the stage . The paleointensities ( Hq ) 
of C-chondrites derived from their NRM analysis are summarized in the following 
table, where To denotes a critical temperature, during a cooling process from 
which the remanent magnetization was acquired. 


meteorite 

1 Ho ( Oe ) 

To ( °C ) 

Inve stigators 

Orgueil ( C-| ) 

0.67 

120 

Banerjee & Hargraves 

Renazzo ( C 2 ) 

2.3 

( 250 ) 

Brecher 

Murray ( C 2 ) 

0.7 

( 250 ) 

Brecher 

Murchison ( C 2 ) 

I 

[ 0.18 i 

90 

( 250 ) 

Banerjee & Hargraves 

1 

[ 0.4 ^ 3.0 

Brecher 



r 1 .09 

130 

Banerjee & Hargraves 

Allende ( C 3 ) 


1 .3 

0. 15^-0. 95 

150 

Butler 

Brecher 

Leoville ( C 3 ) 

0.5 

_ 

Nagata & Sugiura 

0.97 

350 

Nagata & Sugiura 


Achondrites also maintain stable NRM in most cases. Since it is believed that 
achondrites represent the upper mantle part of meteorite parent planets, their 
paleointensity may suggest the magnetic field of parent planets . Paleointensities 
of Yamato-7307 ( howardite ) and Yamato-7401 3 ( diogenite ) are 0.13 and 0,093 
Oe respectively. 

NRM of E-, H-, L- and LL-chondrites are comparatively unstable in most 
cases. However, some of these chondrites maintain reasonably stable NRM, 
from which the paleointensity can be estimated. Estimated paleointensities are 
0.10 -^0.68 Oe for L-chondrites, 0.15 "^ 0.39 06 for H-chondrites and 0.33 
Oe for an E-chondrite ( Stacey et al. , 1961; Gus'kova, 1963, 1970; Nagata and 
Sugiura, 1977 ). However, Brecher and Ranganayaki ( 1975 ) have stated that 
paleointensity of H-, L- and LL-chondrites is much smaller, ranging only from 
0.01 to 0.08 Oe. 

Iron meteorites also have stable NRM in most cases. It seems, hpwever, that 
their stable NRM is due to spontaneous magnetization of (1 1 0} plates nucleated 
and developed along the octahedral /fill] crystallographic plane ( Brecher and 
Albright, 1977 ). 
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THE ELECTRICAL CONDUCTIVITY AND EQUATION OF STATE OF SEVERAL CON- 
STITUENTS OF THE MAJOR AND OUTER PLANETS by Wil liain J. Nellis, Arthur C. 
Mitchell, Maxim I. Kovel*, Ronald S. Hawke and R. Norris Keeler, University 
of California, Lawrence Livermore Laboratory, Livermore, California 94550 

The primary constituents of the major and outer planets are He, Ho, N 2 , 
CH-, NHo and HoO. Equation of state measurements have been made on N 2 to 
600 kb; CH 4 to 400 kb; NHo to 700 kb, and H 2 O to 2.2 Mb. Electrical conduc- 
tivity has been measured for H 2 at 2 Mb; NH 3 to 280 kb; and HoO to 450 kb. 
None of these materials exhibited anomalous behavior within the limits of 
the measurements made with the exception of N 2 , which shows behavior consis- 
tent with transition to an atomic (metallic) phase. The implication of these 
results is that the level of conductivity which exists in hydrogen above its 
transition pressure (2 Mb) does not exist in the most abundant constituents 
of the outer planets at core pressures; therefore, it may be anticipated 
that neither Uranus nor Neptune will have a magnetic field, unless it is 
associated with a small residual metallic core. 

*M. I. Kovel is currently at the Office of the Chief of Engineers (Research 
and Development Branch), Department of the Army, Washington DC 20360 
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Direct observations of the magnetic field of Mercury and its magneto- 
sphere, formed by the interaction of the solar wind, were performed twice by 
the Mariner 10 spacecraft, in March of 1974 and again in March 1975« From 
these data, it is clear that there exists an intrinsic magnetic field of the 
planet, sufficiently strong at present to deflect the solar wind flow around 
the planet and to form a detached bow shock wave in the super Alfvenic solar 
windo Three different methods have been used to analyze the magnetic field 
and derive quantitative values for the description of the planetary field: 

(1) Direct spherical harmonic analysis of the data, 

(2) Modeling of the magnetosphere by an image dipole and infinite current 
sheet in addition to the planetary field, 

(3) Scaling of a mathematical model for the terrestrial magnetosphere. 


The results obtained yield dipole moments \/(gl*^)^+(hl^)^+(gl^)^ ranging 

from 2.4 to S.lxlO^^Gauss cm^, with the lower values associated with certain 

^ n 0 

models using partial quadrupole (g 2 ^) snd octupole (83 ) terms to improve 
the least squares fitting of models to observations. Because the data set 
is incomplete, in the mathematical sense, no unique representation of the 
planetary field multipolar representation can be derived by method (1). The 
use of only 1 of the 5 quadrupole moment terms and 1 of the 8 octupole 
moment terms corresponds to a displacement of the dipole along its axis. 

These terms, used in methods (2) and (3), yield equivalent offsets of the 
dipole by approximately 0.2 Rjyj. The selection of only those higher order 
terms having axial symmetry (m= 0 ) has not been justified on any physical 
grounds. The fields of Earth and Jupiter do not show axial symmetry. Such 
a large offset may reflect the limitations of the models used to represent 
the external current systems. Because of the relatively short radial 
excursion of the data, the and terms may also be spatially aliased 
with the gj^® term. 


Analyses by method (1) of subsets of data from the third encounter, 
taken near closest approach, yield a convergent series of dipole moment 
values which are believed to best represent the intrinsic planetary field. 
These provide a mean moment of 330(+ 18)y = 4,8xlO^^Gauss cm^ at a tilt 

angle of 14°+5 and at longitude of 148°+21°o ttiis means that the surface 
field at Mercury is about 1 % of Earths, while the moment is 6x10“^ of Earths. 


The origin of the field cannot be uniquely '•determined. It may be due 
to an active dynamo, a remanent magnetic field or a combination. Consider- 
ation of remanence as the source leads to some difficulties although 
definitive knowledge of the planetary Interior structure and thermal state 
is lacking sufficient to categorically eliminate this source. Some success 
in attempting to explain the field as due to an active dynamo has encouraged 
these efforts. This paper will critically review the quantitative analyses 
of the Hermean field and discuss its origin. 
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TABLE I 


MERCURY'S MAQilETIC MOLTIPOLE MOMENTS* 


Source 

Data 

0 

h 

ro 

o 

0 : 

§3 

External Terms 

Offset 

Tilt* ** 

Ness et al 
(1974a) 

I 

227 

0 

0 

0 

0.47 


Ness et al 
(1975 

I 

350 

0 

0 

n=2 

0 


Ness et al 
(1976) 

III 

342+15 

0 

0 

n:=l 

0 


Whang 

I & 

266 

0 

0 

2-Diinens ional 

0 


(1977) 

III 

165 

117 

0 

Tail sheet + 

0 

2,3 ° 

i 

1 


166 

75 

48 

image dipole 

0 

0 

Jackson and Beard 

I & 

257 

0 

0 

Scaled terrestrial 

0 


(1977) 

III 

170 

114 

0 

analogue 

0 

0 
o 

1 

l-» 

o 

o 

Ng and Beard 

I & 

190 

0 

0 

Scaled terrestrial 

(0.033, 0.026, 

1.2 ° 

(1978) 

III 




analogue 

0,189) 

0 

Ness (1978) 

III 

330+18 

0 

0 

n=0 and 1 




*Using dipole aligned coordinates. Polarity sense is same as at Earth. 

**Relative to normal to ecliptic or orbital plane, indicated by subscript E or 0. 
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THE METALLIC IRON IN LUNAR SOILS AND BRECCIAS, G.W. Pearce, 
Department of Geology, University of Toronto, Toronto, Ont. MSS lAl, Gan. 

It has for some time been known that the characteristic materials of 
the lunar regoliths — soils and soil breccias — differ markedly from 
local rocks in the amount and grain size of constituent metallic iron, the 
Only important magnetic mineral of the lunar surface. Thus at mare sites 
the regolith material is enriched by factors in the order of 10 in metallic 
iron over the local basaltic rocks (1, 2, 3). Much of the excess iron metal 
is very fine grained (dimensions less than 300 k) as determined magnetically. 
Mare soils are quite similar to each other both in the quantity and grain 
sizes of their metallic iron components. At highland sites the local rocks 
have perhaps, on the average, more metallic iron than the regolith constit- 
uents (4) , However, as with the mare regoliths the metal in the highland 
soils and soil breccias is generally fine grained whereas highland crystal- 
line rocks contain mostly coarse iron (dimensions greater than 300 A) (4) . 

The highland soils are also quite similar in most cases to each other in 
their metal contents. 

The highland and mare soils appear to differ in how they obtained the 
metallic iron they possess. Using a combination of magnetic analysis and 
trace element analysis, Chou and Pearce (5) were able to determine the 
relative contributions of iron from external sources such as meteorites and 
of iron produced by reduction of indigenous silicates and oxides. Thus 
highland soils appear to have obtained almost all metal from meteorites since 
they contain the proper meteoritic abundances of nickel and other trace 
elements. Mare soils, on the other hand, have a considerable contribution 
(approximately one half) from reduced iron. 

In a detailed analysis of grain size separates of an Apollo 15 core (6) , 
it was found that the sub-45 pim fractions are always considerably richer in 
metallic iron than fractions of larger grain size. The concentration of 
metallic iron in these sub-45 ym fractions showed some variation as a 
function of their depth below the regolith surface, but showed remarkable 
similarity in these magnetic properties (such as coercivity and the ratio of 
saturation remanence to saturation magnetization) which are sensitive to the 
grain size distribution of the iron particles themselves. This thus suggests 
that there was a common source during the period of time represented by the 
core or that there was a process acting on the metal modifying its grain 
size towards a constant distribution. The similarity of mare soils to each 
other and of highland soils to each other further suggests that this common 
source or equilibrating process may be widespread geographically as well as 
temporally. 

The major portion of the external source of iron metal is the rain of 
vaporized and melted material arising from meteorite, micrometeorite and 
solar wind particles that hit the lunar surface. Most of the metal in the 
rain would be due to vaporization and melting of meteoritic metal although 
a small component due to vaporization and melting of indigenous metal would 
be present. The rain would consist of the vaporized and melted material 
returning to the lunar surface as solid particles of small but variable size. 
Meteoritic iron whxch remained solid during impact and fallback must be a 
small component of the rain since the iron metal in meteorites is generally 
much coarser than that found in lunar soils . Vaporized metal would now be 
found on the surfaces of small soil particles whereas melted metal might 
sometimes be found inside silicate particles which had simultaneously been 
melted and ejected. If the distribution of mimpacting particles is constant 
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in time and space then the grain size spectrum of the metal in the fallback 
might be expected to be constant from place to place and from time to time. 
Thus there is a possible common source for iron In highland sites. 

In the case of mare sites the indigenous metal must be considered. 
Reduction of iron from the ferrous state on the lunar surface must be 
accomplished through the agents of heat from impacts and reducing gasses 
such as H2 and CO Implanted into regolith particles from the solar wind. To 
get the very fine grained metal particles observed in soils, the reduction 
must take palce under constrained conditions of temperature and heating time 
(7) , since at high temperatures (above 600-700'’C) grain coarsening of the 
iron particles occurs through diffusion (8) . 

Whether the iron metal results from meteoritic addition or reduction 
of indigenous material the latter process of grain coarsening will continue 
to occur as long as impacts produce heating of the regolith. Thus to get 
the uniform distribution of grain sizes of iron particles observed or 
expected, an equilibrium of between the processes of addition of very small 
particles and the grain growth of these particles would have to exist. 
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Magnetic Properties of Asteroids 
John Li. Remo 

Department of Geology, Hofstra University, Hempstead, New York 11550 

Abstract 


The aim of this research is to develop a first order set of magnetic 
models of asteroids based on the magnetic properties of meteorites 
determined from laboratory analysis. Although asteroids have been sus- 
pected as the source of meteorites by many investigators (Opik, 1963; 
Anders, 1964; Arnold, 1965; Wetherill and Williams, 1968), the association 
of asteroid mineralogy with certain corresponding meteorite types has been 
accomplished through telescope spectral data (Me Cord et. al. , 1970; 

Me Cord and Gaffey, 1974). Further studies of the reflectance spectra of 
asteroids and meteorite optical properties has indicated that some meteor- 
ite types have reflectance curves which maybe interpreted in terms of the 
type, composition, relative abundance and distribution of the constituent 
mineral phases. Such matching of meteorite and asteroid surface reflect- 
ivity data suggests the following associations: Vesta -basaltic achondrite; 
Ceres- carbonaceous chondrite; Toro-chondrite; and Juno-iron with a 
silicate admixture. We assume, for the first order models, that the sur- 
face material of these asteroids is representative of their internal constitu- 
tion. 


We consider asteroids composed of one or more layers of meteoritic 
material which possesses its own magnetic moment and which is subjected 
to an external (solar) field. The total magnetic potential, ^ consists of a 
natural magnetic remanence term, and an induction term, 415 , such that 
Laplaces equation is solved in the region where B =jUH and where 
there are no currents. In this region the magnetic field II_ and induction 
are derived from the (scalar) magnetic potential, in regions where ferro- 
magnetic behavior is exhibited, the permeability, , and magnetic induc- 
tion become so large as to permit simplifying assumptions. Our approach 
does not take into account any details of the thermal remanent magnetism 
acquired by a body cooling through the Curie point, but rather uses the 
measured values of the natural magnetic remanence and so the problem is 
simplified. 

We use the laboratory results of Guskova (1972) and Brecher 
( 1975 , 1977 ) who have documented the magnetic properties of numerous 
stony and iron meteorites. Associations between the magnetic properties 
of several meteorite groups has been established by Wassilewski (1974). 


Results based on remanence data indicate homogeneously layered 
stony asteroids possess surface fields from 12 - 2400 gamma at the poles 
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6 - 1200 gamma at the equator. Asteroid models composed of iron 
meteoritic material give unrealistically high values starting from 6000 
gamma to well above 50, 000 gamma at the equator, and double at the poles. 
In addition, the inductance field is vectorially added to the remanence field. 

There are several criticisms of this model; the magnetic remanence 
is assumed to have a single orientation throughout the body, the meteor- 
ites past history affects the remanence properties, the asteroid is 
probably not unformly magnetized -- we may be only observing (asteroid) 
crustal magnetic remanence, and the overall population of meteorite 
classes is not the same as that of observed asteroid spectral data. 
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SOME COMMENTS ON METEORITE MAONETISM: M.W.Rowe and M. Hyman, Dept, of Chem. . 
Te>^as Only., College Station, TX 77843. 

Primordial Magnetic Fteld- During the past decade, the consensus of 
opinion for the origin of the magnetization of stone meteorites has shifted 
from fields associated with meteorite bodies to fields associated with the 
solar nebula during the early period around the formation of the solar 
system fe.g., Stacey and Banerjee, 1974). Evidence toward this conclusion 
rests primarily' on magnetic studies of carbonaceous chondrites (Banerjee and 
Hargraves, 1972; Brecher, 1972 Butler, 1972; Brecher and Arrhenius; 1974). 
Qeochemical evidence from carbonaceous chondrites indicates that they are 
primary condensates which accreted at less than 200“ C and have not been 
heated significantly since accretion (Anders, 1971). Thermomagnetic 
analyses on all known carbonaceous chondrites at the time (Larson, et al . , 
1974; Watson, ^aj^. > 1975; Herndon, ^ a]^. , 1976) indicated that some of 
the carbonaceous chondrites are more suitable for paleointensity determin- 
ation than others. AH Cl chondrites contain only magnetite as a magnetic 
constituent and their magnetic analysis may be more straightforward than the 
carbonaceous chondrites studied to date, which contained troilite. Troilite, 
under the conditions generally used for paleointensity estimates, seriously 
complicates the analysis. Other C2 and C3 chondrites contain magnetite only, 
Viherea^ some others contain primarily magnetite and iron. This information 
should be useful in designing magnetic experiments. Unpublished thermo- 
magnetic analysis of some ureilites are presented here as well. Studies on 
carbonaceous chondrites, especially those on carbonaceous chondrite chondrules 
CLanoix, , 1978) yield information about the fields associated with the 

solar nebula at very early times in the history, i.e., about the first 10 
million years CSonnet, 1978). Estimates of the intensity of this solar 
nebula field range for M-16 0e. 

Later Magnetic Fields - Rowe 0974 a,b) suggested that in some heavily 
shocked meteorites sufficient heating was experienced relatively recently 
in solar system history to have erased all evidence of any magnetization 
acquired prior to the collision. Careful magnetic studies on such 
meteorites may lead to interesting information about magnetic fields of 
more recent origin. 
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Bogard , (1976) have shown by ^^Ar-^^Ar dating that at least 

nine chondrites show ages of < 0.6 X 10^ years which resulted in 89-985^ 

argon loss. Rowe (1974 a,b) previously argued that such loss probably 

requires temperatures in excess of any expected Curie point involved. Of 

particular interest is the Ucera chondrite which Bogard, ^al_. , found to 
40 39 

have a Ar- Ar age consistent with its cosmic ray exposure age. Magnetic 
studies, if properly conducted (Brecher and Ranganayaki (1975), on Ucera and 
other shocked meteorites should yield considerable information on the dura- 
tion of the magnetic field(s) which have been recorded by meteorites which 
have undergone sever shock. 

Magnetic and Troilite Contents of Selected Carbonaceous Chondrites- 
Hyman, ,(1978) have conducted large numbers of measurements on the 

magnetite contents of the Cl chondrites. Those measurements, along with a 
fevt troilite estimates, will be presented here. 
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THE ANCIENT LUNAR DYNAMO; S.K, Runcorn, Institute of Lunar & Plan- 
etary Sciences, University of Newcastle upon Tyne, England. 

The discovery of lunar palaeomagnetism was a surprise because (a) Luna II 
failed to detect a general lunar magnetic field, (b) the mean density of the 
Moon (3.34) was rather generally held to exclude the existence of an iron 
core, and (c) the absence of indications of tectonic movements, which had been 
generally interpreted by a Moon with a cold accretional origin, a minimal 
later differ ention or no internal dynamical history. Some carry over of 
these early ideas has contributed to scepticism concerning the interpretation 
of the facts of lunar palaeomagnetism. As in the early development of 
terrestrial palaeomagnetism, and indeed of any other historical process, 
hypothetical ideas based on limited observations to explain the origin of the 
phenomena are hard to test or to exclude with finality even if physically 
implausible and the correct interpretation is that which explains best the 
sura total of the evidence. Just as terrestrial palaeomagnetism was seen to 
contain a record of the variations in the geomagnetic field, so we think that 
the remanent magnetization of rocks records the presence in the Moon's early 
history of a magnetic field of internal origin. 

Stable remanent magnetization ranging from 10 ^ emu/gra to 10 ^ emu/gm has 
been found in crystalline rocks and high grade brecchias from all the Apollo 
missions, the magnetization being carried by multidomain and single domain 
iron grains. From this random sampling of the Moon, widespread general 
magnetization of its crust must be inferred and fits rather well the 
discoveries of magnetic anomalies from Explorer 35 studies of the solar wind, 
by the magnetometers landed on the surface during the Apollo 12 and 15 
missions and the surveys made by the astronauts during Apollo 14 and 16 
missions, the magnetic surveys made from the sub-satellites of Apollo 15 and 
16 and the magnetic anomalies inferred by the reflected electron technique. 

An essential idea in the interpretation of these anomalies is that their 
magnitude and scale, 1-lOOT and tens of kilometers are consistent with an 
origin in edges, such as are produced by cratering, in horizontal strata 
magnetized with intensities of those of the returned samples. 

The Moon's present dipole moment is extremely small (the surface field 
is 0.05T*). This is strong evidence against a uniform magnetized crust of 
more than negligible thickness; thus this null result excludes magneti- 
zation in a roughly uniform external field such as might be produced by a 
close rapid approach to the Earth or some primeval magnetic field perpendi- 
cular to the ecliptic. The crustal magnetization of the Moon is, therefore, 
either random or is the result of magnetization by a field of internal 
origin, for it has been proved that if a spherical shell became permanently 
magnetized in the same direction and proportional to the intensity of any 
field of internal origin, this magnetizing field later disappearing, then, 
induced magnetization effects being neglected, there is no field whatsoever 
outside it. On this view the lunar anomalies arise from small scale edge 
effects or demagnetization by impact. Recent analyses of the limited region 
of the Moon where magnetic anomalies have been carefully determined support 
the hypothesis of magnetization by an internal field rather than the hypo- 
thesis of random magnetization such as would result from impact magnetization. 

Pa laeoin tensity measurements have been made with apparent success on a 
number of Apollo rocks but in other cases inconsistent results have been 
obtained which suggest that either secondary magnetization or chemical changes 
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The Ancient Lunar D 3 raamo 
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during experimentation prevent the original magnetic field intensity from 
being determined. Fig. 1 shows results by the Thellier-Thellier method and 
the ARM method which appear to be reliable and the range of those found by the 
IRM normalization method. Both the high fields of around 1 gauss for 
4.0 109 yrs ago and the suggested exponential decrease with time are 
significant. This later treni^ if substantiated, is extremely hard to explain 
on the impact hypothesis. The high field values also exclude a primeval 
permanent magnetization of the Moon as being the source of the magnetizing 
field. The hypothesis of a core dynamo is, therefore, the most adequate 
explanation of the whole lunar palaeomagnetism data. 

The evidence .for a present lunar core is not conclusive but comes from 
three lines of evidence; A delay in the P-wave travel time passing through 
the core from a meteorite impact on the far side, electrical conductivity 
measurements indicating a possible high conductivity region in the deep 
interior and recent determinations of the moment of inertia factor as 
0.391 + 0.002. It can be shown, from the absence of changes in the dimensions 
of the Moon since about 4.0 10^ years, that the present core must have formed 
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prior to this time. The formation of a core in the very early history of the 
Moon requires a powerful heat source and the melting of the Moon to the centre 
can not be done by accretion or eddy current heating. The hypothesrs has 

been made that superheavy elements are responsible for the complete melting of 
the Moon. It is likely that these are siderophile and, therefore, could 
provide the heat source which generates the ancient lunar magnetic field. 
Thermodynamic arguments enable .a minimum value of the heat source in a core of 
500 km in radxus necessary to generate the fields observed is too high (at 
4 , lO^ yrs the heat source must be 10^2 be generated by the conventional 

methods suggested for the terrestrial dynamo. Fig. 2 shows the rate of 
fission of SHE nucleii needed for different periods and the results give a 
half life of 10® years (dotted line). 


Fig. 2 

WKI 

^/yr/kg of primeval material 
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THE MARTIAN MAGNETIC FIELD, C.T. Russell, Institute of Geophysics 
and Planetary Physics, University of California, Los Angeles, California 90024 

The first spacecraft to Mars, Mariner 4 in July 1965, carried a magne- 
tometer as well as other particles and fields instrumentation that would have 
detected a Martian magnetosphere had it encountered one. However, at its 
rather distant 'closest' approach of 3.9 Rm, it detected no magnetosphere. 
Mariner 4 did encounter a bow shock similar to the earth' s^. This observation 
led to some controversy over the existence of a planetary field due to the 
uncertainties in determining the size of the obstacle needed to cause a shock 
at the observed position's 3, 

The next set of observations came with the successful orbiting of Mars 
by the Soviet Mars 2 and 3 spacecraft in November and December 1971‘^ and later 
by the Mars 5 spacecraft in February 1974^. These spacecraft had periapsides 
of 1100 km for Mars 2 and 3 and 1800 km for Mars 5. Prom the Mars 2 and 3 
data, Dolginov and colleagues immediately concluded that Mars had an intrinsic 
magnetic field^. They confirmed this conclusion x-zith the Mars 5 data^. The 
identification of a planetary field by Dolginov and colleagues is essentially 
based on three different types of observations: the location of the boxi? shock 
crossings identified by a variety of instruments; a direct penetration (or 
two) into the dayside magnetosphere; and the constancy of the direction of the 
magnetic field in the wake region of Mars. Any one of these observations 
would be very good evidence for an intrinsic planetary field, if it were un- 
ambiguous. However, all three observations have been called into question. 

When the totality of bow shock observations are considered, and the ter- 
restlal scaling between magnetopause and shock location is used to scale be- 
tween Martian shock and obstacle height, the effective obstacle height is in 
the Martian ionosphere^} 7, Some of the distant "shock" encounters may have 
been confused xjith upstream effects on the solar wind8j9. 

If the scaled magnetopause height is in the ionosphere at a height of 
300 km then it is very improbable that Mars 2 and 3 entered the magnetosphere 
proper above the dayside of the planet since their closest approach to the 
planet was 1100 and the claimed encounters were even more distant than this. 

In fact, if the two reputed distant shocks were upstream effects, then the 
magnetopause encounters were in fact simply shock crossings^* 9, Further, 
the field in the putative magnetosphere had the characteristics of shocked 
solar xirind field lines draped around the planet®. 

Finally, the reputed Martian magnetotail observations resemble neither 
those of the Earth nor of Venus. The magnetosheath field points directly into 
the wake of the planet as if there were no obstacle to such penetration as a 
magnetotail would providel®. The field in the wake region Is claimed by 
Dolginov et al. to be independent of the radial component of the interplanetary 
field. However, this is not the component of the solar wind magnetic field 
to which induction effects are sensitive. In fact all tail observations are 
consistent with induction effects^®. 

The estimates of the Martian magnetic moment are all much less than one 
xrould expect if Mars had an active magnetic dynamo^l. In the absence of any 
theory one would expect since the earth and Mars have similar rotation rates 
and since Mercury and Mars are expected to have similar core sizes, the 
Martian moment should be intermediate between those of the Earth and Mercury, 
i.e., between 8 x 1025 anj 3 jj 1Q22 Gauss-cm^. If the dependence on rotation 
rate were linear we would expect the Martian moment, as scaled from Mercury 
to be about 3 x 1024 Gauss-cm®. This estimate is far larger than that reported 
by Dolginov et al.^>3 of 2.5 x 1022 Gauss-cm3, or the upper limit based on 
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rejecting the reported Mars-5 magnetotail encounters^^ of 2.1 x 10^1 
Gauss-cm3. 

If we use the similarity law for planetary dynamos derived by Bussell 
and assume a core of about 1750 kml3 and that internal core conditions are 
similar, we obtain a predicted moment of about 5 x 10^4 Gauss-cm^ a factor of 
200 greater than Dolginov et al.'s estimate and 2000 greater than our esti- 
mated upper limit. We note that Busse’s model when applied to Venus in a 
similar way is only about a factor of 3-4 too high. Another way to approach 
this calculation is to use the similarity law to calculate the dynamo radius. 
Using Dolginov et al.'s estimate we obtain a radius of slightly under 500 km. 
Using our best guess upper limit we obtain a maximum dynamo radius of 250 km. 
Thus, at most only a small fraction of the Martian core, as detemined from 
gravity studies, can be undergoing djmamo action today if Busse's model of 
the geodynamo is correct, or applicable to Mars. Even without the guidance 
of theory, the possible strength of a Martian dynamo is far below expectations. 
Hence, it is most likely, that if Mars ever had an active dynamo, it is not 
operative today. 
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FREEZING THE CORES OF TERRESTRIAL PLANETS BY SUBSOLIDUS WHOLE MANTLE 
CONVECTION. G. Schubert, Dept, of Earth & Space Sciences, Univ. Calif., Los 
Angeles, Calif., 90024, P. Cassen and R. E, Young, NASA Ames Research Center, 
Moffett Field, Calif. 94035. 

It is all too easy for subsolidus whole mantle convection to freeze the 
core of a terrestrial planet thereby destroying its magnetic field. This 
was first demonstrated quantitatively by Schubert and Young (1976) who calcu- 
lated temperatures in a convecting, constant viscosity, internally heated 
fluid model of the Earth's mantle. They showed that the temperature at the 
core-mantle boundary would lie significantly below the iron melting point if 
the mantle viscosity were less than lO^^cm^/s. Cassen et al. (1976) showed 
that convection in the relatively thin mantle of Mercury could freeze its 
core in a billion years or less . 

Here we report additional computations using the Earth model of Schubert 
and Young (1976) which reveal that the low core-mantle boundary temperatures 
are established by the convecting mantle on a time scale less than the age of 
the Earth. If the viscosity of the Earth's mantle is indeed nearly uniform 
with the value of 10 ^cm^/s, as Inferred from glacial rebound data (Cathles, 
1975; Peltier, 1976), how could the outer core still be liquid? The existence 
of a liquid outer core in the Earth places a significant constraint on the 
efficacy of subsolidus convective cooling during the thermal history of our 
planet. While we must find a reason why overly efficient mantle convection 
has not frozen the Earth's core, this problem may not exist for one or more of 
the other terrestrial planets if future seismic observations should reveal a 
solid core or if future magnetic observations should confirm the absence of a 
planetary field. 

For the Earth, there are several ways in which core solidification by 
subsolidus convective cooling can be prevented. One way is to have a signifi- 
cant source of radioactive heating in the core. Another way is to prevent' 
convection from reaching the lower mantle for a portion of the Earth's history, 
particularly during the Initial period of cooling after core formation when 
convection should be especially vigorous. This may be accomplished by 
chemically or viscously stratifying the lower mantle or by hypothesizing that 
the lower mantle geotherm is subadiabatic. In view of both the inference 
from glacial rebound data that mantle viscosity is essentially uniform and 
the argument of Sammis et al. (1977) against large viscosity jumps across the 
major mantle phase transitions, a subadiabatic lower mantle accessible only 
by some form of weak penetrative convection may be the more likely explanation. 
The models of Sharpe and Peltier (1978) rely on an assumed subadiabaticlty of 
the lower mantle to prevent core solidification by solid state convective 
cooling. 

In addition to the numerical results mentioned previously, we present a 
simple model of core cooling by subsolidus convection which has an analytic 
solution and which is therefore particularly well-suited for illustrating the 
efficiency of the solid-state cooling mechanism. The model assumes a vigor- 
ously convecting, essentially isothermal mantle with thermal boundary layers 
at the core-mantle and lithosphere-mantle interfaces. The heat flow across 


each of these boundary layers is by conduction; accordingly, the heat flux is 
directly proportional to the temperature difference across the boundary layer 
and inversely proportional to its thickness. The thermal boundary layer 
thickness is in turn assumed to be proportional to a fractional power (about 
1/3) of the Rayleigh number of the mantle. For constant temperatures at the 
upper and lower boundaries of the mantle, a simple energy balance provides a 
differential equation for the cooling history of the mantle. The model in- 
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corporates a temperature-dependent mantle viscosity which has a controlling 
influence on the thermal evolution. The analytic solutions clearly show that 
subsolldus convection can remove all the latent heat of solidification from 
the core of a terrestrial planet on time scales short compared with geologic 
time. 
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A THERMAL HISTORY MODEL FOR THE TERRESTRIAL PLAHETS WITH PARAMETER- 
IZED CONVECTIOH AHD IMPLICATIONS FOR PLANETARY MAGNETISM 
H.N. Sharpe and W.R. Peltier Dept, of Physics, University of Toronto, 
Toronto, Ontario,, Canada M5S 1A7 

With satellite data now availahle for the magnetic field intensities of 
all the terrestrial planets , it is of interest to consider this data from the 
viewpoint of the thermal histories of the planets. One naturally expects a 
relationship between the nature, and even the presence of an externally observ- 
able magnetic field and the thermal state of the interior. The strong, mainly 
dipolar geomagnetic field is believed to be a consequence of fluid motions in 
the outer, electrically conducting iron core. On the other hand, if strong 
magnetic fields were present in the early solar system, then 'cold’ planets 
could have been magnetized, with the result that if portions of these planets 
have never warmed beyond the Curie Point, a present day remanent magnetization 
may be possible. In this paper we propose to study how the thermal histories of 
the terrestrial planets may 'set the stage* for determining the nature of their 
magnetic fields. 

Traditionally, thermal history models for the planets have neglected the 
convective transport of heat. The resultant 'diffusion' models were characteriz- 
ed by very long thermal time constants and generally proposed conduction 
profiles which were highly unstable to convection. Hence these models should not 
be physically realized. In the present study we have developed a technique which 
parameterizes the energetically important aspects of convection without requir- 
ing an explicit solution for the dynamics. The scheme accounts for the enhance- 
ment of the radial flow of energy by modulating the thermal conductivity in the 
diffusion equation with an estimation for the Nusselt number, and is shown to be 
valid provided viscous dissipation is negligible . This latter requirement is 
satisfied for Mars and Mercury, but may be only marginally valid for Venus and 
the Earth, depending on the extent of convection. A power law dependence of the 
Nusselt number upon the Rayleigh number, as deduced from theory and experiment, 
is employed which maximizes the assumed efficiency of convection. 

To commence the thermal history calculations it is necessary to assume a 
specific model for planetary formation. Assuming the moon accreted homogeneous- 
ly, in an isolated environment, and interpreting the lunar crust as a product of 
late-stage accretional melting, it is possible via a simple energy balance 
argument to compute the minimum conversion efficiency and corresponding rate of 
accretion required to initiate melting of the incoming material. For the model 
which we have employed for the accretion dynamics , this provides sufficient 
information to permit the characteristics of the accreting 'cloud' to be found. 
Assuming ^the terrestrial planets formed from 'clouds' of similar characteristics 
it is then possible to construct the initial thermal states of these planets. 
Such calculations show that virtually all of Venus and the Earth would pass 
through an initially molten stage, while only the outer half of Mars and the 
outer third of Mercury would melt. Additional modifications to the initial 
profile due to adiabatic compression, the heat released during iron-silicate 
separation, and finally the action of fluid convection are also considered. 

It is a characteristic feature of the formation models for Venus and the 
Earth that their interiors are totally depleted in radioactive heat sources by 
the end of accretion, all of these sources being concentrated in an outer 
'crustal' layer. This removes an important but poorly constrained parameter. 
While we do not wish to suggest the interiors actually are depleted in heat 
sources , we may consider that for Venus and the Earth at least , we are studying 
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an important end member in a class of possible thermal history models - namely 
the conseq.uences of 'primordial' heat associated solely with the formation 
process itself. 

Earth : The initial thermal structure of the Earth is essentxally Identical 

to the present day state, so that the study of the Earth's thermal 
evolution reduces to a consideration of the conditions \mder which it is 
possible to preserve the initial structure. This in turn leads to the identifi- 
cation of several 'critical' physical parameters. The most important parameter 
is the mean mantle viscosity, which controls the effectiveness of convection in 
cooling the interior and in particular, the molten outer iron core. Assuming 
that the mantle deforms as a power law fluid, it is possible to find a class of 
models which preserves the dominant features of the initial state until today 
and which, in addition, satisfies the mean mantle viscosity as inferred from 
the glacial rehound data(lO^^ poise) , and the observed mean surface heat flux 
(60m¥/m2). It is interesting to note that the corresponding Rayleigh number of 
lo'^ results in an effective thermal time constant for the entire mantle which 
is actually less than that for the moon, considering heat transport in the 
moon as purely diffusive. 

Venus : Applying the identical thermal model to Venus except for specify- 

ing a surface temperature of 750°K instead, of 300°K, we found that 
the higher mean mantle temperatures resulted in an enhanced cooling of the 
interior due to correspondingly lower viscosities. The model suggests that if 
the presently observed high surface temperature of Venus has persisted through 
most of its histoiy, the enhanced vigour of mantle convection would have 
resulted in core solidification about 2xio^ years ago. At present the Venusian 
geotherm would he colder than that for the Earth and the corresponding mean 
surface heat flux would also be diminished by about 30%. This result could 
explain the observation that while the Earth and Venus have similar physical 
properties, the magnetic field of Venus is nevertheless very weak. 

Mars : The thermal evolution of Mars follows a significantly different 

course than that of the larger terrestrial planets. Shortly after 
accretion the small, molten outer iron core solidifies. The next 2x10 years 
are relatively quiescent, characterized by a general warming of the 'primitive' 
inner half of the planet. At about 2x10? years ago this 'primitive' core would 
have overturned, melted and differentiated following the onset of solid-state 
convection eind an assumed subsequent Rayleigh-Taylor instability. The signifi- 
cant heat pulse associated with the formation of the molten iron core would 
have lead to the resumption of mantle convection, although it is unlikely that 
a form of plate tectonics could have ensued. Nevertheless , the surface would 
have undergone significant tensional deformation at this time. The model 
suggests that core solidification should again have occurred about .5x10^ years 
ago and that at present the mantle may he only weakly convecting. These results 
are in qualitative agreement with the observation that the Martian magnetic 
field is extremely weak. 

Mercury : While the early thermal history of Mercury is similar to that of 

Mars, the assumed absence of potassium, in accordance with 
proposed cosmoehemical models, prevents the interior from warming up to where 
the onset of convection would occur. Consequently the inner 2/3 of Mercury 
would remain as 'primitive', undifferentiated material until the present day. 

In this sense its thermal history would most closely resemble that of the moon. 
This raises the interesting question as to the nature of the observed Mercurian 
magnetic field. We propose that the field is a fossil remanence of au early, 
strong solar system magnetic field, and is carried in the outer 300km of the 
planet. The required intensity of magnetization is such that the remanence 
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could easily be carried by material with a composition similar to that which 
has been observed to carry the lunar remanence. 

In this study we have attempted to understand how the thermal history of 
a planet may determine the nature and even the existence of its magnetic field. 
While the size of the planet is the single most Important parameter, there are 
several other parameters which also enter the evolutionary study in an import- 
ant way. These include the surface temperature, the concentration and distribu- 
tion of radioactive heat sources and the presence and vigour of convection. 

This is only a preliminary investigation and rests on several ad hoc assumptions 
the most important of which surround the formation model. However because of 
the dominating influence which solid-state convection seems to have had on the 
thermal histories of Mars, Venus and the Earth, we feel that the present study 
is a significant step forward in imder standing the energetics of planetary 
evolution. 
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ON THE DETERMINATION OF THE HERMAEAN MAGNETIC MOMENT. 

James A. Slavin and Robert E. Holzer, Department of Earth and Space Sciences and 
the Institute of Geophysics and Planetary Physics, University of California at 
Los Angeles, Los Angeles, California 90024. 

Mariner 10 passed through the near tail portion of Mercury's magnetosphere 
on March 29, 1974 and March 16, 1975. In addition to these encounters, termed 
MI and Mill respectively, a second pass on September 21, 1974 did not bring the 
spacecraft near enough to observe either the magnetosphere or bow shock. Sub- 
sequent studies (7,8,12,13,14,15) have found that magnetic field and particle 
observations to be consistent with Mercury possessing a permanent terrestrial- 
type magnetosphere as modified by both this planet's proximity to the sun and 
large size relative to its magnetospheric cavity. These two factors result 
in the magnetospheric current systems being nearly tantamount to the planetary 
magnetic moment as source of the observed magnetic field during part of the 
Mariner 10 encounters. 

As noted by Ness (8,10) a unique determination of the intrinsic magnetic 
moment by means of spherical hamionic representations of both the internal and 
external scalar magnetic potentials would require the determination of two or 
more components of the field on a surface enclosing the planet. However, the 
Mariner 10 mission returned only *^30 min of magnetic field data from the night- 
side magnetosphere so that it becomes necessary in modeling the observations to 
make assumptions regarding the position and intensity of both the magnetopause 
and neutral sheet currents. This has been done in various ways by Ness et al 
(8,9 ), Whang (17), Jackson and Beard (6), and Ng and Beard (11) with different 
results. When the intrinsic field is assumed dipolar, the values Inferred for 
the magnetic moment range from 2.9-5.2xlO^^G-cm^ depending upon the portion of 
the data considered and the method used. Dipole moments of 2.4-2.9xlO^^G-cm^ 
are found when both MI and Mill are modeled together and the intrinsic field is 
allowed to possess higher order moments and/or offset positions within the 
planet (17,6,11). However, there are serious doubts concerning the physical 
reality of higher order moments inferred from this spatially and temporally 
limited, set of observations (8,10,11). 

An alternate approach to determining the Hermaean magnetic moment is to 
rnfer values of the solar wind stand-off distance, rs, from the positions of 
the Mariner 10 bow shock and magnetopause encounters. The bow shock shape and 
position are determined by the forward portion of the magnetopause surface which 
is expected to have the same blunt shape at both the earth and Mercury. Hence, 
with a knowledge of the shape of the terrestrial bow shock (1) and the ratio of 
the distance to the nose of the shock to rs (3) the solar wind stand-off dis- 
tance may be calculated at the time of each of the 4 bow shock crossings. The 
earth's magnetopause out to sun-earth-satellite angles of MIO^ is known to be 
elliptical with the earth at one focus (4). If the solar wind is indeed de- 
flected around Mercury by an intrinsic dipole moment as found by Ness et al (9), 
then the shape of the magnetopaqse'out to a sun-Mercury-satellite angle where 
the tail field begins to dominate will also be elliptical. Because the Her- 
maean magnetotail current density is greater than in that of the earth (8), 
this angle will be less than 110°. For the purposes of this study the Hermaean 
magnetopause has been assumed elliptical out to sms angles of 'vlOQO so that the 
values of rs were inferred from the two outbound crossings at 83o and 101°, but 
not for the inbound crossings at 1350 and 1210. In this way the solar wind 
stand-off distance has been determined on 6 occasions during MI and Mill when 
the satellite encountered these boundaries. 

The stand-off distances found above are a function of the solar wind dy- 
namic pressure, the amount of magnetic flux transferred to the tail by erosion, 
and the intrinsic planetary magnetic monients. Due to an instrument failing to 



92 


HERMAEAN MAGNETIC MOMENT 
Slavin, J.A., and Holzer,R.E. 


function after launch, no in situ solar wind ion observations were made by Mari- 
ner 10. Solar wind dynamic pressure during MI and Mill may be inferred from 
the electron plasma observations (12,13), the observed magnetic field with 
assumptions to determine the state of compression, and solar wind observations 
made near the earth and scaled back to Mercury at the time of the encounters 
(16), In this study the pressure during MI and Mill is assumed to be 1.3x10“^ 
and 2.6x1 0"^dynes/cm^±30% respectively, on the basis of the magnetic field ob- 
servations (16). Whereas an increase in solar wind pressure produces a com- 
pression of the dayside magnetosphere, erosion transfers magnetic flux to the 
magnetotail and alters the magnetospheric currents in such a way to produce a 
contraction of the dayside magnetosphere. Slavin and Holzer (16) have scaled 
the effects of erosion on a-t the earth to Mercury and found that the magneto- 
pause may be found within 0.2Rni of the Hermaean surface a significant portion 
of the time. In addition, this increase in tail flux enhances the drag on the 
tail by the solar wind which is then balanced by the neutral sheet moving 
closer to the earth (e.g. 2). This results in decreases in the magnetic field 
strength of the near tail regions as evidenced during the outbound portion of 
MI. The stand-off distance at which a balance between the solar wind and the 
magnetic pressure is obtained is also dependent upon the strength of the 
planetary magnetic moments. As the models of Whang (17) and Jackson and Beard 
(6) which attribute significant portions of the fields observed by Mariner 10 
to quadrupole and/or octupole moments are not consistent with the solar wind 
pressures and values of rg found above, it is assumed that the Hermaean mag- 
netic field is dipolar as found ^y Ness et al (9). Pressure balance at the 
stagnation point then requires Mo=2TTpv^Rs®k/f^ where Mo is the dipole moment, 
pv^ is the solar wind dynamic pressure, Rs is the stand-off distance in the 
absence of erosion (i.e. ground state magnetosphere), and k/f^ is a constant. 
With this relationship all 6 values of rs may be scaled to a single solar wind 
pressure as has been done in the table and the dipole moment calculated. Since 
the IMF was observed to be northward only prior to the Mill bow shocks, the 
stand-off distances inferred from these crossings should correspond approxi- 
mately to uneroded stand-off distance, R$. 

Crossing 1*5 As shown the mean stand-off dis- 

tances for MI and Mill are 1.5 

MI Inbound BS 1 .33Rm±14% 1 ‘9Rm respectively. The smaller 

Outbound MP 1.76Rfu±13% values of-rg inferred for MI are 

Outbound BS 1 .43Rm±20% consistent with both the south- 

MIII Inbound BS 2.08Rm±ll% ward IMF and substorm signatures 

Outbound MP 1.62Rtn±15% (15) observed during this en- 

Outbound BS 1.98Rm±17% counter and not during Mill. In 

(All scaled to 6.0xl0-8dynes/cm^) addition, this evidence for mag- 
netic flux transfer during MI 

explains the finding of smaller dipole moments for MI relative to Mill from 
modeling the observations with the magnetopause and neutral sheet currents 
assumed to have the same position and intensity during both encounters. The 
uneroded stand-off distance inferred from the boundary crossings is the larg- 
est value in the table, 2.1±0.2Rm. With this stand-off distance and the pre- 
viously stated assumptions and uncertainties a Hermaean dipole moment of 
6±2xlO^^G-cm^ is found from the pressure equilibrium condition. This value is 
in substantial agreement with the values of 4.7-5.2x10^^ arrived at by Ness 
et al (9). However, final determination of Mercury's magnetic moment must 
await an arbiter mission of sufficient duration to study the time dependent 
contributions of magnetospheric currents and possibly planetary induction cur- 
rents (5). 
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The main properties of Jupiter's intrinsic magnetic field 
are known from radio astronomy observations and from the Pioneer 
10, 11 measurements. Qualitatively, the field is that of a 
reasonably well-centered dipole. A good quantitative^approxima- 
tion is provided by a dipole with moment 4.2 Gauss rj > tilted 
10° relative to Jupiter's rotation axis and displacea from the 
planet's center by 0.1 rj in an equatorial direction- These 
parameters, as well as the longitude of the magnetic pole, are 
consistent with both the remote and the in situ observations. 


However, both sets of observations provide clear evidence that 
the jovian field is actually more complex than that of a simple 
dipole. Spherical harmonic analyses of the spacecraft data have 
been carried out in the region betxvreen 1.6 and 8 rj and reveal 
significant contributions from a centered quadrupole and octupole. 
Non- dipole structure is revealed in the radio astronomy data 
through asymmetric variations in brightness, in the position 
angle of linear polarization, etc. In principle, it should be 
possible to obtain a better description of Jupiter's field from 
both sets of data than from either set taken alone. Furthermore, 


the radio astronomy data acquired over many years also make it^ 
possible to look for secular variations. By a fortunate coinci- 
dence, Sam Gulkis and I have just completed a review of space- 


craft and radio studies of Jupiter's field for a forthcoming 
volume of Annual Review of Earth and Planetary Sciences. This 


talk will be based to a large extent on that review. 
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Ever since the two Pioneer spacecraft provided us with 
spectacular data about the Jovian magnetosphere (1) much effort 
has been devoted to the search for magnetic fields of the other 
giant planets. Although, the larger distances and weaker fields 
of these planets prevent obtaining data of an accuracy and de- 
tail comparable to what we knew about the Jovian field before the 
Pioneer flights, it seems possible to draw certain general con- 
clusions and make comparisons apart from the tentative linear re- 
lationship between magnetic moments and angular momenta of all 
the planets. In this paper the internal fields and their origin 
will be discussed with only parenthetical mention of the external 
fields and of magnetospheres although these provide the main 
clues for the study of the former. The analysis of the origin of 
the magnetic fields of the giant planets is somewhat facilitated 
by the fact that, in contrast to the terrestrial planets, the 
chemistry and structure of the giant planets is relatively . 
simple. Unfortunately the various parameters which enter into 
the necessary conditions (2,3) (no sufficient conditions have as 
yet been unambiguously derived) for the generation of an inter- 
nal magnetic field are so sensitive to subtle details of the mo- 
dels that no rigorous predictions can be made at this time. 

The generation of a magnetic field by the hydromagnetic 
dynamo mechanism (2) requires a sufficiently high magnetic Rey- 
nold's number in which the electrical conductivity of the liquid 
plays the major role. For this reason it was generally accepted 
that on Jupiter the field is generated (A-) only in the metallic 
hydrogen interior ( r-^0.75 Rj) . It turns out, however, that the 
high pressures and temperatures prevailing on this planet lead 
to the presence of metallic conductivity also in the inner part 
(0*75 < 0.92 Rj) of the H£ mantle (5) thereby increasing the 

high conductivity volume by about a factor of two and bringing 
it closer to the surface. The latter aspect seems important 
from the point of view of recent theories ( 6 ) which indicate 
that the field is generated primarily in the outer parts of the 
high conductivity region. It may also help to understand the un- 
usually high quadrupole and octopole components of the field 
just above the surface ( 5 ). 

In contrast to the early negative results recent terrest- 
rial observations ( 7 ) indicate that there is indeed a magnetic 
field on Saturn but that the rings supress the formation of a 
strong radiation belt ( 1 ) so that only hectaraetric pulses and no 
decimetric radiation are observed. Pour recent models of Saturn 
have been examined ( 8 ) and it turns out that the volume of high 
conductivity H 2 is 3 to 6 times larger than that of the metallic 
hydrogen so that it is in this region (0.4 <r <0.7 Rg) that the 
field is generated. The top of this region is, however, more 
than twice as deep below the surface of the planet as on Jupiter 
so that the higher components of the field are probably weak. 
There is also the possibility that the magnetic field of Saturn 
is generated by the thermoelectric effect or by the precession 
currents induced by Tethys ( 3 ). 
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One treads on much shakier ground when the magnetic fields 
of Uranus and Neptune are considered. Although now we know that 
Uranus has a field but probably no internal heat while the oppo- 
site is true for Neptune ( 9 ) one cannot claim that these facts 
are in agreement or not with what we know about the interiors of 
the tvfo planets. The absence of a heat source on Uranus makes 
the hydroraagnetic dynamo mechanism somewhat questionable. The 
thermoelectric mechanism in the layer is a- possibility and one 
should not forget either, that, in contrast to the satellites, 
Uranian rings (especially uhe biggest £-ring) have probably sub- 
stantial inclinations to the equatorial plane and could induce 
the Malkus precessional mechanism ( 3 )* In spite of huge dif- 
ferences between the models of these planets as proposed by 
various authors (10) one can conclude that, unlike Jupiter and 
Saturn, neither one has much, if any, metallic hydrogen and that 
both have large cores made of ‘’rocks" and "ices" although metal- 
lic ammonia has also been suggested (10), On the whole the pres - 
sure and density in Uranus are lower than those in Neptune and de- 
pending upon the equations of state used one might expect much of 
the cores to be solid. The magnetic could be thus generated in 
the liquid H 2 layers and in the liquid parts, if any, of the 
cores. Neptune, having a presumably colder core with lower elec- 
trical conductivity and a thinner layer would be less likely 

to sustain a hydromagnetic dynamo than Uranus. These conclusions 
are clearly based on hindsight rather than on foresight. Best 
estimates indicate that the magnetic moments of Jupiter, Saturn 
and Uranus are in the ration 1000:100:1. 

Finally one should mention that while the Jovian magneto- 
sphere, or rather the curious magnetodisc, has been investigated 
in considerable detail (1) all one can do about Saturn is to 
scale down the Jovian magnetosphere (12) and take into account 
the perturbing role of the rings. Uranus with its unique enor- 
mous tilt of the axis must have an interesting magnetosphere and 
magnetotail (1,8). In fact one may suspect that if there is a 
substantial 'radiation belt then this belt and the radiation 
reaching Earth will show seasonal variation with a period of 
84 years. One does nol; expect significant interaction of this 
magnetosphere with the rings because their width is several 
hundred times lower than that of the rings of Saturn (13). 
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LUNAR THERMAL HISTORY VERSUS LUNAR MAGNETISM, OR 
HOW TO BOX IN A HYPOTHESIS AND BUILD A COFFIN; Sean C. Solomon, 
Dept, of Earth and Planetary Scxences, Massachusetts Institute 
of Technology, Cambridge, MA 02139. 

All theories for internally generated magnetic fields for 
the Moon involve explicitly or implicitly some assumptions about 
the lunar thermal history. It is of interest as a test of such 
theories to explore fully the range of possible thermal histories 
for the Moon consistent with geological, geophysical, and geo- 
chemical data independent of the evidence from lunar paleomag- 
netism. The range of possible thermal models on these grounds 
is surprisingly narrow, and in particular poses severe problems 
for an early (first 0.5 b.y.) lunar core dynamo. 

The essential constraints of a thermal history model are 
(a) that the current heat flow [1] and relative abundances of 
radioactive heat sources be satisfied, (b) that a 70 km thick, 
feldspar-rich crust have formed by crystal-liquid differentiation 
[2] early (first 0.2 b.y.) in lunar history [3,4], (c) that a 
distinct lunar core be no larger than 400 km radius at present 
[5,6], (d) that major surface volcanism have lasted until at 
least 3.0 b.y. ago [7] , (e) that the thermal stress in the lunar 
lithosphere never have accumulated to high enough levels (>1 
kbar) to cause global-scale failure [8] , but (f) that the global 
thermal stress provided an important modulation to the local 
stress due to mare basalt loading in mare basins so as to affect 
the dominance of local extensional or compressive failure [9] . 

Constraint (b) excludes models with a cold initial state 
in favor of models with at least the outer few hundred kilometers 
of the Moon initially at or near melting temperatures . A neces- 
sary condition for (d) , that partial melting persisted in the 
mantle until that time, is satisfied for most models consistent 
with (a) and (b) . Constraint (c) dictates that core formation 
on the Moon is at most a very modest thermal event, releasing 
10^^ J of potential energy as heat, or only a 12 °C equivalent 
mean temperature rise [10] . 

Constraints (e) and (f) are quite powerful. The absence of 
global-scale compressional or extensional tectonics on the Moon 
greatly restricts possible thermal models to those with nearly 
constant volxmies and small accumulated iithospheric thermal 
stresses since the time of heavy bombardment. In particular, 
models with the outer 300 ±100 km of the Moon initially melted 
and the deep interior of the Moon initially cold are possible 
[11] , but models with completely molten initial states [e.g. , 12] 
may be excluded. This conclusion includes consideration of heat 
transfer by subsolidus convection [13] and of models of the type 
advocated by Tozer [14] . The requirement that the change in the 
stress regime near mare basins from one dominated by graben 
formation to one dominated by mare ridge formation at about 3 . 6 
b.y. ago [15] be due to a change in sign of incremental global 
scale thermal stress [9] restricts possible thermal histories 
still further (figure 1) . 
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Thus whxle metal-silicate fractionation in the outer 300 
km may have occurred quite early on the Moon, formation of a 
central core close to final size should have awaited warming of 
the deep interior sufficient for relative flow of mantle and 
core material. Such a time should be close to the time of peak 
lunar volume, or about 1 b.y. after origin [9] . Models for 
magnetic field generation involving a fully developed fluid core 
with an active dynamo 4 b.y. ago [12] are thus unlikely. Models 
with a permanently magnetized, initially cold interior [16, 17] 
are acceptable on thermal and geological grounds . 
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Fig. 1. Summary of temporal relationships among lunar thermal 
history/ global expansion and contraction, and the tectonic and 
volcanic history of a typical mare basin, Serenitatis [9] . Iso- 
therms are shown versus depth and time for a thermal history 
model with initial melting to 300 km depth and an initial central 
temperature of 300°C. Regions of partial melt are shown as 
stippled. 



THE MAGNETIC FIELD OF THE EARLY SOLAR SYSTEM. C.P. Sonett, Lunar and 
Planetary Laboratory, Univ. of Arizona, Tucson, AZ 85721 (*also Dept, of 
Planetary Sciences). 

The strongest evidence of the presence of a magnetic field during the 
formative phase of the solar system comes from observation of the presence of 
ancient thermoremanent magnetization possibly joined with chemical remanent in 
meteorites [1,2]. The requirement for redistribution of angular momentum in 
the proto-solar nebula also infers the existence of a solar spin braking mech- 
anism, which can reasonably be driven by magnetic torques [3,4]. Magnetization 
of the Moon and Mercury, the most critical planetary examples provides less 
conceptual connection to the early field problem because, for example, the 
maria basalt magnetic fields were imprinted between 3.9 and 3.2 aeons ago, 
rather later than primordial. The source for Mercury's magnetic field is 
generally thought to reside in a dynamo; a less favored source is thermorema- 
nence of the crust which could be primordial [5]. 

Recent work by Lanoix [6] on chondrules from Allende shows evidence 

for a primordial magnetic field strength as large as 16 oersted, a figure 
substantially higher than previous values of 0.18-2 oersted. This high value 
may relflect a different magnetic field at the time of acquisition of TRM by 
Allende chondrules than that for the matrix. The primordial magnetic field 
appears most likely to reflect the presence of either a magnetic field extend- 
ed from the early Sun or a nebular dynamo [2]. 
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CONVECTION-DRIVEN DYNAMOS A.M, Soward, School of Mathematics, The 
University, Newcastle upon Tyne, U.K. . 

The fluid motions: responsible for driving the geodynamo are generally 
beleived to be the result of convection, thermal or otherwise. The detailed 
nature of the processes involved is still unknown. Nevertheless owing to its 
relative simplicity theoretical dynamo models have been developed on the basis 
of thermal convection. 

A model, which incorporates the main features of the Earth's liquid core 
is as follows. A self-gravitating, Boussinesq fluid is confined in a 
spherical container and rotates rapidly about an axis with angular velocity 
_TL . The fluid contains ,a uniform distribution of heat sources, which 
leads to an unstable radial temperature gradient p . A dimensionless 
measure of 3 is the modified Rayleigh number 

where c< is the coefficient of expansion, ^ is the acceleration due to 
gravity, L is the core radius, is the thermal conductivity. The critical 
value of the Rayleigh number at which the onset of convection occurs 

depends on the Ekman number E and the Prandtl number <r , where 

and y is the viscosity. A detailed numerical study of this model has been 
made by Roberts (1968). IVhen the fluid is electrically conducting, the 
stability characteristics of the system are altered in the presence of mag- 
netic field, Eltayeb and Kumar (1977) have considered the special case of an 
azimuthal magnetic field, which increases linearly with distance from the 
rotation axis. They find that depends critically upon the sizes of 

r = CS' ^ 

where B/A is the magnetic field gradient, &- is the electrical conductivity, 
is the fluid density and |>- is the magnetic permeability. 

For geophysically relevant values of the parameters, is small and for 
a magnetic field of about 30 gauss P is of order unity. The corresponding 
value of R^. is of order one also and in this case the primary force balance 
is between the Coriolis and Lorentz forces. As the magnetic field is 
decreased ^ increases and when the field vanishes is proportional to & . 

According to the model convection occurs most readily when the magnetic field 
is large enough to release the rotational constraint but not so large that 
strong Lorentz forces inhibit convection (i.e. P = OQ.))* The result 
suggests that the geodynamo may operate in this parameter range but owing to 
the mathematical complexity no theoretical model has been developed. Instead 
attention has been directed mainly towards quasikinematic models for which 
the character of the motion is largely (but not completely) uninfluenced by 
the Lorentz force. 

Busse (1975) has proposed a weak field model C P * 1 of the geodynamo, 
which he developed in an annulus geometry. A dynamo model, which is similar 
in many respects, was also developed by Childress and Soward (1972) in the 
simpler plane layer geometry. For this model Soward (1974) found stable 
hydromagnetic dynamos with very weak magnetic fields, P = 0(E“J-). By con- 
trast Childress (1976) discovered that, when the magnetic field is stronger, 
r = 0(E"2/3)^ -the djmamo process is unstable. This is simply because 
increasing the field strength releases the rotational constraints. On the 
other hand, Basse's model is stable with substantially larger magnetic fields. 
This is explained by the additional constraints imposed by the sloping 
boundaries of the annulus. Nevertheless, for small but moderate field 
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strengths there remains the danger o£ the instability isolated by Childress. 
l£ this occurs, it is reasonable to suppose that th.e dynamo will equilibrate 
with, an order one value of P . 
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CRUSTAL THERMOREMANENCE AND THE ORIGIN OF MERCURY’S PRESENT 
MAGNETIC FIELD; L. J. Srnkaj Lunar and Planetary Institute, 
Houston, Texas 77058. 

Mercury is the only terrestrial planet in addition to 
the Earth which is known unequivocally to have an intrinsic global 
magnetic field. If the origin of this field could be shown to be 
due to dynamo action in a mercurian liquid core, the planet could 
serve as a crucial test for dynamo theories which have been 
developed primarily to understand geomagnetism. 

The discovery from the Apollo program that lunar crustal 
rocks have magnetizations which date from the early history of 
the Moon, and the growing body of data on meteoritic magnetism, 
suggest that many solid objects in space are remanently magne- 
tized. If this remanence is ordered, it may be possible to 
account for substantial planetary magnetic fields without invok- 
ing currently active dynamo processes. Models of this kind have 
been proposed for Mercury by Stephenson (1976) who studied 
crustal remanence acquired in a mercurian dipole field produced 
by a once-active dynamo, and by Sharpe and Strangway (1976) who 
considered the possibility of fossil magnetization in Mercury 
acquired from a strong primordial solar magnetic field. The 
resolution of this present dynamo vs. remanent field question is 
important for the further understanding of dynamo theory and 
(through the physical link of convective processes) the thermal 
evolution of Mercury and the other terrestrial planets. 

The acquisition of global thermoremanent magnetization (TRM) 
in Mercury’s lithosphere is studied here using a planetary model 
which incorporates both the thermal evolution of the planet and 
the possible character! sti cs of a mercurian dynamo. The thermal 
evolution model ofSolomon (1 977) is used to track the evolution 
of the Curie isotherm depth (for Fe) with time. Mercury's 
lithosphere is assumed to consist of an unspecified paramagnetic 
mineral matrix containing dispersed, equidimensional Fe inclu- 
sions which are multidomain and are present in volume concentra- 
tions f of lO”^ 5 f < 10%. A core-generated dipole field is 
assumed to have existed, with the general form; 

P(t) = Pj^t^ cos (2TTt/T) exp (-t/x) 

where T is the reversal period, x is the decay constant for the 
field envelope at late times, and the factor t<^ models the rise 
of the dipole moment during core formation. As a test case, the 
parameters n = 2.47, x = 5 x 10® yr, and 10® 5 T S 10® yr were 
used. These result in a dipole field which peaks at 700 myr after 
the "origin" of the planet, after Mercury's core has fully formed 
(Solomon, 1977). The field is "turned off" 100 myr ago, after 
which time the interior changes very little. 

The results of this calculation are given in Fig. 1, These 
show that if the ancient magnetic dipole moment of Mercury's 
core-generated field did not exceed ~ 5 x 10^® G cm®, if the 
dipole field reversed at least once every 10 myr, and the volume 
fraction of Iron in the planet above the Curie isotherm depth is 
less than 10% by volume, then the present upper limit on Mercury’s 
dipole moment of ~ 4.8 x 10^^ G cm® (Ness, 1978) is too large to 
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be attributed to crustal TRM acquired in an ancient mercurian 
dipole field. On this basis, a currently active dynamo is favored 
as the source for Mercury's field. On the basis of the lunar 
evidence, primordial magnetization in Mercury due to an early, 
strong solar field is unlikely to be the source of its present 
magnetic dipole moment. 
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A STUDY OF THE RIMA SIRSALIS LUNAR MAGNETIC ANOMALY; 

L. 0. Srnka, J. L. Hoyt, J. V. S. Harvey, Lunar and Planetary 
Institute, Houston, Texas 77058, and 0. E. McCoy, NASA/Oohnson 
Space Center, Houston, Texas 77058. 

The discovery of a magnetic anomaly associated with the lunar 
grabenlike nonsinuous rille Rima Sirsalis by Anderson et aZ. (1976, 
1977) is the first positive identification of a lunar surface 
feature which has a local magnetic field. In this paper we study 
various source models for the anomaly which we believe are 
compatible with the geology of the rille and its surroundings. 

The source of this linear anomaly has been modelled using 
exact solutions for the scalar potential and vector field inten- 
sities due to arbitrary uniform remanent magnetization in a finite 
rectangular prism. The earlier work of Anderson et al, (1977) 
used a line dipole model in the form of a circular cylinder of 
infinite length, which has the disadvantage of poorly represent- 
ing the grabenlike structure of the rille and hence the higher 
order magnetic moments of such an anomaly (i.e., fields due to 
sharp corners) near the surface. As a result, the required source 
magnetization has been underestimated. 

Three classes of source models were studied: 1) rectangular 

prisms with magnetizations limited to those found in lunar sam- 
ples; 2) a gap in a uniformly magnetized layer; and 3) a system 
of narrow, parallel dikes intruding the rille and its surround- 
ings. The following general results were obtained: if the field 

source is confined to the width of the rille and the magnetization 
does not exceed that found in lunar samples, the body responsible 
for the anomaly must have a vertical extent of tens of kilometers. 
More specifically, if the field source is a number of thin (100 - 
200 m wide) dikes beneath the rille, the intruded material must 
have a uniform downward component of remance of at least 10"® emu 
g’^ over at least 30 km of depth. We also find that all reason- 
able gap models for this anomaly require surface magnetizations 
of at least 10"® emu g"^ if the gap dimensions are those of the 
rille. All of the source models require a total magnetic moment 
(all orders) of 3 - 8 x 10^® G cm® for this anomaly if the 
Apollo 16 subsatellite magnetometer data (B^ » - O.ly at 155 km 

altitude) is to be matched. This moment is comparable to the 
magnetic dipole moment derived for the Van de Graaff lunar 
magnetic anomaly by Russell et al. (1975). 

If the model magnetization is constrained to those character- 
istic of the returned lunar samples, we favor a massive dike 
swarm as the principal source of the Rima Sirsalis anomaly. A 
system of 20 parallel dikes, each 200 m wide, lying within 5 km 
of the axis of the rille, and extending from near the surface to 
a depth of 50 km or so (the base of the crust) would produce the 
required field at 155 km altitude if the dikes have a downward 
magnetization of 2 - 3 x lO”® emu g~^. The surface fields in such 
a model are irregular, with vertical component peaks of about 
350y over the individual dikes, which is consistent with the 
electron reflection data (Anderson et al.^ 1977), However, this 
remanence is considerably greater than that found in lunar basalts 
(by a factor of 3 or more) so that it is necessary to invoke 
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unusual magnetic properties for this region of the moon. 
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THEORETICAL CONSTRAINTS ON THE MAGNETTSM OF GASEOUS AND ICY PLANETS . 
D.J. Stevenson, Etept .Earth and Space Sciences, U.C.L.A., Cal. 90024, USA 

This paper is in part a conpilation of the present theoretical under- 
standing of the interiors of gaseous and icy planets, and how this knowledge 
relates to the existence of a global planetary magnetic field (c.f. Stevenson, 
1978a). However, the enphasis is on three new results: (i) A global 
constraint on the Jovian dynamo which demonstrates its turbulent character, 
(li) An estimate for the outer boundary of dynamo cessation in these planets 
from a consideration of dynamo generation in a strong conductivity gradient, 
(lii) A consideration of hi^ pressure water and/or ammonia, indicating that 
dynamo generation in icy bodies (Uranus, Neptune; perhaps even Garymede, 

Titan, Pluto ...?) is conceivable. 

The Jovian Eynamo ; The total internal heat flux of Jupiter is Q 5xl024erg/ 
sec. In a dynamo driven by thermal convection, the Ohmic dissipation should 
be of order Ocd/Hn] where Qc is the core flux, d is the thickness of the 
dynamo region and Ht is the temperature scaleheight (Hewitt et al,, 1975)- 
The dissipation could be even larger if the convection is non-thermal, but 
the formation of helium rain, which has been suggested as an energy source 
for Jupiter and Saturn (Salpeter, 1973; Stevenson and Salpeter, 1977b) is 
still a predominantly thermal source because of the small size (''^Icm) of the 
raindrops. Since Qc is almost as large as Q, and d Ht, the Ohmic 
dissipation is conparable to the total heat flux. However, a global field of 
10 Gauss in Jupiter would only dissipate of order lOl^erg/sec for a magnetic 
dlffusivi-ty of 2xl02cm2/sec (Stevenson and Salpeter, 1977a) . It is hi^ly 
improbable that the eight orders of magnitude discrepancy is due to an 
enormous (lOSoauss) toroidal field in Jupiter. Instead, the discrepancy is 
best explained by a turbulent dynamo in which the dissipation is primarily at 
snail lengthscales (i.e. kilometers). The known properties of Jupiter, 
together with a model of turbulent theimal convection (Stevenson, 1978b) 
suggest a magnetic energy spectrum consistent with that expected for helical 
turbulence (Pouquet et al,, 1976). 

Eynamo Generation in a Conductivity Gradient ; These planets are characterized 
by no rigid outer boundary at which dynamo generation ceases (such as the core 
-mantle boundary of the Earth) but by a gradual cessation of dynamc generation 
because of a nagnetlc diffusivity which increases rapidly in the radial 
direction. An exanple of this is molecular hydrogen at megabar pressures, 
where the band gap between valence and conduction electron states becomes 
sufficiently smll for the, electrical conductivity to be almost metallic 
(Smoluchowski, 1975) • A ‘consideration of oi-effect dynamos suggests that the 
outer boundary of dynamo generation can be defined crudely in terms of the 
region in which the toroidal field becomes small relative to the poloidal 
field. This region is characterized by V5-0/x'^^l 'where V is the convective 
velocity, A is the nagnetlc diffusivity and = (dS,nA/(3p)-l is the 
conductivity scalehei^t. Since Ha is much less than the planetary radius, 
dynamo generation ceases at greater depths than is sometimes supposed. The 
outer boundary of dynamo generation is about 0 . 8 of the radius in Jupiter and 
0.5 of the radius in Saturn. These results also have implications for 
possible magnetic secular variations in these planets. 

Eynamo Generation in Icy Bodies ; Infrared observations (Lowensteln et al., 
1977) suggest that Neptune has an intejcnal energy source but Uranus does not. 
Actually, this is misleading because Uranus is closer to the Sun so its 
smaller internal energy source is more difficult to detect (Hubbard, 1977) • 

In fact, the deep interiors of both planets should be warm ('u3000°K) , fluid 
and convective (Stevenson, 1978a) . Water does not become metallic until 
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megabar pressures, but it is already an excellent ionic conductor (with a 
mgnetic diffusivity of order lO^em^sec”^) at about 200kbar (Hamarm and 
Linton, 1966), Addition of NH 3 aids the conductivity further. The 
possibility of metallic ammonium seems unlikely (Stevenson, 1975) but an 
adequate magnetic Reynolds number for dynamo generation seems likely 
(Stevenson, 1978b). The icy satellites and Pluto may also have convective, 
fluid H 2 O-M 3 layers, but the pressures are lower and dynamo generation is 
marginal at best. 
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THE ROIE OP ROT/ITION IN PLANETARY DYNAI^DS 
D.J.Stevensonj Dept.Earth and Space Sciences, U.C.L.A., Cal. 9002^, USA 

Dynamo generation is usually discussed in terms of a magnetic Reynolds 
number Since rotation evidently plays a crucial role, the value of the 
Rossby nuntoer Rq (representing the ratio of inertial to Coriolis effects) 
should also be considered. The basic conjecture in this work is that and 
^ t ogether determine necessary "^d""suf^ci^ conditions for dynamo ekist^ce . 
Rotation also determines the symmetry of the resulting field and its magnitude. 
To support these assertions, a turbifLent thermal convection model is 
constructed in which nominal values of and Rq are directly determined by Pc, 
the convective heat flux. (The basic principles are unaltered if a 
compositional gradient drives the convection.) 

Relevant and Irrelevant Parameters : The flux Pc, convecting fluid depth d, 

rotation rate rragnetic diffuslvity X, boundary conditions and thermodynamic 
properties are sufficient to specify the system. Viscosity v and thermal 
diffusivi-t^ K are excluded since Vd»v,K for any realistic convective velocity 
V. (In the Earth's core, Vd 'V' lO^cm^sec"^, v 'v lO'^cm^sec-l, k 10“lcm2sec~l) . 
Turbulent convection is Inevitable, and any dynamo model which explicitly 
depends on v or fc is either unrealistic or relies crucially on the effects of 
boundary layers. 

The Convective Model (Plasar and Gierasch, 1978; Stevenson, 1978) considers 
the growthrates of linear, Boussinesq convective modes in a planar geometry 
for a superadiabatic fluid with inposed uniform rotation and ma^ietic field. 
Velocity amplitudes are estimated by assuming that these modes are stabilized 
by non-linearities, such as shear instabilities, which cascade the kinetic and 
thermal energies into smaller scale motions. The dominant mode is assumed to 
be that which minimizes the teirperature gradient for a given heat flux. Por 
no field or retation, the predictions are equivalent to those of mixing length 
theory, and agree with observed convective velocities in stellar atmospheres 
(Bdhm-Vi tense, 1977) and thermally-driven turbulence in the earth's 
atmosphere (Monin and Yaglom, 197D • For rapid rotation and no field, the 
preferred modes are rolls aligned along the horizontal component of Q (c.f . 
Busse and Cuong, 1977) • The predicted convective velocity is 

V - (Fcd/pHT)'/^(lfizid)-'/= (1) 

where ^2 is the vertical component of fi, p is the fluid density and Hfp is the 
temperature scaleheight. The wavevector of the preferred mode is essentially 
azimuthal and has magnitude -v RQ-id“T where Rq = Provided Rq < 1, 

the flow possesses a non-zero horizontally averaged helicity < V. V x V > and 
an a-effect is predicted which can transform 'toroidal field to jjoloidal 
field Bj, (c.f. Childress, 1977): 

(^ - = Y X (oBp) (2) 

2 . ^ 

^ ^ Dp^zd (3) 

2Xjf2^1 


where z is the vertical coordinate . This a-effect does not transform Bp back 
to Bm, but this can be done by: 

The Thermal Wind . The convective model predicts a latitude dependent 
tenperature gradient for a given heat flux. For realistic boundary conditions, 
a horizontal temperature gradient is unavoidable. As in the Earth's 
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ROTimON AND PLANETARY DYNAMOS 
Stevenson, D.J. 



Plg.l : Ihe Rossby number - magnetic Reynolds number plane for convective 

dynamos, showing the curve which separates dynamo from non-dynamo. Positions 
of Earth, Jupiter, Saturn, Neptune and the Sun are reasonably well established. 
Ihe dashed error bar for Uranus extends into the non-dynamo region because of 
the large uncertainty in the heat flux. Mercury, Venus and Dfers are shown 
bracketed at the positions they would occupy for Pe “v 10 c.g.s. In fact, Pq 
may be zero in all three . Pluto and the icy satellites lie at the boundary 
and are tests of the theory. 
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atmosphere, a theimal wind is expected (c,f. Eady, 19^9). Ihis provides the 
''oj-«ffect” to conplete an oiw dynamo. 

The Zero-Pield ato-Eynamo ; Ihe model predicts a thermal wind U,^ of order V 
and a~YFL^% where ^ Eynamo onset therefore requires (Moffatt, 

■ U^d , 

(-^) (x> ^ 100 (4) 


(the precise numerical value is unimportant) . Necessary and sufficient 
conditions for dynamo generation are thus % > IORq-Vs and Rq < 1. This is 
illustrated in Pig.l and can be calculated for planets if Fq is known. 
Finite-Pleld Equilibration : If lORo^Ve < % < 10Ro“yih then the field 

equilibrates by reducing a and Ui]^ so as to satisfy eq.(4-). (Tie finite field 
enhances convective efficiency, and reduces convective velocity and 
temperature gradients ) . However, if %[ > IORq" Vk then a bifurcation is 
possible and two convective states (one primarily geostrophic, the other 
primarily magnetostrophic) can coexist . The magnetos trophic solution prevents 
further field anplifi cation, but no truely steady-state can be achieved. 
Poloidal and toroidal fields are conparable, and the bifurcation may be 
related to reversals. Tie ultimate field anplltude is around (8npAi2)'i or 
(4HpV^)^, whichever is larger; but complications arise if % » 1. In Pig.l, 
an atterrpt is made to locate the possible positions of planets and satellites 
on the % - Rq diagram (see figure caption) . 
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ROCK MAGNETISM OF EXTRATERRESTRIAL MATERIALS, D.W. Strangway, Dept, 
of Geology, Univ. of Toronto, Toronto, Ont, Can. M5S lAl 

Access to the Apollo samples has suddenly brought a new focus to rock 
magnetism. In most terrestrial environments we are dealing with the iron- 
titanium oxides and with some of the sulphides. The principles of the type of 
magnetism carried by these have been extensively studied. The lunar samples 
however have required us to examine the magnetic properties of iron and iron- 
nxckel phases and of troilite in particular. The properties of the iron- 
nickel-cobalt phases were studied extensively in the 1930 's and it has been 
a question of reviewing old literature. 

Pure iron has a Curie point of 780 C and is extremely magnetic. The 
Curie point decreases slightly with increasing nickel content but the phase 
change between 2 and o phase is magnetic while the 0 -phase is not so this 
transition becomes the effective Curie temperature. Thus kamacite ('^ 7% 
nickel) becomes magnetic on cooling through a sluggish phase change. Taenite 
(y 40% nickel) on the other hand is ? -phase but at this composition it is 
magnetic below a few hundred °C and becomes magnetic on cooling. These pure 
ferromagnetic substances have their bulk magnetic properties controlled 
largely by grain size and grain shape. 

Extremely small grains <100 A are superparamagnetic and carry no 
remanence. Above 100 A they behave as single domain grains and hence are 
extremely stable magnetically. At grain sizes of ly or so they become multi- 
domain and hence appear to be magnetically soft. 

The lunar igneous rocks contain mainly grains of kamacite in the few 
micron size range while the soils and breccias contain in addition abundant 
super-paramagnetic iron as a result of steady reduction of iron silicates. 

The magnetic properties of meteorites are dominated by kamacite taenite and 
plessite (intergrowths of the above two) with some contributions from mag- 
netite and possibly pentlandite. 

Troilite is antiferromagnetic at 320 C and so its magnetic effects can 
be detected, but they contribute little to the overall remanent magnetism. 
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PALEOINTENSITY STUDIES ON LUNAR SAMPLES, D.W. Strangway and 
N. Sugiura, Dept of Geology, Univ. of Toronto, Toronto, Ont. Can. MSS lAl. 

In spite of considerable efforts to deter'mine satisfactory values of 
the ancient lunar field, there are still very few acceptable determinations. 
It is therefore not yet possible to draw firm inferences about changes in the 
magnetic field during the history of the moon. 

The dterminations have proved to be exceptionally difficult for several 
reasons . These are : 

i) many lunar samples have large, soft components of magnetization 
that make it difficult to Isolate the stable component; ^ 

ii) On heating to simulate the original cooling at about 500 C in 
spite of studies using vacuum and controlled fugacity environ- 
ments. There is more work in progress on this problem in several 
laboratories . 

iii) The intergrowth textures present in the lunar sample magnetic 
minerals lead to complex interaction which reduce the elegance 
of the Thellier-Thellier method; 

iv) Attempts to use non-heating methods depend on an assumption of 
critical ratio relating thermoremanent magnetization (TRM) to 
anhysteretic magnetization (ARM) . These ratios are known to 
vary considerably especially when Interactions are present. 

v) Shock effects can add several complexities. In the absence of 

a field, simple shock will reduce the magnetization preferentially 
removing the soft components. In the presence of a field, 
simple shock may Increase the magnetization observed. In the 
course of shocking a rock, the heating process will be either 
minor (little thermal effect on NRM) , major (if heated over 
800°C, NRM will be totally reset) or intermediate. In the inter- 
mediate case the samples are reheated in a very inhomogenous 
manner so that some parts of a sample will be magnetically 
reset and others will not. It is not possible to fully reproduce 
this latter case to create an artificial magnetization for 
paleointensity determination. 

We conclude that while a few samples do reproduce ancient lunar fields 
it is not yet possible to fully reconstruct the ancient lunar field. 
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MAGNETIC STUDIES ON THE ALLENDE METEORITE, N. Sugiura and D.W. 
Strangway, Dept o£ Geology, Univ. of Toronto, Toronto, Ont. Can. MSS lAl 

Earlier studies on the Allende meteorite showed that the natural 
remanent magnetism (NRM) behaved like a simple thermoremanent magnetism (TEM) 
up to 130 C. This magnetization appears to have been acquired in a field of 
about 1.0 oersted. At this temperature it was found that the remanent 
magnetization showed a sharp change. This effect could be due to one of 
several effects i) reheating in a different field at some stage to 130 C; 
ii) chemical changes due to the laboratory heating or; iii) same property of 
the material involved. We now know that iron with minor impurities undergoes 
sharp reversible changes in its magnetic properties at this temperature 
probably explaining this difficulty. 

Subsequent experiments in our laboratories have been done using 
individual chondrules separated from the Allende meteorite. The chondrules 
vary from strongly magnetic to moderately magnetic to weakly magnetic 
chondrules are the light-coloured ones with presumably less iron content. In 
general, it appears that the most stable magnetic component of the individual 
chondrules is randomly oriented. This suggests that the chondrules themselves 
were magnetized by cooling before they assembled into the meteorite. It 
appears that some of these were magnetized by cooling in fields on the order 
of 10 oersteds. 

The chondrules then have in additit>n, a softer component of magnetization 
and the direction of this component while showing some scatter shows a common 
direction. This is evidence that the meteorite was magnetized at or after 
the assembly of the meteorite. The matrix material of the meteorite is also 
quite magnetic. At least two stages and perhaps three stages of magnetization 
may be present — a) magnetization of the chondrules before assembly; b) 
magnetization of the whole meteorite during or near assembly and c) 
remagnetization, perhaps not uniformly, by subsequent metamorphism. 

In either case we conclude that we have a record of early magnetic 
fields that predate the formation of Allende. 



EVOLUTION AND THE STATE OF THE CORES OF TESRESTRIAL PLANETS; 

M. Nafi ToksOz and Albert T. Hsui, Department of Earth and Planetary Sciences, 
Massachusetts Institute of Technology, Cambridge, Mass. 02139 

The internal structure and temperature regimes of a planet play an 
important role in generating and sustaining its magnetic field. Whether the 
planetary magnetic fields are due to remanent magnetism or are of magneto- 
hydrodynamic origin in their cores, may be discriminated by their temperature 
profiles. The first mechanism requires that the temperature within an 
appreciable portion of the planet remains below the Curie 
temperature throughout its history. The second mechanism requires the 
existence of a conducting liquid core such that dynamo actions can be 
derived either by internal heat sources or by external interactions. In this 
paper we calculate thermal evolution models of the inner planets and compare 
their present day temperatures and the state of their cores. The Earth is 
used as a reference in these studies. 

The Moon : The existence or lack of a lunar core is not known. If an iron- 
rich core does exist, it is unlikely that it will exceed 400 km in radius. 

This upper bound on the lunar core size has been established by density 
constraints, seismic studies^ and electro-magnetic induction studies^. Even 
if the lunar core has a size equivalent to the maximum allowed, it only 
represents about 1% of the total volume. Therefore, such a small core will 
not play any significant role in the thermal history of the Moon. Based on 
our thermal evolution studies, if the Moon has a homogeneous accretion origin, 
its core has probably been molten since its formation about 3.5 billion 
years ago^. Whether this small core is able to Induce the remanent magnetism 
observed in the lunar crust is another matter. If the surface magnetism of 
the Moon is the remanence of a primordial field^, this requires a very 
special thermal evolution model. 

Mercury : Mercury has the second largest magnetic field and the highest iron 

content among the inner planets. Whether this field is generated by the 
memory of a primordial field or by a presently active dynamo is still a 
subject of debate. Most data favor a differentiated model of the planet. 
Thermal evolution models that can produce differentiation favor an initially 
hot planet. An iron-rich core can be separated early in the evolutionary 
history of the planet. In order to have a molten or partially molten core 
at the present time, it is necessary to have a small amount of heating within 
the core. 

Mars : The existence of a weak Martian magnetic field has been reported® 

based on the data obtained by the USSR Mars 2 and 3 spacecrafts. The 
accurate determination of the moment of inertia factor^ and density models® 
favor the existence of an Fe-rich Martian core. The radius of the core varies 
between 1500 and 2000 km depending on the composition of the core material. 

If the Martian core is composed of Fe, the formation of the core will take 
place relatively late in its evolution (approximately 2.5 billion years after 
its formation). However, it will remain molten throughout its history. On 
the other hand, if the Martian core is composed of FeS, core formation will 
take place about 1 billion years after formation. Again, it remains molten 
even at the present time. 

Venus : Venus may actually have a weak magnetic field implied by the 

re-evaluation of the magnetic observations of the planet®. Because of the 
scarcity of data, the internal structure of Venus is not known. However, 
if Venus was formed in the present orbit, it probably has an evolution 
history and structure very similar to the Earth's because of the extreme 
similarity in size and mean density of the two planets. Based on this 
assumption, Venus is modelled to have an Fe-rich core similar in size to the 



STATE OF THE CORES 


M.N. ToksBz and A.T. Hsui 

Earth’s. Our thermal models indicate that the outer core of Venus is molten 
at present. The planet may have a solid inner core. 

A comparison of the present-day models of the internal structure 
of the inner planets is shown in Figure 1. It should be noted that these are 
neither the unique nor the definite models of the interior of the inner 
planets. They only represent models favored by the available data and by 
earth comparisons . There appears to be some pattern between the intensity 
of planetary magnetism and the size and state of the cores of the Moon, 

Mars, Mercury and the Earth. Venus stands out as an exception, indicating 
the importance of other factors and the complexity of the origin of 
planetary magnetic fields . 
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MAGNETIC PROPERTIES OF FeNi ALLOYS - MICROSTRUCTUKE AND REMANENT 
MAGNETIZATION MECHANISMS. 

P. J. Wasllewskij NASA/ Goddard Space Flight Center, Astrochemistry Branch, 
Laboratory for Extraterrestrial Physics, Greenbelt, MD 20771. 

Remanent magnetization in lunar samples and meteorites is carried by 
FeNi alloys. The lack of understanding of the NSM in these samples is 
directly attributable to the absence of an experimental data base which would 
provide basic magnetization data in terms of Ni content and thermal history. 
This need in terrestrial rocks has been recognized early by Nagata and 
colleagues and 25 years later, active research on synthetic FeTi oxides 
continues. 

Results of the first systemmatic study, of weak field remanence acquis- 
ition, hysteresis properties, AF and thermal demagnetization, cryogenic 
cycling in zero field of room temperature acquired remanence, and magnetic 
viscosity, of the complete FeNi system are presented. The magnetic results 
are correlated with recognizable micro structure which is composition depen- 
dent and related to thermal history. Metallographic recognition criteria 
are presented which derive from optical and electron optical observations. 
Each mechanism of magnetization, active during a given thermal or dynamic 
history has a discrete imprint which is easily identified. These are based 
on sound physical metallurgy principles. The use of various carry over 
terms, from oxide experience such as SRM (shock remanent magnetization), 

CRM (chemical remanence) and pTRM( partial thermo remanence), in describing 
lunar and meteoritic NRM are incorrect and can lead to serious errors in 
evaluating the NRM. These objections will be explained in terms of the 
experimental results. 

Specimens of FeNi alloys were acquired as spheres which were solidified 
in free fall using techniques developed for the NASA-zero-g solidification 
program. Therefore macroscopic demagnetization factors were constant and 
the serious problem of shape anisotropy was eliminated by the use of spheres. 
The FeNi system is divided into three groups: (a) the bcc FeNi alloys 
(Ni ^ 28%) , (b) the fee FeNi alloys (30% ^ Ni ^ 50%) , and (c) the fee 
ordering alloys (50% ^ Ni ^ 90%) . These groups are easily dealt with in 
terms of mode of occurrence in natural samples, thermal history and shock 
response, and magnetic behavior. Magnetic hardness is related to thermal 
and shock histories and is determined by microstructure induced. If an alloy 
has experienced a specific history which is unknown, but the event is 
discrete, then examination of the microstructure describes what mechanism 
of magnetization is operative and therefore specifies the analysis technique. 
This is possible because unlike oxides the alloy structures are indicative 
of thermal history, dynamic history, and subsequent thermal history. 

For example, in the bcc FeNi alloys we show that the deformation twin 
density reflects the shock level but beyond the pressure transit ion- twinning 
is absent and a special shock structure results (this indicates a demagneti- 
zation has taken place) and as the peak impact levels increase the adiabatic 
decompression produces shock heating and ’self’ anneal is possible. At the 
highest levels we have a case of a quench from ^ 1000°C. All of these shock 
levels have been studied and the magnetic effects are different for each. 
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Indications are that the 'paleointensity’ of a given metal sample 
in a .given test field is- a function of. its- prior ^history.* For example, 
even though the a e ^ a-j^ transformation demagnetizes’ then remagnetizes 
a sample on a microsecond time scale in a given external field the NRM 
acquired ^may be''l-2 orders < o.f- magnitude less -intense.r than that acquired, 
by simple cooling. The ^mis leading ;^p.oint ‘about the; remanence .is that -it> . 
appears stable, since the shock is a hardening process. Specific natural 
and laboratory examples, of -thermal and j'shocfc effects *;covering the' entire 
FeNi .system are .presented'. 

■ The ord.ering fee alloys' are discussed, in. s.ome detail as they are found 
in the or,dered state in serpentinites and meteorites; Since thermomagnet ic 
and slip induced anisotropy, both induced effects, are dominant anisotropies 
in these alloys, this research is of geophysical importance particularly 
since ordered and disordered states have such vastly different magnetic 
propertxes. 

Sources of magnetic hardness in FeNi alloys are discussed and examples 
presented, e. g., twinning, shock structures, duplex a + y structure, 
quench of shock heated structures. Results of anal.ysis' of remanent 
magnetization and hysteresis properties of chondrules from Allende, Chainpur, 
Bjurbole, Alleghan, Ochansk and Pultusk are presented. The results indicate 
that metal size, shape, and distribution are different in different classes 
of chondrite meteorites. Results of thermal demagnetization in zero field 
of the NRM in iron meteorites and chondrites and the FeNi alloys are 
presented to demonstrate the behavior of the NRM in samples with different 
thermal and shock histories. 

This review outlines the previous experimental research on FeNi alloys, 
presents new calibration for the FeNi alloys, and indicates how to 
analyze samples which contain FeNi alloys. 
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SHOCK INDUCED MAGNETIC EFFECTS IN FINE PARTICLE IRON. 

Peter Wasllewski, NASA/ Goddard Space Flight Center, Astrochemistry Branch, 
Laboratory for Extraterrestrial Physics, Greenbelt, MD 20771. 

Since fine particles not observable by optical microscopy are difficult 
to characterize, since particle shape is extremely difficult to handle 
theoretically and experimentally, and since preparation of dispersed 
synthetic systems are near impossible we have adopted a model system to handle 
these problems. The system, is based on the Copper (iron) precipltabion alloy 
Properly handled a narrow particle size distribution with a reasonably 
well defined interparticle spacing can be obtained. The iron precipitates 
are initially spherical in shape as long as d ^ 1200 S. and maintain the 
antiferromagnetic fee high temperature state when cooled to room temperature. 
Therefore, we have a metastable non-magnetic array of spherical iron particles 
waiting for the plastic deformation in the shock pulse to transform them to 
the ferromagnetic bcc state. Since the fee bcc transformation is a shock 
induced first order magnetic transition, an argument by analogy would suggest 
that the hep bcc transition is similar, (i.e., the remagnetization step 
in the bcc ->■ hep bcc' transitions when bcc Fe is shocked at levels 
> P^ (130 Kb) . Therefore we can study magnetic effects during rapid magnetic 
transitions in controlled external fields. 

The model system enabled us to establish the criterion for particle size 
and shape discreteness, i.e., when the magnetized ensemble is brought from 
saturation through E = 0 to the field [h| = (the remanent coercive force) 
and then reversed to pass through H = 0 to an eq.ual and opposite field 
1 h| = the magnetization hits on the upper portion of the major loop. 

This effect is not observed in coarse mD materials, Lodestones, packed 
1 pm 16304 powder, coarse grained basalts j iron meteorites, chondrites, and 
Disko Island and Mt. Fuji basalts containing Fe. This test appears to be 
a necessary condition for SD behavior in rocks or for dispersions where 
magnetization reversal is ideally similar from particle to particle. In fact, 
for |h 1 in excess of Hp, say Hp, a return through H = 0 to H = Hp should put 
M on the upper porition of the major loop and so on for H2 etc. 

In the spherical Fe samples this holds until d ^ 1200 S and true MD 
behavior is described, i.e., the path of M is somewhere within the upper and 
lower parts of the major loop. With electron microscopy we have observed 
substructure in particles from 250 S. to 1200 S. Maximum magnetic stability 
is observed in the 400 S. to 600 particles and 800 S. to 1100 S particles, 
those 200 - 350 S particles are least stable to AF demagnetization. The 
minimum in the Rg (Rp = Hr/Hc) vs 6 (angle) curve is at 90 , i.e., parallel 
to the shock direction with the maximum perpendicular to the shock direction. 

In the model system particle size, shape and interparticle spacing are 
controlled by varying the conditions of the precipitation anneal. The amount 
of transformed fee iron, of a given size increases as the pressure increases 
and for a given pressure the amount transformed depends on the particle size. 
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Wasxlewski, P. 

The maximum anisotropy for a given pressure is observed in the 400-600 R 
particle. When transformed, the sign of the axial remanent vector always 
obeys the sign of the external field indicating that on the microsecond 
time scale of the transformation the magnetic systems sees the external 
field. The amount of remanent magnetization is considerably less than 
would be expected for TRM in the same material. The ratio of remanence 
after the shock to the saturation remanence is < 0.002. 

In the second generation experiments, we will evaluate directly the 
properties of shock induced remanence compared to thermo remanence, and 
evaluate the effects of recrystallization. These combined data should 
provide us with the experimental basis for evaluating shock effects in 
lunar samples and meteorites. 



MODEL MAGNETOSPHERE OF MERCURY BASED ON MARINER 10 DATA 

Y.C. Whang, Catholic University of America, Washington, D.C. 2006^} 

From the first and the third encounter of Mariner 10 with Mercury's mag- 
netosphere, a combined set of thirty minutes of directly observed magnetospher- 
ic magnetic field data was obtained. Direct observations also provided data 
on some boundary conditions of the magnetosphere such as the size of the mag- 
netos-phere and the magnitude of the magnetic field at the stagnation point of 
the magnetosphere. This information is invaluable and unique in studying the 
planetary magnetic field of Mercury. The analysis and interpretation of mag- 
netic field data from Mariner 10 have been performed in two chases. Phase One 
Analysis used certain selected subsets of the magnetic field data, and usually 
disregarded boundary conditions of the magnetosphere. These boundary condi- 
tions are satisfied in Phase Two Analyses of Mercury’s magnetic field. The 
confinement of planetary magnetic field lines inside a magnetosphere- 1 ike 
boundary requires a proper description of the magnetospheric external field. 
Prior to Mercury III encounter, the entire set of Mercury I data was used by 
Whang and Ness to study the model magnetosphere of Mercury using an image di- 
pole representation for the external field. The result was rather surprising 
in that the inclusion of a magnetosphere- 1 i ke boundary revised the estimated 
planetary dipole moment of Phase One Analyses by a factor of 2. After Mer- 
cury III encounter, further studies of Mercury’s model magnetoshpere using the 
entire combined set of Mercury 1 and MI data revealed that the planetary 
field is very much distorted from a simple centered dipole, the planetary 
quadrupole and octupole are quite significant in magnitude compared with its 
dipole moment as shown in Table 1, The presence of a quadrupole moment repre- 


Table 1. Planetary Magnetic Field Calculated 
From Image Dipole Models 


Planetary Field 

Data Base 

Result 

Centered Dipole 
(Whang S Ness, 

Mercury 1 

D=2.3xl0^^G-cm^ 

March 1975) 



Dipole S Quadru- 

Mercury 1 & 1 1 1 

D=2.4x 10 G-cm^ 

pole S Octupole 
(Whang, March 
1977) 

< 

D:Q;0=1 :0.4:0.3 


sents an axial offset of the apparent planetary dipole center by a distance 

d'x-0.5 ((i/D)x planetary radius 

northward from the planetary gravitational center. The presence of an octupole 
moment indicates that the average radius of the source current System inside 
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Whang, Y.C. 


the planet is a significant fraction of the planetary radius. A spherical 
harmonfc representation for the external field has also been used to study 
the model magnetosphere of Mercury. Table 2 summarizes various Phase Two 
models Based on the entire combined set of Mercury I and III data, and their 
calculated deviations between model field and observed field. 


Table 2. Models Based on All Mercury I £ III Data 


Planetary Field 

Resul t 

Deviation 

/ill 

“V A +*^2 

Dipole & Quadru- 
pole £ Octupole 
(Whang, March 1977) 

D=2.4xlO^^G-cm^ 
D:D:0 =1 :0.4:0.3 

a =11. 4y 

<?^=10.3y 

0^=13.07 

20y 

Dipole £ Quadrupole 
(Jackson £ Beard, 
July 1977) 

D=2.'44xl0^^G-cm^ 

D:a=l :0.66 

a =12. Oy 

a^= g.ify 
0^=1 8. Oy 

24y 

Offset Dipole 
(Ng £ Beard 1978) 

D=2.37xl0^^G-cm^ 
0ffset=0.24 Rm 

CT =13. ly 

a^= 8.6y 

0^=19. 2y 
z 

25y 


Whang's quantitative three-dimensional model of Mercury's magnetosphere 
can Identify the various regions of the magnetosphere, for example, the region 
of closed field lines, field lines connected with the tail current sheet, and 
field lines connected with polar regions. The trajectories of Mercury I and 
Mercury III at different regions exhibit very different plasma characteristics 
Of particular interest is that the observation of high-energy charged parti- 
cles are found to be along the field lines connected with the inner edge of 
the tail current sheet. 

In Mercury's magnetosphere, the external field represents a significant 
fraction of the total field. Along the first half of Mercury I ‘s trajectory, 
the total field was observed to increase from 40 to 98 gammas; whereas, the 
external field estimated by Whang's model is in the order on 30 gammas. Thus 
an accurate model representative of the external field becomes very important 
in data interpretation. The location and strength of the two current sources 
for the external field, the magnetopause and the tail current sheet, vary all 
the time responding to conditions of the solar wind and the magnetosheath. 

This temporal variation causes the fluctuation of the external field in the 
magnetosphere particularly in the region where the gradient of the external 
field is large. Making use of the model magnetosphere, the magnitude for the' 
gradient of the external field is calculated on the trajectory plane of 
Mercury I and Mercury III. The area of large gradient positively correlates 
with the large fluctuations in magnetic field observed on Mercury 1. 
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THE LUNAR SURPAOE MAGNETIC SURVEY ALONG THE ROUTE OP 

LUNOKHOD- 2 

Ye. G. Yeroshenko, IZMIRAN, Moscow, USSR. 

The magnetic field measurements on the southern area of 
the Bay Le Monnier along the route of Lxm.okhod-2 were performed 
in Pehruar - April 1973 [l] • Lunokhod-2 vehicle was not spe- 
cially adapted to the magnetic measurements. The significant 
interferences were generated by the moving details of the vehic- 
le, by the electromagnetic deviation and by the thermocurrents 
as well. As a result of the carefull processing of all the body 
of data it suceeded to eliminate the deviation effects at least 
for two horizontal components. The measurement accuracy of the 
horizontal component was 3+6 gammas and of the vertical one - 
-10 gammas. The relative measurements were made with the reso- 
lution of 1 gamma. The contribution of external fields was taken 
into account as well under analysis of the lunar surface fields. 

The magnetic field intensity is on the average 25+30 gam- 
mas on the observed area of the Bay Le Monnier. The horizontal 
component is directed generally to the south-west. The transi- 
tion region between the floor and the slope of the Le Monnier 
crater is characterized by the abrupt field gradients and the 
both field components are markedly decreased at the joint of the 
floor with the crater rim [2]. It is supposed this anomaly is 
caused by the edge effect of the steep contact of the magnetized 
bazalt thickness with the non-magnet ic (or slightly magnetized) 
rocks of the crater rim. The edge effect was distinctly obser- 
ved at the edges of the tectonic break Straight Rille as well 
(fig. 1). 

The magnetic anomalies which were observed when the Liuao- 
khod crossed the small craters (with sizes of 100+450 m) have 
typically maximum or minimum of the vertical component over the 
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crater center and maximum of the horizontal component over its 
edges (fig. 2), It is revealed the trend to increasing of the 



anomaly amplitude with the increasing of the crater diameter as 
well. That kind of magnetic field features may be explained by 
the shock-explosion origin of the crater anomalies. 

The magnetic anomalies caused by the edge effects were 
used for the estimation of the base rocks magnetization provided 
the uniform magnetization of the selenological structures [2]. 
The calculations of the edge effects of the ledges on oblique 
magnetization have indicated that the theoretical curves were 
in good agreement with the observed data if the bazalt magneti- 

_C O A 

zation was of 10 ^ G*cm-^*g and the magnetization vector was 
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directed down at the angle 30+45° to horizon. 

References 

1. Sh.Sh. Dolginov et al*, In "The Soviet -American Conference 

on Cosmochemistry of the Moon and Planets" Ed. by J.H. Pomeroy 
and N.J. Hubbard, HASA, 1977, p« 433* 

2. Sh.Sh. Dolginov et al.. The Moon, v.15, 1976, p. 3. 



IHE GECmGNETIC SECULAR VABIftTIONS IN HISTORIC TIMES 


T. Yukutake, Earthquake Research Instxtate, University of 
Tbkya, Bunl<^-ku, Tokyo, 113 Japan 

In order to maintain a dipole field by a self exciting dynamo process 
a proper magnitude of non-axisymmetric motions with non-zero helicxty are 
reqiiired in addition to the axisymmetric flows predominant in the liquid core. 
Although the concept of helicity has been originally exploited in the 
turbulent dynamo theory, it has been recently revealed that fluid motions 
with sufficient helicity are capable of regenerating the original axi- 
synitietric field irrespective of the scale of motions. At present no 
preference seems theoretically attributable to either type of fluid motions, 
turbulence or a large scale motion, to maintain the geomagnetic dipole field. 
The type of the fluid motions might possibly differ according to the planeta- 
ry bodies. This paper attempts to provide basic data, by reviewing the 
gecmagnetic secular variations, for investigating what kind of fluid motion 
is most plausible within the earth's core. 

The asyitmetry of the fluid motion acccmpanies the asyrmetric magnetic 
field through induction process within the axisymmetric field. If the fluid 
motion is of turbulent -type, rapidly changing fields with small scales could 
be observed in the gecmagnetic field. Mien the dipole field is subtracted 
from the observed fields, there ronain irregular fields, called non-dipole 
field. If the irregular fields change their patterns in a very shiort time, 
the fluid motions in the core may also be regarded as changing their forms 
as quickly. In fact, there was a time vhen the non-dipole field was 
considered to decay in a hundred years or so and therefore it was regarded 
as the product of turbiulent flews near the surface of the core. 

During the past tenyears, spherical liantonic analyses iiave been conducted 
by several different groups for the historical data going back to the 17th 
oentury. Reconstruction of tiie non-dipole charts from these analyses disclosed 
that the life time of the non-dipole field was much longer than had ever been 
considered. The North American continent viiich is now covered with an intense 
positive anomaly in the vertical e x s tp onent was already under the positive 
anomaly in 1600 AD. This suggests that the American anomaly has life time 
longer than 400 years. Examination of the non-dipole charts indicates two 
distinct types of non-dipole anomalies. One is the migratory anomaly and 
the otter is the anomaly staying almost at the same location. A negative 
anomaly existing off the xsest coast of Africa belongs to the first type. 

It was located m the Indian Ocean in the 17th century. Since then it has 
drifted westwards steadily with an average velocity of 0.28degrees/year and 
reached the present location. On the other hand, an intense positive ancmaly 
that now covers the Asian continent with its. center around Mongolia is the 
second type ancmaly. Until 1650, no ancmaly was recognizable in Mongolia. 

In 1700 a small ancmaly was created and continued to grow to the largest 
non-dipole ancmaly in the same place. An average rate of increase in the 
intensity aimounted to 50 nT/year. The rapid growth, of the ancmaly oould be 
interpreted by a development of localized fluid flows in the core. However, 
this phenoitenon is understandable in a ccmpletely different way based on 
steady field models. 



SBCOIAR VARIATIONS IN HISTORIC TIMES 


Yiikutake, T. 


By making use of the spherical harmonic coefficients, it is possible to 
decorpose the asyimetric parts of the geomagnetic field into two types of 
field, the drifting field and the standing field, with the constant intensity 
during the period of analyses. The decomposition was made of the coefficients 
since the 17 th century up to degree 4. Reconstruction of the geomagnetic 
fields from these two types of the fields approximates the observed secular 
variations fairly well. It was demonstrated that not only the steady drift 
of the African anomaly but also the rapid growth of the Mongolian anomaly 
could be v«ll accounted for ly a simple superposition of the drifting and 
the standing fields. These may be an indirect proof of the long persistence 
of the taro types of the fields, at least over several hundred years. 

One of the differences noticeable betaroen the taro types of fields is 
lack of the sectorial term in the standing field for the degree taro. i.e. 
n=2 and itf=2. The aitplitude of the standing field is only about 700 nT, \drlle 
the drifting field has the anplitude of about 2000 nT. A similar feature is 
also seen in the sectorial term of degree 3, The standing omponent of n=4 
and itt=4 is neither large. The question vdiether or not the standing field 
generally lacks tte sectorial terras except for ilie equatorial dipole 
coiponent, which has the largest anplitude among the standing carponents, 
will be solved by examination of higher terms in future. Though the data 
available are insufficient yet, the renarkable lack of the sectorial terms 
for the degrees of 2 and 3 in the standing field might provide an iirportant 
clue to the investigation of the physical prcx:ess of creating the taro kinds 
of fields. 

The several centuries length is not sufficient to identify \^ether 
"the standing field" is strictly standing at the same locality or it is also 
subjected to a slow drift that may not be detected for several centuries 
observation. Regarding the equatorial coiponent, however, observation 
evidently deviates from the drifting and the standing model even within 
several hundred years period. It was shown recently that introduction of 
a ccitponent drifting eastwards, Instead of a standing field, better 
appraximated the observed secular variations of the equatorial dipole. 

Axisymmetric fields expressed by zonal terms of the spherical harmonics 
are known to hae undergone considerable fluctuations in the past. On the 
other hand, there is no explicit data yet suggesting the variations in the 
intensity of the drifting and the standing fields. If they have changed 
their intensities, the observed secular changes in the hantionic ccmponents 
are suposed to deviate from the model cciiputed on the steady field assunption. 
The deviations around the coiputed for seme of the cenponents amount to 
about 20 to 30 percent of the anplitude of the drifting field. This seotis 
to place an upper bound on the amplitude variation of the drifting field. 

Examination of the geomagnetic secular change during the past several 
hundred years indicated that the asymmetric fields that can be expressed by 
lower spherical harmonics had been stable. This might be favourable for 
considering the existence of considerably stable fluid motions with laterally 
large scale in the core. 
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